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This  report  consists  of  two  volusies: 

Volume  I  -  Generalised  Hake  Geometry  and  Hake  Effect  on  / 

Rotor  Airloads  and  Performance 

Volume  II  -  Hake  Geometry  Charts 


SDIMARY 


An  analytical  investigation  is  described  which  was  conducted  to  gener¬ 
alize  the  wake  geometry  of  a  helicopter  rotor  in  forward  flight  and  to  demon¬ 
strate  the  influence  of  including  wake  deformation  in  the  prediction  of  rotor 
airloads  and  performance. 

A  first  level  generalized  wake  model  was  developed  for  a  helicopter  rotor 
operating  in  steady  level  forward  flight  based  on  theoretically  predicted  wake 
geometries.  The  generalized  wake  model  consists  of  wake  geometry  equations  in 
which  tip  vortex  distortions  are  generalized  as  displacements  from  the 
classical  undiatorted  tip  vortex  geometry.  Tip  vortex  coordinates  are  gener¬ 
alised  in  terms  of  tip  vortex  age,  blade  azimuth,  rotor  advance  ratio,  thrust 
coefficient,  number  of  blades,  and  rotor  disk  attitude  for  a  selected  repre¬ 
sentative  blade  design.  Emphasis  was  placed  on  the  region  of  the  wake  near 
the  rotor  and  the  predominant  axial  tip  vortex  distortions  normal  to  the  rotor 
disk.  An  exploratory  study  of  the  sensitivity  of  tip  vortex  geometry  to  rotor 
operating  condition  and  design  parameters  was  conducted  to  guide  the  gener¬ 
alized  wake  formulation.  The  generalized  wake  equations  were  programmed  in  a 
computer  module  which  can  be  used  to  expeditiously  generate  wake  coordinates 
for  use  with  computer  programs  for  rotor  airflow  and  airloads  prediction.  A 
set  of  wake  geometry  charts  are  presented  in  Volume  II  to  provide  a  con¬ 
venient,  readily  accessible  source  for  approximating  rotor  forward  flight  wake 
geometry  and  identifying  wake  boundaries  and  locations  of  blade-vortex 
passage. 

Tip  vortex  geometry  and  the  sensitivity  of  predicted  rotor  airloads  and 
related  rotor  performance^and  blade  bending  moments  to  various  rotor  inflow 
and  wake  geometry  models  are  demonstrated  herein  for  selected  low  to  moderate 
high  speed  flight  conditions  of  a  representative  rotor  and  the  H-34  rotor. 

The  inflow/wake  models  consist  of  uniform  rotor  inflow  and  variable  inflow 
based  on  classical  undistorted  wake  geoswtry,  predicted  distorted  wake 
geometry,  and  the  newly  developed  generalized  wake  geoswtry.  Relative  to  the 
undistorted  wake  model,  the  forward  and  lateral  sides  of  the  wrte  are 
distorted  toward  the  rotor  in  forward  flight.  This  results  in  the  prediction 
of  a  greater  nueber  of  close  blade-vortex  passages  and  larger  local  gradients 
in  blade  airloads  with  the  distorted  wdw  models.  The  influence  of  the  wake 
on  blade  airloads  is  related  to  the  degree  of  proximity  of  the  tip  vortices  to 
the  rotor  and  the  nui^r  of  close  blade-vortex  passages.  Tor  steady  level 
flight,  this  influence  is  predicted  to  generally  increase  with  decreasing 
advance  ratio,  decreasing  rotor  disk  nose  down  tilt,  and  increasing  nuaber  of 
blades.  Increasing  thrust  coefficient  increases  the  undistorted  wake  deflec¬ 
tion  angle  away  from  the  rotor  but  siaatltaneously  increases  the  amplitude  of 


the  distortions  of  the  tip  vortices  toward  the  rotor.  Tip  vortex  distortions 
relative  to  the  undistorted  tip  vortex  geowetry  are  predicted  to  be  insensi¬ 
tive  to  rotor  disk  attitude. 

Although  additional  experiaental  data  are  required,  liwited  data  frooi 
previous  investigations  tend  to  support  these  theoretical  observations.  The 
occurrence  of  close  blade-vortex  passages  and  the  high  sensitivity  of 
predicted  blade  airloads  to  small  deviations  of  blade  to  tip  vortex  distance 
dearanstrate  the  requirement  for  an  improved  blade-vortex  interaction  model. 
Although  the  need  for  further  analytical  refinement  is  evident,  the  degree  to 
which  distorted  wake  models  (generalised  or  predicted)  improved  the  correla¬ 
tion  of  the  measured  B-34  airload  and  bending  moment  distributions  is 
encouraging  and  indicates  that,  with  future  emphasis  on  blade-vortex  interac¬ 
tion  modeling,  distorted  wake  smthodology  has  the  potential  to  provide  a  more 
accurate  tool  for  blade  airload  prediction. 
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INTRODUCTION 


One  of  the  most  fundamental  factors  affecting  the  performance  of  rotor- 
craft  systems  is  the  induced  influence  of  the  main  rotor  wake  on  rotor 
airloads  and  related  performance.  The  taak  of  predicting  the  flow  field 
induced  by  a  helicopter  rotor  continues  to  be  of  primary  importance  for 
providing  and  evaluating  improved  rotor  deaigna.  The  continuing  development 
of  more  accurate  aerodynamic  methods  for  predicting  the  temporal  and  spatial 
variation  of  the  induced  flow  velocities  at  the  r  tor  blades  and  in  the  rotor 
flow  field  will  result  in  significant  advancea  in  rotor  aerodynamic 
methodology  important  to  nearly  every  area  of  rotary  wing  design.  These  areas 
include  performance,  vibrations,  stability,  structural  reliability,  and 
acoustics.  Improved  capability  to  predict  the  detailed  rotor  flow  character¬ 
istics  will  facilitate  the  prediction  of  more  accurate  blade  airloads  and  a 
more  rational  design  of  the  blade  atructure.  This,  in  turn,  can  lead  to 
reduced  aircraft  development  and  maintenance  costs  along  with  improved 
performance  because  of  Che  more  optimum  design.  Also,  the  capability  to 
predict  rotor-induced  flow  velociciea  away  from  the  rotor  is  important  for  the 
calculation  of  aerodynamic  interference  effects  at  aircraft  components  such  as 
Che  fuselage,  wing,  tail  rotor,  and  empennage  surfaces. 

The  accurate  prediction  of  helicopter  rotor  blade  airloading  has  been  a 
continuing  goal  in  Che  helicopter  industry.  Significant  progress  has  been 
made  as  advanced  computational  methods  accounting  for  wake,  aeroelastic,  and 
fuselage  interference  effects  have  been  developed.  However,  the  helicopter 
designer  cannot  yet  rely  on  airloads  methods  to  provide  him  with  sufficient 
accuracy  to  determine  performance,  loads,  vibrations,  stability,  noise,  and 
interference  effects  over  Che  range  of  flight  conditions.  Of  the  numerous 
phenomena  influencing  rotor  airloads  that  require  further  research,  a  primary 
one  is  wake  interaction.  Fur  many  flight  regimes  the  tip  vortices  from  Che 
main  rotor  blades  pass  sufficiently  close  Co,  or  even  impact,  following  blades 
such  as  to  produce  a  complex  flow  environment  which  influences  both  the  local 
blade  loading  and  Che  vortex  characteriatics. 

Classically,  the  wake  shed  from  a  rotor  blade  in  forward  flight  has  been 
modeled  as  a  rigid,  undeforawd  helical  sheet  of  vorticity  of  varying  strength, 
with  a  trajectory  defined  by  rotor  attitude  and  flight  condition.  Generally, 
these  vortex  sheets  are  then  asideled  by  assuming  that  their  influence  can  be 
approximated  by  uaing  diacrete  vortex  filaments  of  appropriate  atrength  posi¬ 
tioned  within  the  plane  of  the  vortex  sheet  (Fig.  1).  These  wake  models  are 
then  used  to  calculate  the  induced  inflow  at  the  rotor  blade  by  applying  the 
Biot-Savart  law  to  the  particular  rotor-wake  geometry. 
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In  recent  years,  Che  demand  for  improved  rotorcraft  performance  and 
airload  prediction  has  led  to  Che  developsMnt  of  srare  detailed  wake  models.  A 
major  characteristic  of  Che  wake  geosKtry,  currently  ignored  in  most  forward 
flight  analyses  in  general  use  in  industry,  is  the  deformation  of  the  wake 
from  Che  rigid  helical  shape.  In  particular,  the  rollup  of  the  sheet  into  a 
deformed  tip  vortex  whose  position  in  space  is  not  consistent  with  the 
classical  geometry  is  shown  in  Figs.  2  and  3.  Any  methodology  «dtich  does  not 
consider  this  deformation  phenomena  may  be  in  question  for  certain  applica¬ 
tions. 

The  major  technological  advance  in  recent  years  to  rotor  induced  velocity 
methodology  has  been  Che  improvement  in  the  mathematical  modeling  of  the  rotor 
wake.  This  has  becosie  possible  with  the  advent  of  high-speed  digital 
computers.  Many  efforts  have  been  directed  toward  providing  this  methodology 
as  sho%m  in  the  survey  presented  in  Ref.  1  and  other  examples  listed  in  the 
References.  As  a  major  part  of  a  long-term  effort  to  advance  the  technology 
associated  with  the  aerodynamic  and  aeroelastic  behavior  of  helicopter  rotor 
blades,  the  United  Technologies  Research  Center,  UTRC,  has  been  conducting 
studies  on  Che  influence  of  the  rotor  wake  on  rotor  airflow,  airloads,  and 
dependent  rotorcraft  characteristics.  Some  of  Che  UTRC  research  programs 
conducted  in  this  area  are  discussed  in  Refs.  1  through  10.  The  scope  of  this 
effort  has  included  both  analytical  and  experimental  investigation  programs, 
and  the  development  of  distorted  and  undisCorted  wake  prediction  analyses  as 
applied  Co  single-  and  multi-rotor  configurations  and  compound  helicopters 
operating  under  hovering  and  forward  flight  conditions.  Methods  developed 
have  been  provided  to  several  organisations  within  the  Covernment  and 
industry.  The  UTRC  Prescribed  Wake  Rotor  Inflow  Analysis,  the  Rotorcraft  Wake 
Analysis,  and  the  Rotor  Wake  Geometry  Analysis  are  described  in  Refs.  1  Co  10 
and  brief  descriptions  will  be  presented  in  a  later  section  of  this  report. 

The  functions  and  features  of  Che  Rotorcraft  Wake  Analysis  are  shown  in  Fig. 

4.  A  major  feature  of  these  analyses  is  the  provision  for  various  wake 
modeling  options. 

In  hovering  flight,  the  need  to  accurately  predict  the  position  of  the 
rotor  wake  is  well  documented  (Refs.  1  to  13).  The  use  of  lifting  line  or 
surface  models  with  discrete  trailing  vortex  systems  modeling  the  wake  of  each 
blade  is  state-of-the-art  in  the  industry  for  general  hovar  performance 
predictions.  Accurate  hover  performance  analyses  use  prescribed  wake  models, 
whose  tip  vortex  and  inboard  sheet  geometry  have  been  determined  from 
empirical  studies.  Originally  the  vtke  geosMtry  in  hover  was  predicted 
numerically,  but  this  was  found  to  be  of  insufficient  accuracy  and  computa¬ 
tionally  costly,  and  had  to  ha  calculated  for  each  flight  condition.  This 
method  was  replaced  by  empirical  mathoda  describing  the  wake  geometry  from 
flow  visualisation  data  which  was  generalised  at  UTRC  for  a  wide  range  of 
parameters  (Refs.  9  and  10).  This  generalisation  wae  performed  for  a  range 


of  hovering  rotor  geometries  and  flight  conditions  Co  obtain  analytical 
expressions  defining  Che  hovering  wake  geomeCry  in  a  compact  form.  Further 
expansion  of  the  empirical  wake  geometry  general ixat ion  for  hover  was  reported 
in  Ref.  11.  A  validated  lifting  line  hover  performance  analysis  is  the 
Circulation  Coupled  Hover  Analysis  Program  (CCHAP)  developed  at  Sikorsky 
Aircraft  (Ref.  8).  This  analysis  is  an  extension  of  Che  Ref.  9  prescribed 
wake  hover  analysis.  CCHAP  adjusts  the  generalized  hover  wake  geometry  by 
coupling  the  wake  geometry  and  blade  airloading  in  a  consistent  manner.  Other 
hover  method  developments  have  included  a  lifting  surface  representation  of 
the  blade  (e.g..  Refs.  11  Co  13).  The  importance  of  accurately  modeling  the 
tip  vortex  geometry  in  forward  flight  is  not  as  extensively  documented  as  for 
hover,  but  the  need  has  been  recognized  by  several  investigators  (e.g.,  Refs. 
4,  14  to  17). 

A  sample  computer  plot  of  the  deformed  tip  vortices  for  a  two-bladed 
rotor  at  a  30  kt  flight  condition,  as  calculated  by  the  UTRC  Wake  Geometry 
Analysis,  is  presented  in  Fig.  3.  The  predicted  wake  geometry  for  this  condi¬ 
tion  is  compared  with  that  of  a  hover  condition  and  the  predicted  distorted 
and  undistorted  tip  vortex  patterns  of  a  15  kt  condition  in  Fig.  5,  which  is 
reproduced  from  Ref.  6.  The  tip  vortex  geometry  predicted  by  the  same 
analysis  for  one  blade  of  a  two-bladed  rotor  at  89  kt  is  shown  Fig.  6, 
which  is  reproduced  from  Ref.  5.  In  several  UTRC  studies  suitaarized  in  Ref. 

6,  experimental  and  theoretical  techniques  were  applied  to  investigate  rotor 
wake  and  airflow  characteristics  for  low  speed  flight.  Flow  velocities  were 
measured  with  a  laser  velocimeter  and  wake  features  were  determined  from  flow 
visualizat  on  techniques.  For  example,  measured  wake  boundaries  and  tip 
vortex  locations  from  Ref.  6  are  presented  in  Figs.  7  and  8  which  show  that 
the  path  of  the  tip  vortices  goes  above  and  then  down  through  the  rotor  at  low 
speed,  and  blade-vortex  impact  occurs.  The  application  of  laser  velocimetry 
to  locate  the  tip  vortex  passages  under  the  advancing  blade  is  discussed  in 
Ref.  18.  A  comparison  of  calculated  and  measured  airloading  for  that  test 
program  is  presented  in  Ref.  19.  Photographs  of  a  model  rotor  wake  in  a 
water  tunnel  are  presented  in  Ref.  20.  In  the  Ref.  5  study,  the  water  tunnel 
data  of  Ref.  20  and  the  UTRC  Wake  Geometry  Analysis  were  utilized  to  explore 
the  forward  flight  wake  problem  and  the  potential  of  forward  flight  wake 
generalisation  was  identified.  For  forward  flight,  experimental  wake  data  is 
extreswly  limited,  and  theoretical  wake  predictions  have  been  limited  to 
isolated  cases.  A  systematic  determination  of  wake  geometry  and  airload 
sensitivity  to  wake  deformations  over  a  range  of  operating  conditions  was 
required.  The  theoretical  study,  reported  herein,  systematically  explored  the 
characteristics  and  influence  of  the  forward  flight  wake. 

A  primary  reason  for  the  limited  use  of  distorted  wake  geometries  in 
rotor  analyses  is  the  computer  cost  and  complexity  involved  in  wake  geometry 
computations.  Hake  coordinates  are  calculated  by  large  computer  progrema. 


such  as  the  UTRC  Rotor  Wake  Geometry  Analysis  described  in  Ref.  4,  in  which 
the  computer  time  required  is  typically  many  times  the  combined  requirement 
for  the  calculation  of  all  of  the  other  rotor  aerodynamic  and  performance 
characteristics. 

Applications  of  the  Rotor  Wake  Geometry  Analysis  at  UTRC  indicated  that 
the  potential  existed  for  generalizing  rotor  tip  vortex  geosietry  for  forward 
flight  conditions.  From  information  generated  during  the  Ref.  5  study,  the 
development  of  techniques  for  generalizing  forward  flight  tip  vortex  geometry 
as  functions  of  primary  parameters  appeared  feasible.  It  was  recognized  that 
wake  generalization  would  facilitate  the  development  of  a  set  of  generalized 
equations,  charts  and/or  tables  «rtiich  would  define  rotor  wake  geometry  in  a 
form  suitable  for  use  in  computer  analyses  for  predicting  rotor  airloads  and 
related  performance,  structural,  and  dynamic  response  characteristics.  Such 
generalizations  could  be  prograssned  in  stodular  form  to  make  wake  geometries 
imsiediately  available  in  such  rotor  analyses. 

Considering  the  above,  this  investigation  was  conducted  with  the 
following  objectives: 

1.  To  analytically  demonstrate  the  significance  of  rotor  tip  vortex 
distortions  on  the  aerodynamic  environment  of  a  helicopter  rotor  in  forward 
flight  with  particular  emphasis  on  rotor  airload  distribution  and  associated 
rotor  performance. 

2.  To  generalize  rotor  tip  vortex  coordinates  in  terms  of  primary  vari¬ 
ables  in  equation  form  for  incorporation  in  a  computer  subroutine  to  provide  a 
first  level  wake  geometry  generalization  for  the  calculation  of  rotor  inflow 
and  airload  distributions  required  for  helicopter  rotor  performance  and  blade 
aeroelastic  response  programs. 

3.  To  perform  an  initial  assessment  of  the  developed  generalized  wake 
representation  and  its  application  for  analytical  rotor  flight  simulation. 

4.  To  provide  a  set  of  graphical  wake  geometry  charts  for  undistorted 
and  distorted  tip  vortex  coordinates,  blade-vortex  passages,  and  wake 
boundaries  for  a  comprehensive  range  of  primary  rotor  wake  parameters:  number 
of  blades  (2  and  4),  advance  ratio,  thrust  coefficient,  rotor  attitude,  blade 
azimuth  and  vortex  age. 

The  investigation  is  described  in  the  following  sections  of  this  report 
which  are  entitled  (1)  Technical  Approach,  (2)  Description  of  Analyses,  (3) 
Selection  of  Analytical  Parameters  and  Procedures,  (4)  Generalized  Wake 
Geometry,  and  (5)  Influence  of  Wake  Geometry  on  Rotor  Airloads  and  Perfor¬ 
mance.  A  set  of  wake  geometry  charts  is  provided  in  Volume  II  of  this 
report . 


TECHNICAL  APPROACH 


In  accordance  with  the  objectives  stated  in  the  Introduction,  a  technical 
program  was  performed  «diich  was  comprised  of  the  following  tasks: 

1.  Determination  of  the  significance  of  wake  distortions  on  rotor 
airloads  and  related  rotor  performance  and  blade  bending  moments. 

2.  Hake  generalisation. 

3.  Development  of  a  generalised  wake  computer  module. 

4.  Assessment  of  the  wake  generalisation. 

5.  Provision  of  a  set  of  rotor  wake  geoaietry  charts. 

The  technical  approach  used  in  the  performance  of  each  of  these  tasks  is 
summarised  below. 


Significance  of  Hake  Distortions 

The  computer  program  for  the  UTRC  Rotor  Hake  Geometry  Analysis  was  used 
in  conjunction  with  the  UTRC  Prescribed  Hake  Rotor  Inflow  Analysis  (Single 
Rotor  Version  of  the  UTRC  Rotorcraft  Hake  Analysis),  and  the  Sikorsky  General¬ 
ised  Rotor  Performance  Analysis,  as  shown  in  Fig.  9,  to  determine  the  signif¬ 
icance  of  main  rotor  tip  vortex  distortions  on  the  aerodynamic  environment  and 
performance  of  a  helicopter  rotor  in  steady  forward  flight.  The  wake 
geometry,  blade  loading  distribution,  and  rotor  performance  characteristics 
were  calculated  for  varying  combinations  of  rotor  design  and  flight  condition, 
and  the  results  based  on  a  distorted  wake  geometry  were  compared  with  the 
results  based  on  both  an  undistorted  wake  geometry  and  uniform  inflow.  The 
calculated  wake  distortions  consisted  of  self-induced  distortions  of  the  tip 
vortices  of  an  isolated  rotor,  and  the  distorted  wake  geometry  used  in  the 
inflow  and  performance  calculations  included  these  tip  distortions  along  with 
an  undistorted  inboard  wake  representation.  The  time-varying  circulation 
strength  of  the  rolled-up  tip  vortices  of  the  distorted  wake  was  based  on  the 
time-varying  peak  of  the  blade  radial  loading  distribution  as  calculated  in 
the  Prescribed  Hake  Rotor  Inflow  Analysis.  Flight  speeds  of  approximately  20- 
130  knots  (0.05-0.3  advance  ratio)  were  investigated  to  emphasise  the  low  to 
moderate  speed  range  where  wake  effects  are  most  influential.  The  effects  of 
tip  vortex  distortions  on  airloads  and  overall  performance  were  determined  as 
a  function  of  rotor  radius,  asimuth,  and  flight  condition  for  several  flight 
conditions  and  design  variations  for  a  representative  rotor.  In  addition. 
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three  B-34  rotor  operating  condition  coad»inationa  Mere  selected  on  the  basis 
of  available  test  data,  and  for  the  above  co^inations  calculated  airloads, 
controls,  and  bending  sioaients  were  coapared  with  the  test  data.  An  available 
Sikorsky  normal  modes  aeroelastic  rotor  analysis  was  used  to  predict  bending 
moments. 


Wake  Generalisation 

Using  results  from  the  UTRC  Rotor  Wake  Geometry  Analysis  as  used  in  con¬ 
junction  with  the  UTRC  Prescribed  Wake  Rotor  Inflow  Analysis  and  the  Sikorsky 
Generalised  Rotor  Performance  Analysis,  techniques  were  developed  for  general¬ 
ising  rotor  tip  vortex  coordinates  in  terms  of  primary  parameters.  The  wake 
generalisation  findings  reported  in  Ref.  5  were  expanded  and  investigated 
analytically  for  additional  rotor  designs  and  flight  conditions  to  establish  a 
data  base  from  tfhich  to  generalise  the  tip  vortex  distortions,  develop  gener¬ 
alisation  techniques,  and  determine  the  range  of  values  of  the  variables  over 
which  the  generalisation  techniques  apply.  The  computer  programs  were  used  to 
compute  distorted  and  undistorted  wake  geometries  for  several  combinations  of 
significant  parameters  such  as  advance  ratio,  number  of  blades,  thrust  coeffi¬ 
cient,  rotor  disk  attitude,  blade  twist,  and  blade  aspect  ratio.  A  first 
level  wake  generalisation  was  successfully  developed.  Emphasis  in  the  gener¬ 
alisation  was  placed  on  rotor  advance  ratio  (0.05  to  0.3),  number  of  blades  (2 
and  4),  thrust  level  and  rotor  disk  attitude. 

It  is  noted  that  it  was  not  the  intent  of  this  investigation  to  neces¬ 
sarily  develop  a  wake  generalisation  ^ich  is  all  inclusive  for  every  param¬ 
eter  and  the  complete  range  of  each  parameter.  Recognising  that  this  wss  the 
first  systematic  attempt  to  generalise  forward  flight  wake  geometry  and  the 
prohibitive  costly  computer  requirements  to  establish  an  all  inclusive  data 
base  in  this  initial  effort,  the  intent  was  to  develop  a  first  level  general¬ 
isation.  The  scope  of  this  generalisation  was  determined  by  the  cga^Axity 
and  range  of  applicability  of  the  generalised  techniques  that  were  developed, 
and  the  predefined  computer  utilisation  limits  of  the  project. 

In  order  to  reduce  the  large  number  of  values  required  to  document  the 
coordinates  of  the  time-varying  tip  vortex  geometries  for  a  specific  flight 
condition,  application  of  several  generalisation  techniques  were  explored. 
These  included  the  developawnt  and  investigation  of  the  adequacy  of  numerical 
curve  fit  equation  techniques  and  the  generalisation  of  coordinates  for 
varying  rotor  positions  and  wake  revolutions.  Also,  the  adequacy  of  undis¬ 
torted  wake  coordinates  for  certain  wake  regions  was  investigated.  The  wake 
geometries  were  generalised  in  equation  form  with  tabulated  coefficients  for  a 
representative  blade  design,  which  is  dependent  on  the  number  of  blades, 
advance  ratio,  thrust  coefficient,  and  disk  attitude. 


DevelopBcnc  of  a  Generalised  Wake  CoapuCer  Module 

The  results  of  the  wake  generalisation  were  prograaaed  in  a  computer 
subroutine  and  incorporated  in  the  UTRC  Prescribed  Wake  Rotor  Inflow  Program. 
The  computer  subroutine  was  also  provided  as  a  module  to  the  computer  program 
for  the  Rotorcraft  Wake  Analysis  (Ref.  2).  Wake  coordinate  equations  and 
tabulations  with  interpolation  procedures  were  coded  in  the  subroutine  to 
generate  the  wake  distortions  corresponding  to  a  desired  rotor  design  and 
operating  condition.  The  resulting  generalised  wake  module  coupled  with  the 
Prescribed  Wake  Rotor  Inflow  Program  or  the  Rotorcraft  Wake  Analysis  provides 
the  capability  for  calculating  distorted  wake  induced  velocities  in  a 
computationally  efficient  sianner. 


Assessment  of  the  Wake  Generalisation 

The  generalised  wake  module  was  used  in  conjunction  with  the  UTRC 
Prescribed  Wake  Rotor  Inflow  Program  and  the  Sikorsky  Generalised  Rotor 
Performance  Program  to  perform  an  initial  assessment  of  the  generalised  wake. 
Comparisons  of  wake  coordinates,  airloads,  and  associated  performance  results 
based  on  use  of  the  generalised  wake  module  were  compared  with  the  following: 
(1)  corresponding  coordinates,  airloads  and  performance  based  on  distorted 
wake  geometries  obtained  directly  from  the  Rotor  Wake  Geometry  Analysis,  (2) 
corresponding  coordinates,  airloads,  and  performance  based  on  undistorted  wake 
geometries,  and  (3)  airloads  and  performance  based  on  constant  rotor  inflow. 
Comparisons  were  made  for  a  representative  rotor  and  for  three  flight  condi¬ 
tions  of  an  H-34  full-scale  rotor  for  which  test  data  were  available. 


Rotor  Wake  Geometry  Charts 

Wake  geometry  charts  were  developed  using  computer  graphics  techniques. 
These  charts  provide  the  necessary  information  to  determine  the  relative  posi¬ 
tion  of  the  tip  vortices  trailed  from  helicopter  rotor  blades  for  a  range  of 
parameters  sufficient  for  most  conventional  steady  forward  flight  operating 
conditions.  Both  the  classical  undistorted  and  generalised  distorted  wake 
models  are  included.  The  charts  can  be  used  for  a  variety  of  applications 
which  require  the  geometric  relationship  between  the  tip  vortices  and  spatial 
locations  relative  to  the  helicopter.  In  particular,  locations  of  blade- 
vortex  interactions  and  wake  boundaries  beneath  the  rotor  can  be  rapidly 
estimated  using  these  charts.  A  series  of  tip  vortex  patterns  showing  top, 
side,  rear,  end  isometric  views  are  also  included.  The  wake  geometry  charts 
wera  prepared  as  Volume  II  of  this  report. 
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Th«  abov*  technical  acCivitiaa  and  tha  raaulta  tharafroai  are  described  in 
detail  in  the  following  aactiona  of  thia  report.  Brief  daacriptioos  of  the 
coaputer  analyaaa  uaad  in  thia  invent igat ion  are  praaented  in  the  next 
section. 


DESCRIPTION  OP  ANALYSES 


Several  coaputer  analyeee  were  weed  during  tbia  inveatigation  to  predict 
rotor  wake  geoawtry,  inflow,  eirloeda,  perforiMnce,  and  bending  aowenta. 
Included  in  tbia  aection  of  tbe  report  ere  brief  deacriptiona  of  tbeae 
enelyaea  along  with  a  deacription  of  tbe  aolution  procedurea  required  for 
tbeir  uae. . 


Tbe  Preacribed  Wake  Rotor  Inflow  Prograa 

Tbe  UTRC  Preacribed  Wake  Rotor  Inflow  Progra  (F389SR)  ia  a  aingle  rotor 
veraion  of  tbe  Rotorcraft  Wake  Analyaia  (RWA)  deacribed  in  Ref.  2.  A  chart  of 
tbe  operational  atatua  and  capabilitiea  of  RHA  ia  abown  in  Fig,  4.  Tbeae 
analyaea  are  uaed  in  conjunction  with  other  analyaea  to  provide  either  inflow 
velocitiea  at  the  rotor  diak  for  rotor  perforaance  and  blade  reaponae  predic- 
tiona  or  blade  circulation  diatributions  for  uae  in  the  Rotor  Wake  Geoawtry 
Analyaia  (Ref.  4).  A  flow  chart  of  tbia  type  of  coupling  to  fon  a  coaplete 
"rotor  analyaia"  ia  preaented  in  Pig.  9.  The  detaila  of  tbe  aingle  rotor 
veraion  of  the  inflow  progran  uaed  in  tbia  inveatigation  are  preaented  in  Ref. 
21.  A  brief  technical  deacription  for  tbia  aingle  rotor  veraion  followa. 

Technical  Deacription 

The  UTRC  Preacribed  Wake  Rotor  Inflow  Prograai  (hereafter  referred  to  aa 
the  rotor  inflow  prograa  or  inflow  analyaia)  conputea  tbe  bound  circulation 
and  induced  inflow  dietributiona  along  the  rotor  bladea  baaed  on  a  preacribed 
wake  aodel  and  aodeling  of  tbe  aerodynaaica  at  each  blade  aection.  Aa  abown 
in  Fig.  9,  tbe  aection  operating  conditiona  are  generally  preacribed  froa 
blade  aotion  and  control  inforaation  obtained  froa  a  aeparate  blade  airload 
and  reaponae  prograa.  The  wake  aodel  ia  either  generated  internally  in  the 
coaputer  prograa  (undiatorted  w^cc)  or  preacribed  froa  a  aeparate  prograa. 

Briefly,  tbe  aatheaatical  aodel  in  tbe  rotor  inflow  prograa  conaiata  of 
the  repreaentation  of  each  blade  by  a  aegaented  lifting  line,  and  the  helical 
wake  of  tbe  rotor  by  diacrete  aegaented  vortex  filaaenta  conaiating  of 
trailing  vorticity  which  reault  froa  tbe  apanwiae  variation  of  bound  circula¬ 
tion.  Figure  10  illuatratea  the  coordinate  ayataa  uaed  in  tbia  analyaia  for 
the  rotor  and  wake  geoaetry.  The  circulation  of  the  wake  for  each  blade 
changea  with  aaiauth  poaition  and  ia  periodic  for  each  rotor  revolution.  The 
bladea  are  divided  into  a  finite  nubber  of  radial  aegaenta,  and  tha  induced 
velocity  at  tbe  center  of  each  aelected  blade  aegaent  ia  coaputad  by  auaaing 
the  contributiona  of  each  bound  and  trailing  wake  aegaent  (Fig.  11).  Tbe 
contribution  of  each  vortex  aegaent  ia  obtained  through  uae  of  the  Biot-Savart 
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equation  trttich  expresses  Che  induced  velocity  in  teras  of  the  circulation 
strength  of  the  vortex  segment  and  its  geoaetric  position  relative  to  Che 
blade  segment  at  which  Che  induced  velocity  is  desired.  The  bound  circulation 
distribution  is  determined  by  relating  Che  wake  circulations  to  Che  bound 
circulations,  expressing  the  wake  induced  velocities  in  terms  of  the  unknown 
bound  vortex  strengths  by  weans  of  the  Biot-Savart  equation,  and  developing  a 
set  of  simultaneous  equations  in  which  each  tera  relates  the  bound  circulation 
and  local  blade  angle  of  attack  at  each  blade  segment  and  axiauth  position. 
These  equations  thus  involve  Che  known  flight  condition,  wake  geometry,  and 
blade  motion  and  control  parameters  and  Che  unknown  bound  circulation  values. 
Solution  of  these  equations,  to  include  measured  airfoil  section  character¬ 
istics,  yields  the  desired  bound  circulation  values,  which,  when  combined  with 
Che  appropriate  geometric  relations  in  the  Biot-Savart  law,  produce  the 
required  rotor  induced  velocity  distribution.  Many  of  the  fundamental  assump¬ 
tions  are  based  on  those  of  Ref.  22. 

The  rotor  inflow  program  is  based  on  the  following  formulation  and 
assumpt ions; 

1.  Each  blade  is  represented  by  a  lifting  line  located  at  the  blade 
quarter  chord  line  (bound  vortex)  and  divided  into  a  finite  number  of  segments 
(blade  segments)  each  having  a  different  circulation  strength  (Fig.  11).  The 
aerodynamic  characteristics  at  Che  centers  of  each  segment  are  assumed  to  be 
representative  of  the  entire  segment  over  a  finite  asimuCh  interval.  The  use 
of  lifting  line  rather  chan  lifting  surface  theory  permits  the  use  of 
empirical  airfoil  data  for  nonlinear  aerodynamics  such  as  stall. 

2.  The  wake  is  represented  by  a  finite  number  of  vortex  filaments 
trailing  from  the  blade  segment  boundaries.  Each  filament  is  divided  into 
straight  segments,  the  lengths  of  which  are  determined  by  a  specified  wake 
asimuCh  interval  which  is  equivalent  Co  a  specified  blade  azimuth  interval 
(see  Fig.  11).  The  circulation  strength  of  each  trailing  vortex  segment  is 
constant  along  its  length,  and  is  equivalent  to  Che  difference  between  Che 
circulation  values  of  its  adjacent  bound  vortex  segments  in  accordance  with 
the  Helmholtz  laws  of  conservation  of  vorticiCy.  The  circulation  strengths  of 
different  vortex  segments  along  a  vortex  filament  vary  in  accordance  with  the 
variation  of  Che  bound  circulations  with  azimuth  position. 

3.  Viscous  dissipation  effects  on  the  wake  circulation  strengths  are 
neglected  in  that  the  circulation  of  a  given  wake  segment  is  constant  with 
tisw.  However,  the  number  of  wake  revolutions  retained  in  the  analysis  can  be 
limited  to  evaluate  an  abrupt  dissipation  of  the  wake.  Also,  the  viscous 
rollup  of  the  tip  vortex  can  be  approximated  by  combining  tip  filaments  beyond 
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■  prescribed  rollup  azimuth  interval.  In  addition,  a  vortex  core  size  is 
assigned  to  each  vortex  filament.  Potential  theory  is  assusted  to  apply 
outside  the  vortex  core.  Inside  the  vortex  core  zero  flow  or  rigid  body 
rotation  is  assumed.  Other  vortex  core  models  can  be  readily  incorporated. 

4.  It  is  assumed  that  the  rotor  is  operating  in  steady-state  flight. 

The  inflow  and  wake  from  each  blade  is  assumed  to  be  periodic  with  blade 
spacing.  That  is,  the  inflow  and  wake  geometry  are  the  same  for  each  blade 
when  at  a  specific  azimuth  position. 

5.  Shed  wake  segments  (segments  normal  to  the  trailing  wake  segments 
mentioned  above)  arising  from  the  tisw  rather  than  radial  variations  of  bound 
vorticity,  are  not  included  in  the  wake  model.  Although  this  omission 
technically  violates  the  Helmholtz  law,  it  is  believed  that  a  more  accurate 
representation  of  the  shed  wake  effects  is  obtained  through  the  use  of  experi¬ 
mental  unsteady  airfoil  characteristics  in  the  analysis.  This  implies  that 
the  primary  effects  of  the  shed  vorticity  are  those  associated  with  the  shed 
wake  immediately  behind  the  blade,  and  thus  can  be  approximated  by  those  of  a 
fixed-wing  type  of  wake.  Miller  (Ref.  23)  shows  that  this  is  generally 
reasonable.  This  approach  not  only  permits  a  factor  of  two  reduction  in 
computer  time  and  prevents  unrealistic  results  associated  with  a  discrete 
filament  wake  model,  but  more  importantly,  permits  nonlinear  unsteady  stall 
effects  to  be  included  in  an  empirical  manner  in  the  airloads  prediction 
program. 

6.  The  wake  geometry  is  prescribed  from  analytical  or  experimental 
results.  Various  options  are  available  for  selecting  the  representation  of  the 
wake  model. 

7.  The  airfoil  at  the  blade  is  assumed  two-dimensional  (radial  velocity 
components  are  neglected).  For  the  linearized  circulation  solutions  and 
associated  calculations,  a  set  of  lift  curve  slopes,  stall  angles  of  attack, 
and  angles  of  zero  lift  are  provided  which  vary  with  Mach  number  and  blade 
thickness.  These  may  be  based  on  unsteady  aerodynamic  data,  if  available. 

Below  stall  the  lift  curve  slope  is  assumed  constant.  The  blade  section 
circulation  is  limited  to  a  constant  value  for  angles  of  attack  above  stall. 

An  optional  nonlinear  solution  procedure  can  be  used  which  employs  measured 
tabulated  airfoil  data  to  account  for  nonlinear  aerodynamics  including  stall. 

8.  In  the  blade-wake  geometrty  calculations  the  blades  are  assumed 
rigid.  However,  flexibility  effects  may  be  included  in  the  circulation  solu¬ 
tion  by  providing  the  necessary  noninduced  velocity  at  each  segment  associated 
with  flexible  blade  motions  and  Controls  from  a  blade  response  program.  Also, 
an  anhedral  or  swept  tip  may  be  prescribed. 


9.  T«at*»ti«l  *04  radial  iadoead  aalaeity  caapoaaata  aia  aaglaetad  aad 
aull  aa(la  aaaaaytioaa  arc  iacludad  ia  tba  liaaac  circalatiaa  aalatiaa. 

Thaca  aaaaapciaaa  «ra  cliaiaacad  ia  the  aytiaacl  aoa-liaear  aalatiaa. 

10.  The  aaradyaaaic  iatarfereacc  cffaeta  af  the  rotor  hvA,  £ueelevc»  etc. 
ere  mtmmllj  a^lected.  Bowever,  thaea  can  he  eccoaewdatad  ia  the  eolutioa, 
if  deairad,  ae  Mverievoeed  flow  velocity  ^reaeate.  The  rotor  ia  aaewed  to 
he  aperatiag  oat-of-groaad  effect. 

Wehe  CeeewtrT 

The  rotor  iaf  low  yrotrae  requirea  thtf  the  r<^or  wake  geoaetry  he 
apecified  ia  order  for  circulationa  aad  iadaced  velocitiea  to  be  dMeraiaed. 
Generality  regarding  the  apecificatioa  of  the  wAa  geaaetry  ia  provided  ia  the 
coapater  prograa.  Thia  ia  accoapliched  by  allawing  the  wake  caordiaatea  for 
the  wake  aegaeat  end  poiatc  to  be  input  fraa  aa  eateraal  aaurca  (tape,  carda, 
etc.)  for  coaputatioaal  purpoaea,  if  deairad.  At  the  eaaet  ot  thia  iaveatiga-' 
tioa  there  were  two  prograa  ^tiona  avail^le  far  wake  geaaetry:  (1)  a 
claaaical  uadiatorted  w^e  aodel  (Fig.  1),  available  iateraally,  aad  (2)  an 
analytical  w^  aodel  (diatorted  wake.  Fig.  3),  obtaiaad  froa  a  aeparate  wake 
geaaetry  prograa.  Xa  addition,  a  general iaed  wake  aodula  developed  during  the 
eourae  of  thia  iaveatigatioa  haa  aince  been  iacluded. 

Generally  the  claaaical  uadiatorted  wake  aodel  haa  been  weed  for  perfor- 
aance  prediction.  It  haa  the  leant  coaplex  wake  geoawtry  that  ia  aia^ly  a 
function  of  the  flight  eoadition  aad  aoaeatua  inflow  velocity.  The  coord- 
iaatea  of  thia  helical  wake.  Which  ia  deewad  ia  forward  flight,  are  generated 
iaternally  in  the  prograa.  A  brief  deacriptioa  of  the  claaaical  aad  diatorted 
wake  geoawtriea  ia  provided  below. 


£i^ij|tOT^ed  1 

la  ita  aii^leat  fom,  the  watsa  can  be  aaawawd  to  be  a  claaaical  ua¬ 
diatorted  rtMwed  helical  abcet  of  vortieity  defined  froa  aoaaatw  eaaaidera- 
tioaa  (hereafter  referred  to  aa  the  claaaical  wake  or  undiatorted  wake).  A 
aaaple  undiatertod  wake  repreaeatation  ia  ahewe  in  Pig.  1.  The  ceordiaatea  of 
the  aadiatorted  wwha  repreaeatatioh  we  ^ninad  ftton  the  rotor  advaaee  ntie, 
w,  thrwot  oeefficioat,  G^,  Ohd  ooglo  of  nit4^  relative  to  the  tig  path 
ploae,  dm*  For  cnaaglo,  the  top  view  of  the  tip  vortan  filoaenta  ia 
obtai^  directly  froa  Ao  kelicOidal  path  ^  the  btndd  tip  on  It  traaolotea 
at  the  voloeity  foooa^  aad  retotoa  at  the  velocity  fli.  Wha  oide  view  ia 
depoBdaat  on  the.  <ipdui.  .iiKow  aaglo  .(eaiMA!  ..batwaoa;  ooraal  to  the  tetor  diA  aad 
wake  bowadary),  ^ieh  ia  aanially  dofinod  fraa  aaaaetwa  coaaidwotiaao. 
Iqaotieao  1^  chi  nadiotertod  w0ka  oeerdiaotM  will  he  pceaoatai  ia  a  later 
oaetiea  ei  the  report. 
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Dia^t  or£e^  W^e_Jlode^l 

A  distorted  wake  aodel  way  be  incorporated  in  the  analysis  via  tape 
input.  A  rotor  wake  geoaetry  analysis,  such  as  that  developed  at  UTRC  and 
described  in  the  following  section,  can  be  used  to  generate  wake  coordinates 
on  tape  for  use  in  the  rotor  inflow  analysis.  Figure  3  is  an  example  of  such 
a  distorted  wake  aodel. 


The  Rotor  Wake  Geoawtry  Analysis 

The  UTRC  Rotor  Wake  Geoaetry  Analysis  predicts  the  distorted  tip  vortex 
geoaetry  of  a  rotor  in  steady  forward  flight.  A  detailed  description  of  the 
development  and  fundamental  concepts  used  in  this  analysis  along  with  sample 
applications  are  presented  in  Ref.  4.  Because  this  analysis  is  used  to 
predict  the  distorted  tip  vortex  geoaetry,  and  it  is  these  predicted 
geometries  which  determine  the  influence  of  wake  distortions  on  the  predicted 
rotor  performance  for  this  investigation,  a  description  of  the  technical 
approach  is  presented  which  provides  more  detail  than  the  descriptions  of  the 
other  analyses  used  in  the  investigation.  The  basic  technical  approach  and 
numerical  procedure  used  in  this  analysis  follow,  and  the  major  assumptions 
and  approximations  used  in  the  analysis  are  noted. 

Technical  Description 

The  wake  model  used  in  the  calculation  of  the  wake  distortions  is 
composed  of  trailing  and  bound  vortex  filaments.  The  trailing  vortex  fila¬ 
ments  are  divided  into  finite  length  straight  vortex  elements  whose  end  points 
are  free  to  convert  at  velocities  equal  to  the  vector  sum  of  the  free  stream 
velocity  and  the  velocities  induced  by  the  trailing  and  bound  vorticity.  The 
lengths  of  the  wake  eleswnts  are  related  to  the  aximuth  increment  chosen.  The 
wake  is  truncated  at  a  specified  number  of  revolutions  downstream  of  the 
blade.  Generally,  the  numerical  accuracy  of  the  analysis  is  increased  by 
decreasing  the  aximuth  increment  and  increasing  the  number  of  wake  revolu¬ 
tions.  The  exact  choice  of  these  parameters  is  dependent  upon  the  flight 
condition  under  study  and  the  trade-off  between  required  accuracy  and  accept¬ 
able  computing  cost.  Details  of  this  selection  of  parameters  is  noted  in 
another  section  of  this  report.  The  accuracy/computing  cost  trade-off  also 
influenced  the  wake  interaction  model  and  the  selection  of  interacting 
elements  as  described  below. 

The  tip  vortex  elements  distort  over  successive  time  increments  from  an 
initially  prescribed  geometry  (generally  the  clessicel  undistorted  wake 
geometry)  under  the  influence  of  the  induced  flow  distribution  velocities 
contributed  by  each  of  the  vortex  elements  of  the  blade  bound  vortex,  the 
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inboard  trailing  wake  filaments,  and  the  tip  vortex  filaments.  A  regionalisa- 
tion  procedure,  which  distinguishes  between  "near"  and  "far"  elements  to 
reduce  computer  cost,  described  in  detail  in  Ref.  4,  will  be  described  briefly 
at  the  end  of  this  report  section. 

The  bound  vorticity  of  each  blade  is  represented  by  a  single  segmented 
vortex  filament  located  at  the  quarter  chord  (lifting  line).  In  the  general 
case,  the  strength  of  this  vorticity  varies  with  both  radial  position  and 
azimuth  angle.  Up  to  ten  radial  values  of  bound  vorticity  can  be  accounted 
for  in  the  program  at  each  azimuth  angle. 

The  inboard  wake  elements,  representing  the  inboard  vortex  sheet  of  each 
blade,  are  generally  prescribed  as  undistorted  and  their  influence  at  the  tip 
vortices  is  generally  based  on  the  undistorted  wake  geometry. 

The  wake  elements  generated  by  azimuthal  variations  in  blade  bound 
circulation  have  not  been  included  in  the  wake  model.  This  shed  wake  repre¬ 
sents  one  of  the  unsteady  effects  arising  from  the  azimuthal  variation  of 
blade  bound  circulation.  It  was  reasoned  during  earlier  evaluation  studies  of 
undistorted  wake  variable  inflow  methods  that  the  number  of  discrete  shed  wake 
filaments  that  could  be  accommodated  in  any  practical  analysis  was  so  limited 
that  the  use  of  discrete  filament  shed  wake  model  was  of  questionable 
accuracy,  particularly  in  the  important  region  near  the  blade.  Although  it  is 
recognized,  from  a  rigorous  theoretical  viewpoint,  that  the  Helmholtz  laws  of 
vortex  continuity  have  been  violated  as  far  as  the  wake  geometry  problem  is 
concerned,  the  omission  of  the  shed  wake  effects  on  the  distortions  of  the 
trailing  vortex  filaments  is  assumed  to  have  only  a  secondary  effect.  This 
assumption  is  partially  based  on  the  fact  that  the  circulation  strength  of  the 
shed  wake  is  significantly  less  than  that  of  the  trailing  filaments,  particu¬ 
larly  that  of  the  tip  vortices. 

The  input  circulation  strengths  assigned  to  the  wake  elements  are 
generally  transferred  from  the  Prescribed  Vake  Rotor  Inflow  Program  (Fig.  9). 
In  such  a  calculation,  a  steady  flight  condition  is  assumed  and  the  circula¬ 
tion  of  a  given  wake  element  is  related  to  the  bound  circulation  of  the  blade 
section  at  the  time  the  wake  element  was  created.  Viscous  dissipation  effects 
on  wake  circulation  are  currently  neglected. 

A  provision  is  included  in  the  analysis  for  eliminating  the  irrotation- 
ality  assumption  (i.e.,  applicability  of  the  Biot-Savart  law)  in  computing 
induced  velocities  at  points  located  within  an  arbitrary  distance  from  any 
vortex  element.  This  distance  is  considered  to  be  the  vortex  core  radius. 
Irrotational  potential  flow  induced  velocities  are  calculated  at  points  out¬ 
side  the  vortex  core  with  a  limiting  value  at  the  vortex  core  boundary.  The 
results  presented  herein  were  obtained  assuming  a  vortex  core  diameter  of  1 
percent  of  the  rotor  radius  (approximately  0.1  to  0.15  chord  length).  Due  to 


the  current  limitations  of  the  state-of-the-art  concerning  rotor  vortex  core 
structure,  the  detailed  rotational  induced  effect  for  points  vithin  the  core 
and  changes  in  core  radius  consistent  with  distortions  of  the  lengths  of  the 
elements  were  neglected.  Removal  of  these  assumptions  can  be  easily  accoomo- 
dated  when  more  data  on  the  structure  and  induced  influence  of  the  vortex  core 
for  helicopter  rotor  applications  become  available.  Finally,  interference 
effects  from  sources  external  to  the  rotor  blades  such  as  the  fuselage,  hub, 
and  wings  were  not  included  in  the  analysis  for  this  isolated  rotor  study. 

Numerical  Procedure 

The  procedures  followed  in  the  numerical  computations  of  the  distorted 
tip  vortex  geometry  are  described  below.  The  procedure  involves  two  basic 
steps  —  the  use  of  the  Biot-Savart  law  to  compute  the  wake  distortion 
velocities  produced  by  a  given  wake  geometry  and  the  integration  of  these 
velocities  over  a  small  time  increment  (azimuth  increment)  to  establish  a  new 
tip  vortex  geometry.  These  steps  are  successively  repeated  until  a  converged 
periodic  wake  geometry  results. 

The  computation  is  initiated  by  specifying  the  wake  geometry  and  circula¬ 
tion  to  be  those  corresponding  to  the  classical  undistorted  wake  solution. 

The  three  components  of  distorting  velocities  induced  by  each  wake  eleoient  at 
the  end  points  of  the  trailing  wake  elements  are  then  computed  using  Biot- 
Savart  relationships  (Ref.  4).  The  resulting  changes  in  the  coordinates  of  a 
given  element  are  given  by  the  following  relations,  which  also  include  the 
effect  of  the  flight  velocity,  V,  the  rotor  disk  attitude,  ot^pp*  blade 

tip  speed,  OR. 

Ax  ■  (V£  +  Vcosojpp)  At/R 

Ay  ■  Vj  At/R 

y 

Az  ■  (V£  +  Vsino^p)  At/R 
s 

At  -  RAt/OR 

where  all  distances  are  nondimensionalized  by  the  rotor  radius,  R. 

The  equations  are  baaed  on  the  assumption  that  the  tip  vortex  distortion 
velocities  remain  constant  during  the  time  step.  At.  Also,  the  length  of  any 
vortex  element  is  allowed  to  vary  as  the  wake  distorts  and  the  blades  are 


allowed  to  rotate  to  new  positions.  The  rotation  of  the  blades  leads  to  the 
generation  of  new  wake  elements;  however,  for  convenience  and  computer  cost 
efficiency,  the  total  number  of  wake  elements  retained  in  the  computation  is 
held  fixed  by  discarding  elements  from  the  ends  of  the  trailing  filaments 
farthest  from  the  blades.  Once  the  new  wake  position  has  been  established, 
the  procedure  is  repeated  and  the  process  continues  for  a  sufficient  number  of 
time  steps  to  achieve  convergence  of  the  periodic  distorted  tip  vortex 
geometry. 

The  development  of  efficient  computation  techniques  and  simplifying 
approximations  to  the  basic  technical  approach  has  been  emphasized  at  UTRC  to 
provide  a  practical  computer  program  for  engineering  applications.  Without 
these  simplifications,  the  cost  of  computing  the  distorted  geometry  of  a 
typical  rotor  wake  becomes  prohibitive  for  routine  applications  —  even  with 
the  speed  of  current  computers.  This  is  especially  true  for  computations 
involving  inboard  vortex  filaments,  large  numbers  of  blades,  large  numbers  of 
wake  revolutions,  and  small  azimuth  increments. 

The  number  of  wake  elements  involved  in  the  basic  technical  approach 
outlined  above  is  essentially  equal  to  the  product  of  the  number  of  blades, 
number  of  vortex  filaments  per  blade,  number  of  wake  revolutions,  and  number 
of  wake  elements  per  revolution  (which  is  fixed  by  the  azimuth  increment).  It 
is  therefore  essential  that  simplifying  assumptions  and  efficient  computing 
techniques  be  used. 

To  reduce  the  computational  complexity  of  the  basic  technical  approach 
described  above,  an  approximate  method,  described  in  Ref.  4,  was  developed 
based  on  the  rapid  reduction  of  the  velocity  induced  by  a  wake  element  with 
increasing  distance  from  the  element.  During  the  development  of  this 
analysis,  it  was  initially  considered  that  elimination  of  interactions  between 
wake  elements  spaced  greater  than  a  specified  distance  apart  might  be  a 
reasonable  simplification  to  the  analysis.  However,  it  was  reasoned  that 
although  the  individual  contribution  of  any  one  wake  element  to  the  induced 
velocity  at  a  point  in  the  wake  might  be  insignificant,  the  elements  of  a 
rotor  wake  system  could  be  oriented  in  such  a  way  that  the  integrated  effect 
of  a  group  of  such  elements  could  be  significant.  As  a  result,  the  following 
approximate  approach  was  developed.  The  wake  surrounding  each  tip  vortex 
element  end  point  was  divided  into  a  near  and  far  wake  region.  The  grouping 
of  elements  into  near  and  far  wake  regions  is  determined  by  specifying  a 
limiting  value  on  the  distortion  velocity  induced  by  a  single  wake  element. 

If  the  element's  induced  velocity  is  greater  than  the  limit  value, it  is 
classified  as  a  near  wake  element.  Thus,  the  near  wake  region  contains  the 
elements  which  individually  contribute  significantly  to  the  distortion 
velocity  at  a  wake  point  under  consideration.  The  induced  effect  of  each  of 
these  near  wake  elements  is  therefore  recalculated  for  each  time  step  in  the 
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distortion  iteration.  In  contrast,  the  distortion  velocities  at  the  wake 
point  induced  by  the  far  wake  elements  are  only  computed  once  —  for  the 
initial  time  step.  The  summation  of  these  individual  distortion  velocities 
represents  the  integrated  induced  effect  of  the  far  wake  region  on  the  element 
end  point  under  consideration  and  is  suintained  constant  for  succeeding  time 
steps.  In  effect,  this  assumes  that  the  distortion  velocities  induced  by  the 
elements  in  the  far  wake  region  can  be  based  on  the  initial  undistorted  wake 
model.  The  neglect  of  the  distortions  of  the  far  wake  elements  when  computing 
their  induced  effects  at  a  given  wake  point  is  reasonable  because  the  induced 
velocity  is  insensitive  to  small  changes  in  the  position  of  a  point  relative 
to  the  vortex  when  the  point  and  the  vortex  are  widely  separated.  The  reten¬ 
tion  of  the  far  wake  effect  by  the  approximate  method  described  above  results 
in  a  satisfactory  compromise  between  accuracy  and  computer  cost  requirements. 


The  Generalized  Rotor  Performance  Analysis 

The  Sikorsky  Generalized  Rotor  Performance  (GRP)  Analysis,  a  method  for 
predicting  rotor  performance,  is  described  in  Ref.  24,  and  has  been  refined 
and  improved  since  its  inception.  The  analysis  is  a  rigid,  articulated  blade 
response  analysis  that  is  a  significant  refinement  of  the  classical  Wheatley- 
Bailey  blade  element  method  (Ref.  25)  for  rotor  performance  predictions.  The 
siethod  used  in  the  analysis  represents  a  numerical  stepwise  (in  the  direction 
of  rotation)  solution  in  the  equations  of  motion  of  rigid  blades.  Tabulated 
airfoil  data  are  used.  Details  of  the  numerical  procedure  are  available  in 
the  above  references.  The  latest  version  of  this  analysis  can  be  coupled  with 
the  UTRC  Prescribed  Wake  Inflow  Analysis  to  obtain  compatible  inflow  and  rotor 
performance  predictions. 


The  Normal  Modes  Rotor  Aeroelastic  Analysis 

To  calculate  the  flexible  blade  response  and  bending  moments  for  this 
investigation,  a  version  of  the  Sikorsky  Normal  Modes  Rotor  Aeroelastic 
Analysis  (Y200)  was  used  (Refs.  26  to  29).  The  version  chosen  for  this 
investigation  was  the  T200  program  version  which  had  been  transmitted  to  the 
government  by  Sikorsky  Aircraft  for  use  in  the  Aerodynamical ly  Comformable 
Rotor  Reaearch  Program,  Refs.  28  and  29.  A  brief  discussion  of  the  analysis 
is  presented  below.  A  detailed  description  of  the  equations  of  motion  used  is 
given  in  Ref.  26. 

The  computer  analysis  constitutes  a  numerical  method  for  solving  the 
differential  equations  governing  the  fully  coupled  flapwise,  edgewise,  and 
torsional  responae  of  flexible  rotor  blades.  A  normal  mode  approach  combined 
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with  a  finite  difference,  atep-by-atep  integration  of  the  equationa  ia 
employed.  The  former  ia  widely  uaed  throughout  the  aircraft  induatry  and 
aaaumea  that  the  deflectiona  of  the  rotor  blade  can  be  repreaented  by  auaana- 
tiona  of  the  blade'a  natural  vibratory  modea.  Although  an  elaatic  body  haa  an 
infinite  number  of  auch  modea,  the  normal  mode  approach  owea  ita  aucceaa  to 
the  fact  that  only  a  relatively  amall  number  of  modea  ia  generally  required  to 
adequately  approximate  the  elaatic  deformationa  of  intereat.  The  uae  of 
normal  modea  alao  reaulta  in  certain  aimplificationa  of  the  equationa  of 
motion  (auch  aa  elimination  of  certain  dynamic  coupling  terma)  thereby 
facilitating  the  aolution  of  the  equationa  throu^  numerical  integration 
techniquea  on  a  digital  computer. 

The  veraion  of  the  analyaia  aelected  for  thia  inveatigation  haa  the 
following  capabilitiea  (recent  veraiona  have  increaaed  capabilitiea) : 

1.  The  program  includea  the  effecta  of  up  to  4  flapwiae,  3  chordwiae, 
and  2  toraional  elaatic  modea  aa  well  aa  the  flapping  and  lagging  rigid  body 
modea  for  an  articulated  rotor. 

2.  The  chordwiae  and  toraion  modea  can  be  auppreaaed  to  permit  appli- 
cationa  for  which  only  the  influence  of  the  uncoupled  flapwiae  ax>dea  are  of 
intereat. 

3.  The  computer  program  can  be  applied  equally  well  to  bladea  having 
many  typea  of  root-end  boundary  conditiona  by  aimply  providing  aa  input  the 
appropriate  mode  ahapea  and  natural  frequenciea.  Any  degree  of  hinge 
reatraint  from  fully  articulated  (with  flapping  hinge  offaet)  to  rigid  in 
either  the  flapwiae  or  chordwiae  direction  can  be  atudied  and  both  lag  hinge 
damping  and  control  ayatem  flexibility  can  be  included. 

4.  The  elaatic  deflectiona  of  the  rotor  blade  and  the  aerodynamic  forcea 
generated  by  it  are  fully  coupled.  Aerodynamic  forcea  in  the  blade  equationa 
can  be  conaidered  to  be  compoaed  of  (1)  componenta  independent  of  blade  motion 
(forcing  functiona),  (2)  componenta  dependent  on  blade  deflectiona  (apring 
forcea),  and  (3)  componenta  dependent  on  the  rate  of  change  of  blade  deflec¬ 
tiona  (damping  forcea).  Analyaea  which  do  not  couple  blade  deflectiona  and 
blade  loada,  therefore,  conaider  only  the  firat  of  theae  three  componenta.  At 
high  forward  apeeda  the  contributiona  of  the  coupled  aerodynamic  loada  become 
increaaingly  important. 

5.  The  computer  program  can  be  uaed  to  compute  time  hiatoriea  of  blade 
motion  even  when  auch  motion  ia  unatable.  Such  conditiona  of  inatability  or 
near-inatability  can  exiat  at  high  advance  ratio  flight  conditiona  where 
either  retreating  blade  toraional  divergence  or  blade  flapping  inatability 
impoaea  a  limit  on  maximum  apeed. 


6.  Dynamic  cnupling  batwean  the  flapwise,  chordwiae,  and  toreional 
degrees  of  freedom  have  been  included. 

7.  Airfoil  lift,  drag,  and  pitching  moment  data  at  angles  of  attack  from 
0  to  360  deg  and  at  Mach  numbers  from  0  to  0.9  are  employed. 

8.  High  frequency  blade  response  in  the  primary  modes  can  be  predicted. 
This  can  be  particularly  important  at  high  forward  speeds  and  low  rotor  speeds 
when  the  entire  blade  enters  the  reversed  flow  registe. 


SELECTION  OF  ANALYTICAL  PARAMETERS  AND  PROCEDURES 


Numerical  Prediction  Procedures 

In  this  investigation  the  analyses  were  used  to  predict  rotor  wake 
geometry,  inflow,  airloads,  performance,  blade  response,  and  bending  sioments. 
It  is  therefore  necessary  to  explain  the  calculation  procedure  that  was  used 
to  predict  this  information.  Details  of  the  individual  analyses  used  in  the 
procedure  have  been  outlined  in  the  preceding  section.  The  procedure  used  to 
calculate  this  information  using  combinations  of  these  analyses  is  explained 
below. 

The  procedure  consisted  of  running  the  separate  analyses  in  a  sequential 
and  iterative  manner  to  provide  compatible  solutions.  A  flow  chart  of  this 
procedure  is  shown  in  Fig.  9.  In  this  chart,  either  the  Generalised  Rotor 
Performance  Analysis  (GRP)  or  the  aeroelastic  analysis  (Y200)  is  indicated  as 
the  blade  airloads  and  response  program  to  provide  rotor  controls  and  blade 
siotion  information  for  the  variable  inflow  analysis.  The  airloads  program  was 
first  run  with  a  uniform  inflow  model  using  a  parabolic  tip  loss  factor  to 
provide  initial  input  as  noted  above  for  the  inflow  analysis.  The  inflow 
analysis  was  then  run  with  the  desired  wake  model  (classical  or  distorted)  to 
calculate  the  induced  inflow  and  blade  circulation  distribution.  The  airloads 
progrmn  was  then  rerun  using  the  variable  induced  inflow  as  input  to  calculate 
new  rotor  controls  and  blade  response.  This  procedure  was  repeated  until  the 
resulting  inflow  solution  and  rotor  controls  and  blade  response  were  compat¬ 
ible.  To  obtain  the  distorted  wake  geometry,  the  circulation  distributions 
were  first  used  from  the  classical  wake  solution  as  input  to  the  wake  analysis 
and  the  resulting  wake  geometry  was  calculated.  This  wake  geometry  was  used 
as  input  to  the  inflow  analysis  to  obtain  the  distorted  wake  circulations  and 
induced  inflow,  and  the  iterative  procedure  noted  above  was  used  to  obtain  the 
resulting  rotor  performance  and  blade  response  solution  compatible  with  the 
distorted  wake  inflow.  This  procedure  could  have  been  repeated  by  using  the 
distorted  wake  circulation  distributions  as  input  to  the  wake  analysis  to 
recompute  a  new  wake  geometry  and  the  complete  cycle  repeated,  however  past 
experience  has  shown  that  it  is  generally  not  required  for  high  flight  speed 
conditions  and  can  result  in  numerical  convergence  problems  for  flight  condi¬ 
tions  with  close  blade-vortex  encounters  using  existing  blade-vortex  inter¬ 
action  models. 


Input  Parameter  Selection  for  the  Wake  Geometry  Analysis 


The  UTRC  Wake  Geometry  Analysis  can  be  run  with  several  input  options 
which  can  affect  the  speed  of  computation  and  numerical  accuracy  of  the 
results  for  a  given  condition.  These  are:  (1)  the  number  of  wake  revolu¬ 
tions,  (2)  the  blade  azimuth  increment,  (3)  the  number  of  blade  stations,  (4) 
the  number  of  time  step  intervals,  and  (5)  the  tolerance  for  selecting  near 
and  far  wake  regions  for  the  computations  of  the  self-induced  velocity  field 
(wake  distortion  velocities). 

Because  of  the  significant  cost  of  running  the  Wake  Geometry  Analysis  and 
the  large  number  of  computer  cases  required  for  this  study,  the  need  was 
recognized  for  determination  of  the  appropriate  combination  of  these  param¬ 
eters  considering  cost  effectiveness,  numerical  convergence,  and  accuracy. 

The  number  of  blade  stations  is  basically  determined  by  the  solution  resolu¬ 
tion  required  in  the  rotor  inflow  calculation  program,  so  this  input  was  kept 
consistent  with  the  requirements  of  that  program.  Nine  (9)  stations  were 
used  for  this  investigation.  Knowledge  of  the  number  of  wake  revolutions 
required  for  a  given  flight  condition  had  been  gained  from  previous  experience 
with  the  analysis.  The  number  of  wake  revolutions  used  in  this  investigation 
varied  from  5.5  to  7.  Initial  runs  were  made  with  the  distorted  wake  analysis 
for  different  blade  azimuth  increments,  number  of  time  steps,  and  far  field 
tolerance  criteria.  From  these  runs  it  was  determined  that  a  blade  azimuth 
increment  of  15  degrees  would  be  used  for  both  two-  and  four-bladed  rotors. 
Thus,  throughout  these  investigations,  15  degrees  was  used  as  the  blade 
azimuth  increment  for  computing  the  tip  vortex  distortions.  The  use  of  one 
value  (15  deg)  for  all  cases  also  allowed  for  the  generation  of  a  consistent 
number  of  wake  geometry  data  points  for  each  condition  for  the  wake  general¬ 
ization  portion  of  the  investigation.  The  number  of  time  steps  (incremented 
blade  azimuth  positions)  required  for  distorted  tip  vortex  geometry  conver¬ 
gence  is  dependent  on  the  number  of  blades  used,  with  the  four-bladed  cases 
requiring  more.  The  values  used,  idiich  were  determined  to  give  reasonable 
convergence,  were  48  and  72  time  steps  for  the  two-  and  four-bladed  configura¬ 
tions  respectively.  The  study  of  the  induced  velocity  tolerance  to  determine 
near  and  far  wake  regions  for  each  wake  point  led  to  the  final  value  of  0.5 
feet  per  aecond  as  a  reasonable  value  for  this  rotor.  Smaller  values  made  the 
computation  time  per  time  step  too  large  without  gaining  significant  increases 
in  numerical  convergence  and  speed  of  convergence,  and  larger  values  did  not 
give  reasonable  convergence  or  convergence  rates  although  the  computation  time 
per  time  step  is  reduced.  It  should  be  noted  here  that  the  question  of 
convergence  is  an  engineering  judgment  question,  since  the  solution  obtained 
with  an  analysis  using  a  finite  length  wake  may  never  completely  converge. 
Also,  it  may  be  cost  prohibitive  to  run  the  analysis  until  far  wake  conver¬ 
gence  is  obtained.  Emphasis  was  placed  on  convergence  of  the  wake  geometry 
near  the  rotor. 
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Selection  of  a  Representative  Rotor 


To  provide  a  baseline  rotor  design  to  investigate  rotor  wake  effects  on 
the  airloads  and  perfornance  of  helicopter  rotors  in  forward  flight,  and  to 
develop  a  large  data  base  of  distorted  wake  geometries  for  the  first  level 
forward  flight  wake  generalisation,  a  physical  description  of  a  representative 
rotor  system  was  required.  In  order  to  select  a  rotor  blade  design  which  is 
representative  of  existing  rotor  systems,  the  physical  characteristics  of  nine 
helicopter  rotors  were  obtained  and  tabulated  in  Table  I.  Pour  two-bladed 
rotors  and  five  four-bladed  rotors  tiere  selected  to  establish  a  representative 
rotor.  This  information  was  primarily  provided  from  data  furnished  by  the 
government  for  this  contract  activity.  The  selection  of  these  nine  rotors  was 
based  on  the  desire  to  include  existing  military  aircraft  and  newer  generation 
aircraft  currently  entering  the  isarket  place.  Only  two-  and  four-bladed 
systems  were  considered  for  the  definition  of  a  representative  rotor  because 
the  wake  generalisation  process  was  planned  as  a  limited  first  level  effort 
which  could  be  expanded  at  a  future  date.  Included  in  the  table  are  the  rotor 
tip  speed  and  aircraft  gross  weight  used  to  establish  representative  values  of 
thrust  coefficient.  With  this  information,  the  characteristic  parameters  of  a 
representative  rotor  were  selected,  and  are  included  in  Table  I.  The  repre¬ 
sentative  rotor  is  a  two-  and  four-bladed  rotor  with  rotor  radius  and  constant 
chord  chosen  to  provide  rotor  solidity  and  aspect  ratio  values  which  are 
inclusive  of  both  the  two-  and  four-bladed  rotor  configurations  presented  in 
Table  1. 

The  values  presented  in  Table  I  are  the  baseline  values  for  the  repre¬ 
sentative  rotor.  Variations  of  these  parameters  in  the  course  of  the  study 
were  made  about  these  basline  values.  In  particular,  thrust  level  variations 
of  about  IS  to  20  percent  were  made,  resulting  in  thrust  coefficient  ranges  of 
0.0026  to  0.0039  for  the  two-bladed  configuration  and  0.00S2  to  0.0077  for  the 
four-bladed  configuration.  For  the  solidities  chosen,  the  upper  thrust 
coefficient  bounds  are  close  but  below  stall  for  the  MACA  0012  airfoil.  The 
selection  of  flight  speed  conditions  are  not  shown  in  Table  I;  however,  the 
selection  of  speeds  was  made  to  define  a  reasonable  operating  range  for  the 
rotorcreft.  The  low  to  auderate  flight  speed  range  was  selected  to  investi¬ 
gate  conditions  where  wake  effects  are  most  significant  under  normal  study 
flight  conditions.  The  advance  ratio  range  selected  was  0.05  to  0.3,  and  the 
range  of  rotor  tip  path  plane  angles  was  0.0  to  -6.0  degrees  consistent  with 
typical  flight  attitudes.  Limited  variations  in  blade  linear  twist  and  aspect 
ratio  were  also  made  over  a  range  of  -6  to  -14  degrees  and  11  to  16,  respec¬ 
tively. 
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The  rotors  presented  in  Table  I  are  representative  of  the  current  config¬ 
urations  in  the  rotorcraft  industry.  One  of  the  basic  differences  between 
several  of  the  rotors  presented  in  Table  I  is  the  Lock  nuaber,  the  ratio  of 
air  forces  to  mass  forces.  Since  the  intent  of  this  investigation  was  to 
study  the  wake  defornations  and  their  effect  on  airloads  and  performance,  the 
selection  of  a  representative  rotor  for  the  study  required  a  Lock  number  irtiich 
was  representative  of  the  rotors  in  use  today.  Due  to  the  large  differences 
between  the  Lock  numbers  of  the  rotors  in  Table  I,  a  small  sensitivity  study 
was  perforiwd  to  establish  the  effects  of  varying  the  Lock  number  on  perfor¬ 
mance  for  the  representative  rotor.  The  representative  rotor  presented  in 
Table  I  was  run  for  two  different  Lock  numbars  (4  and  10)  at  an  advance  ratio 
of  0.2  for  a  two-bladed  configuration  with  a  uniform  inflow  model,  classical 
wake  inflow  model  and  a  distorted  wake  inflow  model.  Since  the  analysis  used 
for  the  performance  predictions  was  run  trimswd  for  aero  flapping,  the  effect 
of  Lock  number  on  performance  should  be  ssull.  This  was  found  to  be  the  case 
and  the  results  are  presented  in  Table  II.  As  anticipated,  the  only  signif¬ 
icant  differences  occur  in  the  blade  coning  angle  and  blade  cyclic  controls 
required  for  aero  flapping  trim.  From  this  study  and  the  data  of  Table  I  it 
was  determined  that  a  Lock  number  of  10  would  be  representative  and  not 
significantly  influence  the  performance  results. 


Selection  of  Trim  Procedure 

A  preliminary  activity  was  performed  to  determine  the  relationship 
between  the  performance  predictions  for  an  articulated  rotor  with  various 
flapping  responses  at  the  same  tip-path  plane  trim  attitude.  This  study  was 
necessary  because  it  was  anticipated  that  rotor  tip-path  plane  attitude  would 
be  a  primary  parameter  in  the  wake  generalisation  procedure. 

The  Sikorsky  Generalised  Rotor  Performance  Analysis  (GRP)  was  used  for 
this  initial  activity  using  a  uniform  inflow  model.  Variations  in  shaft 
angle,  flapping  with  respect  to  the  shaft  and  flapping  hinge  offset  were 
investigated  for  one  flight  condition  of  a  UM-ID  rotor.  The  tip-path  plane 
attitude  was  maintained  constant  using  different  rotor  shaft  attitudes  and 
flapping  response.  Prom  this  study  it  was  confirmed  that  the  use  of  (GRP)  for 
an  articulated  rotor  with  sero  flapping  gave  performance  results  similar  to  an 
articulated  rotor  with  non-aero  flapping  where  the  appropriate  adjustments  to 
blade  cyclic  controls  were  made  to  run  at  the  same  rotor  disk  trim  attitude. 
Thus,  an  articulated  rotor  trimmed  for  aero  flapping  was  used  for  the  course 
of  this  first  level  investigation.  It  is  recognised,  however,  that  nonlinear 
effects  occur  at  low  speed  (e.g.,  y  ■  0.05)  with  nonuniform  inflow  (Ref.  32). 
Table  III  is  a  sumsMry  of  the  results  of  this  brief  trim  study. 
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Selection  of  Wake  Geometry  Parameters 


The  major  parameters  to  consider  for  rotor  wake  geometry  are  listed  in 
Table  IV.  The  identification  of  primary  parameters  which  affect  the  wake 
distortions  was  both  a  preliminary  activity  and  an  ongoing  one  during  the 
course  of  the  investigation.  This  was  because  the  determination  of  the  param¬ 
eters  which  were  significant  in  terms  of  forward  flight  wake  distortions  was 
an  objective  of  this  project.  Initial  activities  consisted  of  the  selection 
of  several  parameters  to  concentrate  on  in  the  investigation.  The  parameters 
considered  wre:  thrust  coefficient  (Cj);  rotor  advance  ratio  (u);  blade 
number  (b);  rotor  tip  path  plane  angle  (Ojpp);  blade  twist  (6^;  and  blade 
aspect  ratio  (AR). 

It  should  be  noted  chat  for  this  investigation,  rotor  disk  attitude, 

Ojpp,  was  considered  as  an  independent  parameter.  To  eliminate  the  function¬ 
ality  of  disk  attitude  in  rotor  advance  ratio,  p,  the  advance  ratio  is  defined 

as: 


V  V 
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This  effectively  assumes  a  small  angle  for  disk  attitude.  Angles  less  than  10 
degrees  are  typical  of  level  flight  for  the  low  to  moderate  airspeeds,  at 
which  wake  effects  are  important. 

In  Che  Ref.  5  investigation,  a  brief  study  for  a  single  flight  condition 
was  conducted  Co  indicate  the  sensitivity  of  wake  geometry  to  several  of  the 
parameters.  Primary  wake  parameters  and  potential  wake  generalisation 
techniques  were  identified  on  a  preliminary  basis.  Specifically,  wake  sensi¬ 
tivity  in  terms  of  thrust  level,  blade  flexibility,  rotor  disk  attitude,  blade 
cyclic  pitch,  and  viscous  and  compressibility  effects  on  the  airfoil  charac¬ 
teristics  were  investigated.  From  this  study  and  other  results  reported  in 
Ref.  5,  the  following  primary  wake  parameters  were  identified:  rotor  sdvance 
ratio,  thrust  level,  tip-path  plane  angle  and  blade  number. 

To  provide  a  first  level  wake  generalisation,  the  parameters  selected  for 
emphasis  in  this  investigation  were  Cj,  u,  Oppp,  and  the  number  of  blades.  By 
selecting  a  representative  rotor  for  the  wake  generalisation,  typical  values 
of  blade  design  parameters  were  used.  The  sensitivity  of  tip  distortions  to 
two  of  Che  blade  design  parameters,  blade  linear  twist  and  aspect  ratio,  was 
investigated,  but  these  parameters  were  not  included  independently  as  general¬ 
ized  wake  parameters.  It  was  decided  to  limit  the  investigation  to  two-  and 
four-bladed  rotors  since  these  are  the  more  common  current  configurations. 
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TABLE  IV.  WAKE  GECMETRY  PARAMETERS 


FUMDiMEMTAL  PARAMETERS  (UNDISTORTED  WAKE  PARAMETERS) 

NUMBER  OF  BLADES,  b 
ADVANCE  RATIO,  fi 
THRUST  COEFFICIENT,  Cj 
ROTOR  DISK  ATTITUDE,  a^pp 


AZIMUTH  PARAMETERS 

BLADE  AZIMUTH, 
WAKE  AGE, 


TRIM  PARAMETERS 


CYCLIC  PITCH,  Ais,  B^g 
FLAPPING,  bj^g 


ROTOR  TRIM  FORCES 
AMD  MOMENTS 


BLADE  DESIGN  PARAMETERS 


BLADE  LINEAR  TWIST,  6, 

BLADE  TAPER  RATIO,  TR^ 

BLADE  ASPECT  RATIO,  AR 
AIRFOIL  SECTION 

TIP  DESIGN  (SWEEP,  ANHEDRAL,  NON-LINEAR  TWIST,  ETC.) 


The  first  computational  efforts  were  directed  mainly  at  the  two-bladed  config¬ 
uration  to  keep  computational  costs  down  until  a  preliminary  understanding  of 
the  relationships  between  the  parameters  under  study  and  the  wake  distortions 
were  obtained.  The  range  of  thrust  coefficient  (Cj.)  to  be  considered  was 
determined  by  choosing  blade  loading  conditions  (Cj/o)  which  would  not 
involve  stalled  operation  of  the  rotor  blades  and  would  represent  realistic 
variations  in  this  parameter.  Nominal  values  of  blade  loading  chosen  were 
.06,  .08  and  .1,  with  .08  chosen  as  a  baseline  condition.  The  range  in 
advance  ratio  chosen  for  this  investigation  was  between  .05  and  .3.  This  was 
due  in  part  to  the  confirmation  during  Che  course  of  the  investigation  that 
the  distortions  of  the  tip  vortices  become  smaller  and  the  number  of  close 
blade-vortex  interactions  decrease  with  increasing  advance  ratio.  However,  it 
is  recognized  that  the  effects  of  the  wake  distortions  at  high  advance  ratio 
may  be  important  for  some  conditions  even  though  the  wake  distortions  are 
small.  Limiting  the  advance  ratio  range  also  avoided  retreating  blade  stall 
conditions  which  were  considered  beyond  the  scope  of  this  initial  investiga¬ 
tion.  Selecting  conditions  at  an  advance  ratio  as  low  as  .05  allowed  trending 
of  the  wake  geometry  distortions  with  advance  ratio  by  comparison  with  known 
tip  vortex  distortions  in  hover. 

The  selection  of  rotor  attitudes  was  based  on  test  data  available  for 
two-  and  four-bladed  geometries.  The  rotor  attitude  for  a  given  aircraft  with 
a  specific  center  of  gravity  changes  with  rotor  advance  ratio  and  thrust 
level.  To  determine  representative  tip  path  plane  attitudes  associated  with 
current  helicopters,  the  following  study  was  conducted. 

From  simple  balance  of  force  relationships  for  a  single  rotor  helicopter 
in  forward  flight,  with  no  hinge  offset,  no  trimming  horizontal  stabilizer, 
and  assuming  small  angles,  it  can  be  shown  chat  the  rotor  tip  path  angle  is 


Ojpp  ■ 
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where  D  is  the  drag  on  the  fuselage,  assumed  to  act  through  the  center  of 
gravity  and  GW  is  Che  gross  weight  of  the  helicopter.  For  any  particular 
aircraft  operating  at  a  fixed  flight  speed  with  a  fixed  center  of  gravity, 
this  relationship  implies  that. 


°TPP^T  *  constant 


IJ 
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for  small  angles;  and  thus,  for  a  given  gross  weight  the  aircraft  must  operate 
at  a  fixed  rotor  tip-path  plane  attitude.  To  check  the  validity  of  this 
relationship,  data  were  gathered  from  several  sources  for  various  aircraft. 
Some  of  the  data  came  from  flight  test  data  and  some  from  analytical  predic¬ 
tions.  Figure  12  presents  the  data  of  this  study  in  the  form  of  a  plot  of 
Cij-otipp  versus  flight  speed.  From  this  figure  it  can  be  seen  that  for  each 
selected  aircraft,  the  data  could  be  represented  by  a  smooth  single  valued 
curve  which  is  a  function  of  flight  speed.  For  similar  aircraft  (similar 
gross  weight)  the  curves  are  similar.  Also  the  effect  of  center  of  gravity  on 
the  data  presented  in  this  manner  appears  to  be  small.  Normal  operating 
ranges  of  the  rotor  tip-path  plane  attitude  with  advance  ratio  were  determined 
for  the  representative  thrust  levels  from  the  test  data  available.  From  this 
information,  values  for  the  rotor  tip-path  plane  were  chosen  tdiich  were 
compatible  with  the  rotor  configuration  (two  or  four  bladed)  and  the  operating 
conditions. 


Selection  of  Representative  Rotor  Conditions 

The  approach  taken  in  this  investigation  to  acquire  theoretical  informa¬ 
tion  for  wake  generalisation  was  to  build  a  data  base  of  performance  predic¬ 
tions  and  distorted  wake  geometries  for  representative  rotor  flight  conditions 
in  a  sequential  manner,  lAereby  each  case  was  run  and  the  results  analysed 
before  the  next  case  was  run.  This  technique  allowed  the  generation  of  a  data 
base  of  information  which  was  constructed  in  a  logical  and  consistent  manner, 
fihere  the  selection  of  the  following  case  could  be  influenced  by  the  most 
current  results  of  the  previous  cases  as  they  became  available.  Table  V  is  a 
listing  of  the  representative  rotor  cases  studied  under  the  current  investiga¬ 
tion. 
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GENERALIZED  WARE  GEOMETRY 


The  amplest  geaeral  last  ion  of  forward  flight  wake  geometry  is  the 
classical  undistorted  wake  model.  This  is  described  first,  for  it  forms  the 
base  for  the  newly  developed  generalised  wake  geometry  aiodel  in  which  the  tip 
vortex  distortions  relative  to  the  undistorted  tip  vortex  coordinates  have 
been  generalised. 


Classical  Ondistorted  Wake  Model 

Most  forward  flight  variable  inflow  methods  that  are  currently  used  in 
the  helicopter  industry  assume  that  the  rotor  wake  geometry  is  undistorted  and 
defined  from  momentum  considerations.  The  forward  flight  wake  is  modeled  by 
undistorted  skewed  helical  vortex  filaments.  An  example  of  what  will  here¬ 
after  be  referred  to  as  the  classical  undistorted  wake  model  was  presented  in 
Fig.  1.  The  equations  for  the  undistorted  wake  coordinates  of  a  point  on  a 
tip  vortex  filament  are  presented  in  Table  VI. 

A  diagram  of  the  rotor-wake  coordinate  system  is  presented  in  Fig.  10. 

The  origin  of  the  coordinate  system  is  defined  as  the  center  of  the  rotor  disk 
in  the  tip-path  plane.  The  coordinate  axes  are  oriented  with  respect  to  the 
tip-path  plane.  The  azimuthal  wake  coordinate,  ^  (or  hereafter 

referred  to  as  "wage  age"  or  "vortex  age",  is  the  age  of  a  point  on  a  vortex 
filament  measured  in  terms  of  the  blade  rotation  angle.  That  is,  it  is 
equivalent  to  the  blade  azimuth  travel  relative  to  the  blade  azimuth  at  which 
the  wake  point  was  shed  by  the  blade,  can  exceed  360  degrees  (2ir).  For 
example,  is  720  degrees  for  a  wake  point  that  is  two  wake  revolutions 
from  the  blade  that  generated  it.  It  is  noted  that  a  wake  age  parameter, 
t(>^,  was  used  for  the  hovering  rotor  wake  generalization  in  Refs.  9  and  10. 
Considering  a  vortex  filament  model  of  the  wake,  such  as  that  shown  in  Fig. 

10,  the  wake  point  is  a  point  on  a  specified  vortex  filament.  For  this  study, 
in  which  the  focus  is  on  tip  vortex  geometry,  the  wake  points  will  refer  to  a 
point  on  the  tip  vortex  filament.  For  the  distorted  wake  model,  the  tip 
vortex  filament  represents  a  rolled  up  tip  vortex  in  idiich  the  vorticity,  shed 
from  the  outer  portion  of  the  blade,  is  concentrated  in  the  core  of  a  single 
vortex  filament. 

The  coordinates  of  the  undistorted  wake  representation  are  obtained  from 
the  rotor  advance  ratio,  p,  thrust  coefficient,  C^,  and  angle  of  the  tip- 
path  plane,  o^pp*  For  example,  the  top  view  of  the  tip  vortex  filaments  is 
obtained  directly  from  the  helicoidal  path  of  the  blade  tip  as  it  translates 
at  the  velocity  Vcosojpp  and  rotates  at  the  velocity  (JR.  These  two  velocity 
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components  ere  contained  in  the  advance  ratio,  p.  The  side  view  of  the 
classical  undistorted  wake  is  defined  by  the  wake  deflection  angle,  X^ppi 
which  determines  the  skewed  linear  path  of  each  wake  point,  and  thus  also  the 
wake  boundaries.  The  wake  deflection  angle,  Xxpp»  *  function  of  advance 
ratio  and  inflow  ratio.  Inflow  ratio  is  dependent  on  the  momentum  value  of 
induced  velocity  normal  to  the  rotor  disk. 

Formulation  of  the  Generalized  Wake  Model 

A  primary  task  of  this  study  was  to  investigate  forward  flight  wake 
geometry  to  gain  insight  into  the  behavior  of  wake  distortions  with  variations 
in  selected  parameters.  The  objective  was  to  develop  a  first  level  general¬ 
ized  wake  Budel  which  could  form  the  basis  for  future  refinement,  expansion, 
and  validation.  In  this  report  section,  the  wake  distortion  characteristics 
are  demonstrated  along  with  the  relationship  of  the  distortions  to  rotor 
design  and  flight  parameters.  Tip  vortex  distortions  are  modeled  mathe¬ 
matically  and  the  distorted  wake  coefficients  of  the  wake  model  are  general¬ 
ized  as  a  function  of  the  identified  wake  definition  parameters.  This  results 
in  the  development  of  a  first  level  generalized  wake  model  which  is  coded  to 
provide  a  computer  module  for  use  in  conjunction  with  other  computer  codes 
requiring  rotor  wake  coordinates.  An  important  feature  of  the  generalized 
wake  module  is  chat  its  use  requires  insignificant  computer  time. 

Initial  activity  toward  formulating  the  mathematical  relationship 
required  to  model  Che  distorted  wake  geometry  consisted  of  studying  available 
wake  geometry  predictions  and  experimental  data.  Predicted  wake  geometries 
and  flow  visualization  techniques  applied  in  model  rotor  tests  conducted  in  a 
water  tunnel  (Ref.  20)  and  a  wind  tunnel  (Ref.  6)  have  Ishown  chat  forward 
flight  wake  geometry  is  significantly  distorted  from  the  classical  (undis- 
Corted)  geometry.  In  Fig.  5,  reproduced  frost  Ref.  6,  15  kt  and  30  kc 
distorted  tip  vortex  geometries,  predicted  at  DTRC  for  a  rocket  interference 
investigation,  arc  compared  with  Che  hovering  and  undiscorted  15  kt  tip  vortex 
geosMtries.  Rey  features  of  Che  predicted  tip  vortex  distortions  were 
confirmed  by  smoke  visualization  and  laser  vclocimeter  measurements.  The 
predicted  tip  vortex  geometry  from  Ref.  5  for  a  OH-ID  rotor  at  89  kt  is 
compared  with  the  undistortcd  geometry  in  Fig.  6  further  illustrating  the 
features  of  the  wake  distortions.  Photographs  from  Ref.  20  of  the  tip 
vortices  from  a  two-bladed  model  rotor  et  two  simulated  flight  speeds,  35  and 
90  kts,  were  reduced  to  coordinate  form  and  compared  to  coordinates  obtained 
by  applying  the  OTRC  Wake  Geometry  Analysis  to  predict  the  distorted  wake 
geometry.  These  comparisons  are  reported  in  Ref.  5,  and  although  the  test 
data  is  limited,  the  key  features  of  the  wake  geometry  appear  to  be  well 
predicted. 


Figures  13  and  14  illustrate  the  conparisons  of  the  experisiental  axial 
tip  vortex  coordinates  with  the  distorted  wake  predictions  for  the  two  advance 
ratios,  0.074  and  0.19,  corresponding  to  35  and  90  kts,  respectively.  These 
comparisons  show  that  the  analysis  predicted  the  tip  vortex  distortion 
msplitudes  and  distribution  with  wake  age,  <|i^,  reasonably  well.  However,  it 
will  be  shorn  later  in  this  report  that  the  remaining  differences  can  be 
significant  for  rotor  airload  applications. 

At  this  point,  it  is  necessary  to  define  the  coordinates  for  wake  distor¬ 
tions  as  used  in  this  report.  The  axial  wake  coordinate,  z,  is  the  coordinate 
of  a  point  on  a  vortex  filament  representing  the  distance  from  the  rotor  tip- 
path  plane.  The  lateral  and  longitudinal  wake  coordinate  axes,  x  and  y,  are 
In  the  tip-path  plane  (Fig.  10).  Distortion  is  defined  as  the  change  in  a 
wake  coordinate  from  the  classical  undistorted  wake  definition  of  that 
coordinate.  In  particular,  Che  axial  distortion  (Az)  of  a  point  on  the  tip 
vortex  is 


Az  -  z  - 


where  X  and  are  Che  inflow  ratio  and  wake  age  as  defined  earlier. 

A  graph  of  the  axial  coordinates  of  the  tip  vortex  for  varying  wake 
azimuth  position  (wage  age),  such  as  presented  in  Figs.  13  and  14,  displays 
the  distance  to  the  rotor  disk  of  all  the  points  on  the  distorted,  skewed, 
helical  tip  vortex  for  several  helical  revolutions  from  its  origin  at  the 
blade.  The  tip  vortex  coordinates  are  dependent  on  the  particular  blade 
azimuth  selected,  such  as  shown  in  Figs.  13  and  14  for  the  0,  90,  180  and  270 
degree  blade  azimuths.  The  axial  coordinates  of  the  undistorted  tip  vortex 
lie  along  a  straight  line  in  these  figures  with  slope  equivalent  to  the 
momentum  inflow  ratio,  X,  as  shown  in  Table  VI.  The  axial  distortion  of  a 
point  on  a  tip  vortex  filament  in  these  figures  is  simply  the  Az  displacement 
from  Che  undistorted  tip  vortex  line. 

In  the  predicted  and  experimental  wakes,  exemplified  in  the  figures 
presented,  some  common  features  of  tip  vortex  geometry  are  consistently 
evident.  Several  of  these  characteristic  distortion  features  are  listed 
below: 

1.  Relative  to  the  undistorted  geometry,  large  tip  vortex  distortions 
occur  in  the  axial  direction.  The  distortions  near  the  rotor  are 
significantly  smaller  in  the  in-plane  (longitudinal  and  lateral) 
directions . 
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2.  The  magnitude  of  the  distortions  generally  decrease  with  increasing 
flight  speed. 

3.  The  forward  portion  of  the  rotor  wake  tends  to  be  displaced  toward 
the  rotor,  the  aft  portion  is  displaced  downward,  and  the  sides  are 
displaced  upward. 

4.  Portions  of  tip  vortex  filaments  can  be  at  or  above  the  plane  of  the 
rotor.  Upward  tip  vortex  distortion  increases  with  decreasing 
advance  ratio,  decreasing  thrust  coefficient,  and  decreasing  rotor 
disk  attitude  (nose  down  tilt). 

The  first  task  in  generalizing  the  tip  vortex  coordinates  x,  y,  and  z  was 
to  develop  analytical  expressions  for  each  of  these  coordinates  in  terms  of 
the  blade  azimuth  position,  t|i^,  and  wake  age,  This  was  essentially 

the  task  of  developing  curve  fit  equations  for  these  coordinates  as  functions 
of  and  i|i^.  From  these  equations,  the  coordinates  of  any  point  on  a  tip 
vortex  filament  could  be  calculated  from  its  wake  age,  and  the  azimuth 
position  of  the  blade  from  tdtich  the  tip  vortex  originated,  In  the 

equations,  coefficients  of  terms  containing  these  parameters  would  be 
dependent  on  the  rotor  design  and  flight  condition  parameters. 

From  the  limited  flow  visualization  data  from  the  water  tunnel,  it  was 
observed  that  the  tip  vortex  distortions  in  the  longitudinal  and  lateral 
directions,  x  and  y,  are  generally  small  near  the  rotor  compared  to  the  axial 
distortions.  This  is  shown  in  Fig.  15.  Predicted  wake  geometries  have  shown 
agreement  with  this  observation  at  advance  ratios  above  approximately  0.1.  It 
was  thus  decided  to  concentrate  the  activity  for  the  first  level  wake  general¬ 
ization  on  the  axial  tip  vortex  distortions,  and  initially  generalize  the  top 
view  of  the  wake  (longitudinal  and  lateral  coordinates)  by  the  undistorted 
wake  coordinates  as  defined  in  Table  VI.  This  generalizes  the  x  and  y 
coordinates  in  simple  equations  in  terms  of  advance  ratio  and  number  of 
blades.  An  alternative  generalization  which  considers  the  distortions  of  the 
lateral  coordinate  will  be  described  later. 

As  can  be  seen  in  Figs.  13  and  14,  the  variation  of  the  distorted  axial 
tip  vortex  coordinate  with  wake  age  is  blade  azimuth  dependent.  Tracking  each 
axial  coordinate  peak  with  blade  azimuth,  through  one  blade  revolution,  it  can 
be  seen  that  the  peak  follows  a  consistent  pattern.  The  peak  grows  in 
amplitude  at  approximately  a  constant  rate  away  from  the  blade  that  generated 
it.  This  behavior  is  better  illustrated  if  the  axial  coordinate  distributions 
for  the  various  blade  azimuth  positions  are  superimposed  on  one  another  and 
phase  adjustments  are  made  for  blade  azimuth  position.  Figure  16  is  a  plot  of 
such  a  superposition  for  the  predicted  wake  geometries  of  Fig.  13.  From  this 
figure  the  behavior  of  the  tip  vortex  geometry  with  rotor  position  is  more 
clearly  evident.  Phasing  the  wake  age  by  the  azimuth  angle  of  the  blade  that 
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generated  the  tip  vortex,  as  in  Fig.  16,  aligns  tip  vortex  points  shed  from 
the  sane  blade  asimuth.  For  example,  the  first  set  of  vertically  aligned 
negative  peaks  in  Fig.  16,  points  A,  mere  all  shed  from  a  blade  when  it  was 
near  the  downstream  position  (tp^o  ~  330*).  In  fact,  the  aligned  positive 
peaks  in  Fig.  16,  spaced  nearly  one-half  revolution  apart,  generally  represent 
those  portions  of  the  tip  vortex  that  were  shed  when  the  blade  was  near  the 
lateral  azimuth  locations  (4’bo  *  270*),  points  B  and  C  in  Fig. 

16.  The  aligned  negatively  directed  peaks,  also  spaced  by  nearly  one-half 
revolution,  represent  portions  of  the  tip  vortex  that  were  shed  when  the  blade 
was  near  the  downstream  and  upstream  locations  “  0*  and  180*  respec¬ 

tively).  However,  the  much  lower  peaks  from  the  downstream  location,  point  D, 
are  distinguished  from  the  upstream  locations  (E)  which,  although  negatively 
directed,  remain  near  the  rotor  disk.  The  above  is  consistent  with  the  wake 
distortion  observations  noted  earlier  (up  on  sides,  down  in  rear,  near  rotor 
disk  in  front).  The  phasing  of  the  wake  coordinates  as  in  Fig.  16  effectively 
results  in  platting  them  versus  a  phased  wake  age,  4^  ~  4^  which  is  desig¬ 
nated,  This  was  found  to  be  a  useful  parameter  to  generalize  the  shape  of 

the  plotted  axial  coordinate  curve.  It  is  noted  that  4",,  ~  4^  also  appears  in 
the  undistorted  wake  equations  in  Table  VI.  is  shorn  in  Fig.  10  to  be  the 
wake  azimuth  relative  to  the  downstream  directed  x  axis. 

Recognition  of  these  consistent  wake  features  was  an  important  step 
leading  to  a  mathematical  formulation  for  the  tip  vortex  distortions.  It  was 
noted  that  the  peaks  in  the  distortions  (relative  to  the  undistorted  wake 
line)  in  Fig.  16,  and  in  similar  plots  for  other  conditions,  grow  smoothly  in 
amplitude  with  wake  age.  Graphical  analysis  showed  that  the  growth  in 
amplitude  resembled  an  exponential  function. 

Observation  of  the  above  mentioned  wake  features  was  instrumental  in  the 
selection  of  the  fundamental  mathematical  relationships  for  the  representation 
of  the  axial  distortion  coordinate.  It  was  decided  to  introduce  two  func¬ 
tions,  an  aaplitude  scaling  function,  E,  and  a  geometric  shape  function,  G, 
which  %fhen  multiplied  together  equal  the  tip  vortex  axial  distortion,  6z.  As 
will  be  discussed,  E  is  acquired  from  an  exponential  curve  fit  of  the  varia¬ 
tions  of  the  distortion  peaks  with  wake  age,  and  thus  represents  an  “envelope 
function”.  The  fundasiental  mathematical  relation  for  the  nondimensional  axial 
distortion  is  simply. 


Ax/R  -  E  •  G 


where  B  ■  envelope  function 
G  ■  shape  function 
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and  the  axial  coordinate  relation  is  retained  as, 

s/R  ■  ^  ♦  As/R 

The  intent  of  the  envelope  function  is  to  generalise  the  amplitude  of  the 
distortions  with  wake  age,  and  the  intent  of  the  shape  function  is  to  general¬ 
ise  the  characteristic  asimuthal  distribution  of  the  distortions  with  wake 
age. 

An  analytical  study  was  performed  in  which  several  mathematical  expres¬ 
sions  for  the  envelope  function  and  the  shape  function  were  postulated  and 
evaluated  to  develop  accurate  curve  fit  relations  for  the  tip  vortex  axial 
distortion.  The  study  was  based  on  distorted  wake  predictions  for  several 
actual  and  representative  rotors  and  operating  conditions.  The  study  resulted 
in  the  generalised  wake  equations  listed  on  the  following  page  in  Table  VII. 
These  equations  with  defined  coefficients  make  available  a  procedure  for 
mathematically  describing  and  calculating  the  tip  vortex  axial  distortions  at 
all  combinations  of  blade  asimuth  and  wake  age  for  selected  rotor  designs  and 
operating  conditions.  Sample  plots  of  the  axial  tip  vortex  coordinate  varia¬ 
tion  with  wake  age  and  the  corresponding  envelope  function  and  shape  function 
are  presented  in  Fig.  17  for  a  sample  rotor  and  operating  condition. 

The  independent  variable  in  the  equation  for  the  envelope  function  is  the 
previously  defined  wake  age,  The  independent  variable  in  the  equation 

for  the  shape  function  is  which  is  defined  as  the  difference  between  the 
wake  age,  of  a  point  on  the  tip  vortex  and  the  asimuth  position  of  the 
blade  from  ^ich  it  was  shed, 

As  indicated  in  Table  Vll,  the  envelope  function  is  an  exponential  func¬ 
tion  of  the  wake  age,  for  the  first  two  tip  vortex  filamnt  revolutions 
and  a  linear  function  of  thereafter.  The  shape  function  is  in  the  form 
of  a  Fourier  series  with  one  set  of  harmonic  coefficients  for  the  first  tip 
vortex  revolution  and  a  second  set  for  the  following  revolutions.  As 
indicated,  the  envelope  and  shape  function  coefficients  are  dependent  on  rotor 
design  and  operating  condition  parameters. 

Generalisation  for  the  Blade  Asimuth  and  Wake  Age  Parameters 

Tip  vortex  geometry  is  both  blade  asimuth,  and  wake  age, 
dependent.  Primary  obJecCivas  in  the  selection  of  the  envelope  and  shape 
functions  were  to  make  the  function  coefficients  independent  of  blade  asimuth 
and  waka  age,  and  generalise  the  functions  over  multiple  wake  revolutions. 
Accuracy  near  the  rotor  was  emphasised.  It  was  recognised  that  the  achieve¬ 
ment  of  these  objectives,  while  retaining  acceptable  eccuracy,  would  greatly 
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•implify  the  generalisation  and  result  in  a  large  reduction  in  the  number  of 
function  coefficients.  As  evident  in  the  generalised  equations,  these  objec¬ 
tives  were  accomplished.  Background  and  justification  for  the  generalisation 
of  the  blade  asimuth  and  wake  age  parasieters  in  the  envelope  and  shape  func¬ 
tions  are  presented  below. 

Envelope  Function 

Following  the  formulation  and  evaluation  of  several  mathematical 
relationships  for  the  envelope  function,  the  equation  in  Table  VII  was 
selected  as  one  that  provided  an  acceptable  scaling  function  for  the 
amplitudes  of  the  tip  vortex  distortions  <dtile  satisfying  the  above  objectives 
for  blade  asimuth  and  wake  age  generalisation. 

To  emphasise  accuracy  near  the  rotor,  an  exponential  function  was  formu¬ 
lated  for  the  first  two  wake  revolutions  and  a  separate  linear 

function  was  formulated  for  the  following  wake  revolutions  >  4x).  The 
independent  variable  in  both  expressions  is  wake  age,  and  both  expressions  are 
blade  asimuth  independent.  The  coefficients  and  A^  were  obtained  by  a 
least  squares  fit  through  all  of  the  positive  axial  distortion  peaks  in  the 
first  three  wake  revolutions  from  all  blade  asimuth  positions,  The 

positive  peaks  (e.g.,  see  Fig.  16)  were  selected  to  increase  the  accuracy  of 
the  tip  vortex  geometry  near  the  rotor.  It  was  found  that  the  negative  peaks 
were  still  approximated  with  this  formulation.  The  tip  vortex  regions  at  the 
negative  peaks  are  sufficiently  removed  from  the  rotor  disk  so  that  the 
influence  of  small  deviations  do  not  generally  have  a  significant  effect  on 
rotor  blade  inflow  calculations.  Exclusion  of  consideration  of  the  negative 
peaks  may  somewhat  compromise  flow  field  velocity  predictions  in  the  immediate 
vicinity  of  the  negative  peak  regions.  The  coefficient  M  for  the  second  part 
of  the  envelope  function  was  obtained  by  using  a  least  squares  fit  through  all 
of  the  positive  axial  distortion  peaks  in  the  second  and  third  wake  revolu¬ 
tions  from  all  blade  azimuth  positions.  The  coefficient  B  was  obtained  by 
matching  the  second  part  of  the  envelope  function,  idtich  is  linear  with  slope 
M,  to  the  first  part  at  two  wake  revolutions  (1*,,  *  Ax).  Least  squares 
averaging  of  the  envelope  function  coefficients  over  all  blade  azimuths 
eliminated  blade  azimuth  as  a  variable  when  the  envelope  function  is 
implemented. 

Shape  Function 

Based  upon  the  consistent  character  of  the  axial  tip  vortex  coordinate 
variation  with  wake  age  idicn  phase  shifted  by  blade  asimuth  (Fig.  16),  it  was 
assumed  that  the  envelope  function  could  be  blade  asimuth  independent,  and  the 
shape  function  would  fully  accomsmdate  the  blade  azimuth  dependence. 
Preliminary  mathematical  formulations  for  the  generalisation  of  the  tip  vortex 
distortions  were  explored  for  the  two-bladed  UH-ID,  3S  and  90  kt  conditions 
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using  Che  predictsd  distorted  tip  vortex  geoaetry  es  the  reference.  In  Fig. 
18,  the  exiel  tip  vortex  distortions  for  the  90  kt  condition  ere  plotted 
versus  wske  ege  for  a  blade  axiauth  position  of  0.0  degrees.  Division  of  the 
distortions  by  the  envelope  function  value  at  each  wake  age  results  in  the 
shape  distribution  shown  in  Fig.  18.  It  is  noted  that  the  indeterminate  shape 
distribution  value  at  the  blade  (at  9^  *  0*)  was  set  to  the  value  at  one 
wake  revolution  (9„  *  360*)  to  provide  a  periodic  formulation  for  the 
geometric  shape  function.  This  is  possible  because  the  tip  vortex  distortion 
at  sero  wake  age  will  be  forced  to  a  sero  value  at  the  blade  when  the  shape 
distribution  is  multiplied  by  the  envelope  function.  From  Fig.  18,  it  is 
evident  that  the  second  revolution  of  the  shape  distribution  is  similar  to  the 
first.  Figure  19  is  a  comparison  of  the  shape  distributions  for  various  blade 
aximuth  positions.  From  this  figure,  the  similerities  in  the  character  of  the 
curves  for  different  blade  aximuths  is  illustrated.  When  these  curves  are 
phased  for  blade  aximuth  position  and  superis^sed,  the  similarity  is  more 
clearly  demonstrated,  as  shown  in  Fig.  20.  This  shape  distribution  phasing, 
accomplished  by  shifting  the  wake  age  variation  by  the  blade  aximuth  value,  is 
similar  to  the  phasing  of  the  axial  tip  vortex  coordinates  shown  in  Fig.  16. 
The  insight  gained  from  Fig.  20,  along  with  similar  results  for  other  rotors 
and  conditions,  led  to  the  inclusion  of  the  phasing  of  the  wake  age  by  the 
blade  aximuth  position  in  the  shape  function  equation  in  Table  VII.  This 
resulted  in  eosibining  the  wake  age,  ♦w*  and  blade  axistuth,  9|,,  into  a 
single  independent  parameter  for  the  shape  function,  -  9^  -  9|,. 

In  order  to  emphasise  the  accuracy  of  the  wake  generalisation  in  the 
vicinity  of  the  rotor,  the  coefficients,  g,  of  the  Fourier  series  for  the 
shape  function  in  Table  Vll,  are  separately  acquired  for  the  first  and 
following  wake  revolutions.  For  each  revolution  the  harmonic  coefficients 
obtained  at  each  blade  aximuth  position  are  averaged.  The  averaging 
eliminates  the  blade  aximuth  dependency  when  the  coefficients  are  implemented. 
Separate  sets  of  shape  function  coefficients  are  noted  in  Table  VII  where  g* 
indicates  those  for  the  first  wake  revolution  and  g"  indicates  those  applic¬ 
able  to  all  revolutions  beyond  the  first.  The  shape  function  coefficients 
acquired  from  second  wake  revolution  information  are  assumed  to  be  applicable 
to  the  second  and  all  subsequent  revolutions.  Twelve  harmonics,  N,  were  used 
to  determine  the  shape  function  coefficients  currently  stored  in  the  computer 
module  for  the  generalised  wake  geometry. 

Assessment  of  General ixat ion  for  the  Blede 
Aximuth  and  Wake  Age  Parameters 

teconstructions  of  the  axial  tip  vortex  geometry  using  tbs  generalised 
wake  equations  are  compared  with  the  original  predicted  geometry  in  Figs.  21 
to  25  for  several  coiAinations  of  rotor  and  ^orating  condition  parameters. 
These  figures  deawnstrate  the  capability  of  the  generalised  wake  equations  to 
curve  fit  the  wdke  age  variationa  of  the  axial  coordinate  for  various  blade 


42 


azimuths.  The  range  of  parameters  in  these  sample  comparisons  includes  two- 
and  four-bladed  rotors,  advance  ratios  from  0.05  to  0.3,  thrust  coefficients 
from  0.0023  to  0.0064,  and  rotor  disk  attitudes,  oijpp>  from  -^.5  to  -4 
degrees.  These  sample  comparisons  are  generally  representative  of  the 
complete  range  investigated  for  these  variables.  The  degree  of  accuracy 
achieved  by  the  wake  generalization  for  wake  age  and  blade  azimuth  was 
sufficiently  encouraging  to  proceed  with  the  detailed  study  of  the 
generalization  for  the  primary  rotor  design  and  operating  conditimi  parameters 
to  be  discussed  in  the  fallowing  section. 


Generalization  for  Rotor  Operating  Condition 
Parameters  and  Number  of  Blades 

Following  the  establishment  of  the  generalized  wake  equations,  the  depen¬ 
dency  of  the  envelope  and  shape  function  coefficients  on  the  operating  condi¬ 
tion  paramaters  (w,  C.j>,  Oppp)  was  investigated,  and  the  generalization  for 
these  parameters  %iaa  performed  for  two-  and  four-bladed  rotors  with  the  blade 
design  of  the  “representative  rotor": 

Blade  twist  -  -10* 

Taper  Ratio  ■  1 

Aspect  Ratio  *  16 
Airfoil:  NACA  0012 

As  will  be  shown  in  following  sections,  the  rotor  disk  attitude,  Oppp* 
can  be  eliminated  as  a  distortion  parameter.  Furthermore,  thrust  coefficient, 
Cp,  will  be  shown  to  be  a  parameter  defining  the  envelope  function  and  not 
Che  shape  function.  Based  on  these  findings,  a  computer  module  was  developed 
which  contains  the  generalized  wake  equations  and  their  coefficients.  The 
coefficients  were  stored  in  tabular  form  for  interpolation  over  a  range  of 
advance  ratio,  v>  snd  thrust  coefficient,  Cp,  values.  The  envelope  function 
coefficients  were  predicted  to  increase  smoothly  with  increasing  thrust 
coefficient,  to  decrease  smoothly  with  increasing  advance  ratio,  and  to  be 
insensitive  to  nominal  rotor  disk  attitude  variations  in  steady  forward 
flight.  The  shape  function  coefficients  were  found  to  be  inaensitive  to 
thrust  coefficient  variations  and  disk  attitude  variations.  This  reduced  the 
interpolation  requirement  in  the  computer  module  to  a  bivariont  parameter 
generalisation  for  the  envelope  function  (p,  Cp)  and  a  single  parameter 
generalisation  for  the  shape  function  (p).  Separate  sets  of  coefficients  were 
included  for  two  and  four  blades,  snd  rotors  with  other  nuobers  of  blades  were 
assigned  the  coefficients  of  the  four-bladed  rotor  os  a  first  level  approxi¬ 
mation.  The  sensitivity  of  the  .:i»efficients  of  the  generalised  wake  fmctions 
to  blade  twist  and  aspect  ratio  was  explored,  thereafter,  using  the  “represen¬ 
tative  rotor"  as  a  baseline  design.  (The  effect  of  taper  ratio  on  wdke 


geometry  is  similar  to  that  of  blade  twist.  Airfoil  section,  if  uniform  in 
the  spanwise  direction,  is  primarily  accounted  for  by  the  thrust  coefficient 
generalization. ) 

These  wake  generalization  activities  and  results  related  to  the  rotor 
design  and  operating  condition  parameters  will  now  be  discussed  in  greater 
detail. 

Rotor  Disk  Attitude 


In  the  Ref.  5  investigation,  the  sensitivity  of  rotor  wake  geometry  to 
rotor  disk  attitude  (tip-path  plane  angle),  Ojpp  was  explored.  The  predicted 
effect  of  rotor  disk  attitude  on  blade  loading  and  the  tip  vortex  axial 
coordinate  variation  with  wake  age  is  shotm  in  Figs.  26  and  27,  which  are 
reproduced  directly  from  Ref.  5.  The  relatively  large  tip  vortex  coordinate 
differences  for  the  two  disk  attitudes  in  Fig.  26  were  found  to  be  primarily 
due  to  the  inflow  component  of  the  forward  velocity  (V  sin  m^pp)  which  changes 
the  undistorted  wake  deflection  angle  defined  in  Table  VI.  The  change  in  wake 
deflection  angle  influences  the  axial  positioning  of  the  undistorted  wake 
relative  to  the  rotor  disk,  as  shown  in  Fig.  26.  However,  it  was  found  that 
the  tip  vortex  distortions  relative  to  Che  undistorted  wake  location  are  not 
significantly  influenced  by  the  variation  in  rotor  disk  attitude.  This  is 
shown  in  Fig.  27.  It  was  concluded  in  Ref.  5  that  disk  attitude  is  an 
important  parameter  for  rotor  wake  geometry.  However,  its  influence  is  mainly 
on  the  fundamental  undistorted  wake  geometry  and  not  the  axial  distortions  of 
the  wake  from  that  geometry. 

In  order  to  confirm  the  conclusion  from  Ref.  5,  the  envelope  functions 
and  shape  functions  were  computed  for  a  disk  attitude  variation  at  several 
advance  ratios.  The  insensitivity  of  the  shape  function  to  disk  attitude  was 
clearly  demonstrated  at  values  of  0,  -3,  and  -6  degrees  for  advance  ratios  as 
low  as  0.05.  Insensitivity  of  the  envelope  function  was  demonstrated  within 
limited  extents  of  disk  attitude  for  fixed  thrust  level.  These  limited 
extents  were  determined  from  Fig.  12,  discussed  earlier,  to  be  within  the 
range  of  typical  values  associated  with  current  helicopters.  That  is,  using 
the  range  of  realistic  disk  attitudes  from  Fig.  12  for  several  advance  ratios 
and  thrust  coefficients,  it  was  found  that  the  effect  of  disk  attitude  varia¬ 
tions  on  the  wake  distortion  is  sufficiently  limited  to  justify  remoiving  disk 
attitude  as  a  parameter  for  defining  the  envelope  function  as  well  as  the 
shape  function. 
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Based  on  these  results,  disk  attitude  was  removed  as  an  independent 
parameter  for  the  axial  tip  vortex  distortion  (Az).  Figure  12  was  used  in  the 
selection  of  disk  attitudes  for  given  thrust  levels  used  in  calculating  the 
wake  distortion  coefficients  for  the  representative  rotor  stored  in  the  gener¬ 
alized  wake  computer  module.  It  is  noted  that,  although  Ojpp  is  not  a  param¬ 
eter  for  the  tip  vortex  axial  distortion,  it  is  a  primary  parameter  for  the 
total  axial  coordinate  in  that  it  appears  in  the  ziomentum  inflow  ratio  used  to 
define  the  undistorted  portion  of  the  coordinate  (see  Tabic  ”1). 

Advance  Ratio  and  Thrust  Coefficient 


With  the  removal  of  the  rotor  attitude  as  a  parameter  for  definition  of 
the  tip  vortex  distortions,  the  task  of  generalizing  the  wake  equation  coeffi¬ 
cients  for  advance  ratio  and  thrust  coefficient  proceeded  next.  First,  a 
preliminary  computer  program  was  coded  to  calculate  and  tabulate  the  wake 
coefficients  for  the  predicted  distorted  tip  vortex  geometries  generated  for 
the  representative  rotor  over  the  investigated  range  of  wake  definition  param¬ 
eters  (0.05  <Vi'<  0.3;  0.0026  <  Cj;  <  0.0039  for  b-2;  0.0052  £  Cf  £  0.0077 
for  b"4).  Envelope  function  and  shape  function  coefficients  were  tabulated 
for  bivariant  interpolation  on  thrust  coefficient  and  advance  ratio  (Cp, 
u').  Many  coefficients  for  the  computer  program  were  obtained  by  interpola¬ 
tion  or  extrapolation  of  the  available  coefficient  data  to  regions  within  the 
bounds  of  the  matrix  of  conditions  studied  for  which  no  data  were  calculated. 
Judgment  for  the  interpolation  and  extrapolation  task  was  based  on  trend 
observations  from  the  available  coefficient  data  and  the  expected  behavior  of 
the  coefficients  from  past  experience  and  other  data  sources.  For  example, 
the  behavior  of  the  coefficients  at  low  advance  ratio  was  influenced  by 
information  available  from  hover  distorted  wake  data.  Figure  29  presents  the 
advance  ratio  variation  of  the  envelope  functions  for  the  two-  and  four-bladed 
conf igur at  ions . 

The  above  analytical  study  resulted  in  a  generalized  wake  computer  module 
which  linearly  interpolates  between  thrust  coefficient  levels  and  advance 
ratio  to  obtain  required  envelope  coefficients.  The  geometric  shape  function 
coefficients  for  Che  first  and  second  wake  revolutions  are  also  bivariantly 
interpolated  for  advance  ratio  (u)  >nd  thrust  level  (Cj).  If  the  required 
Cj  goes  beyond  the  range  for  this  parameter  in  the  computer  module,  the 
boundary  values  are  currently  used  to  define  the  envelope  function  coeffi¬ 
cients.  For  advance  ratios  beyond  0.4,  Che  module  extrapolates  to  aero 
i  distortion  at  an  advance  ratio  of  1.0.  The  computer  module  generates  separate 

^  sets  of  coefficients  for  two-  and  four-bladed  configurations.  It  will  be 

I  shown  Chat  the  coefficients  are  less  sensitive  to  blade  number  variations,  and 

1  thus,  for  three  or  more  blades  this  first  generation  computer  sudule  currently 

I  uses  Che  coefficients  for  four  blades. 
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Figures  30  through  35  contain  plots  of  the  envelope  and  shape  functions 
which  show  the  behavior  of  these  functions  with  advance  ratio  and  thrust  level 
for  two-  and  four-bladed  rotors.  It  can  be  seen  that  the  behavior  of  the 
functions  follow  consistent  and  definable  patterns  with  the  u  and  varia¬ 
tions.  It  is  observed  that  the  functions  are  most  sensitive  to  advance  ratio 
variations.  In  fact,  the  shape  function  is  shown  in  Figs.  34  and  35  to  be 
insensitive  to  thrust  coefficient  variations.  This  is  consistent  with  the 
observation  from  Fig.  28,  reproduced  from  Ref.  5,  from  which  it  is  shown  that 
the  amplitude  and  not  the  shape  of  the  wake  age  distribution  of  the  distor¬ 
tions  is  influenced  by  Cj..  It  was  hypothesized  in  Ref.  5  that  the  tip 
vortex  axial  distortions  are  approxiiaately  proportional  to  thrust  coefficient. 
Considering  the  observed  insensitivity,  thrust  coefficient  was  eliminated  as  a 
dependent  variable  for  the  shape  function,  and  the  generalized  wake  computer 
module  was  modified  to  reflect  this. 

Humber  of  Blades 

Two  advance  ratio  conditions  were  investigated  for  two  and  four  blades  at 
the  same  thrust  level  to  explore  the  sensitivity  of  tip  vortex  geometry  to 
blade  number.  Figures  36  and  37  are  comparisons  of  the  generalized  wake 
formulations  (envelope  and  shape  functions)  for  the  two  advance  ratios,  0.05 
and  0.1,  for  two-  and  four-bladed  rotors.  These  figures  indicate  a  small  but 
significant  difference  in  envelope  function  between  two  and  four  blades.  The 
shape  functions  are  generally  similar  at  values  above  0.5,  but  notable 
differences  are  evident  for  the  lower  peaks  idiich  apply  to  points  away  from 
the  rotor.  In  Fig.  38,  variations  of  the  envelope  function  coefficients  with 
thrust  coefficient  are  presented  for  two  and  four  blades.  Although  large 
differences  in  the  coefficients  for  two  and  four  blades  are  evident  for  the 
hover  condition,  the  differences  appear  to  be  significantly  decreased  in 
forward  flight.  In  fact,  the  agreement  at  the  0.1  advance  ratio  indicated  the 
possibility  of  using  single  curves  for  Ag  and  A^  by  coiri>ining  the  two-  and 
four-bladed  coefficients  at  each  forward  flight  advance  ratio  for  interpola¬ 
tion  on  thrust  coefficient.  In  the  computer  module,  the  two-bladed  results 
are  used  for  lower  thrust  coefficients  typical  of  two-bladed  rotors,  and  four- 
bladed  results  are  used  for  the  higher  thrust  coefficients  typical  of  four- 
bladed  rotors.  If  the  thrust  level  for  a  two-bladed  condition  exceeds  the 
range  used  for  the  generalisation,  the  module  will  use  the  four-bladed 
results.  Similarly,  if  the  thrust  level  for  a  four-bladed  condition  is  toe 
low  for  the  four-bladed  generalisation,  the  module  will  use  the  two-bladed 
results.  It  is  also  noted  from  Fig.  38  that  the  curves  are  sufficiently 
ssnoth  for  reliable  interpolation  on  thrust  coefficient. 


Smooth  curves  for  the  envelope  function  coefficients  for  two  nnd  four 
blades  were  shown  in  Fig.  29,  and,  as  previously  mentioned,  a  bivariant 
interpolation  is  included  in  the  generalized  wake  module  for  thrust 
coefficient  and  advance  ratio.  Separate  sets  of  shape  function  coefficients 
were  retained  in  the  generalized  wake  module  for  two  and  four  blades.  The 
current  generalized  wake  module  thus  separately  accomskodates  two-  and  four- 
bladed  rotors.  The  four-bladed  generalization  is  currently  used  for  other 
numbers  of  blades.  Future  investigations  should  expand  the  generalization 
beyond  four  blades. 

Overall,  blade  number  appears  to  have  a  smaller  influence  on  wake 
geometry  in  forward  flight  than  in  hover,  and  its  influence  appears  to  be  less 
than  that  of  advance  ratio  and  thrust  level.  For  hover,  tip  vortex  geometry 
was  found  to  be  Cp/b  dependent.  For  forward  flight  the  dependence  seems  to 
be  more  on  the  overall  rotor  thrust,  Cp,  than  on  blade  loading,  Cp/b.  The 
rapid  transition  from  blade  loading  to  overall  thrust  dependence  is  indicated 
by  the  coefficient  in  Fig.  29  between  an  advance  ratio  of  0  and  0.05. 

Some  differences  due  to  blade  number  are  expected  considering  that  as 
blade  number  is  increased  there  are  more  tip  vortices  within  any  given  wake 
region,  and  thus,  the  tip  vortices  move  under  a  different  self  inducing  flow 
field.  Also,  the  spanwise  distribution  of  peak  blade  circulations  which 
define  the  tip  vortex  strengths  are  predicted  to  be  different  for  two  versus 
four  blades  for  the  same  integrated  thrust  level.  That  is,  the  peak  tip 
vortex  circulation  strengths  do  not  directly  scale  by  blade  number,  but  are 
also  dependent  on  other  criteria. 


Sensitivity  to  Blade  Design  Parameters 


Aspect  Ratio 

The  sensitivity  of  tip  vortex  geometry  to  blade  aspect  ratio  was  investi¬ 
gated  for  the  two-bladed  rotor  at  one  advance  ratio,  0.1,  and  at  a  constant 
thrust  coefficient,  0.0039.  The  aspect  ratio  was  varied  from  10.5  to  16.  The 
corresponding  rotor  solidity  (o  •  bc/«R  «  b/«  AR)  variation  was  0.06  to  0.04. 
Figure  39  presents  the  comparison  of  the  envelope  function  and  shape  function 
for  this  independent  aspect  ratio  variation.  From  this  figure  it  is  evident 
that  the  envelope  functions  are  essentially  identical  and  the  shape  functions 
are  similar.  In  Fig.  38,  Che  A^  and  Aj  envelope  function  coefficients  for 
this  single  condition  are  also  shown  to  be  similar  for  this  aspect  ratio 
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variation.  Although  limited  in  scope,  these  results  indicate  that  aspect 
ratio  is  at  most  a  secondary  wake  geoswtry  parameter  in  forward  flight.  It  is 
noted  that  aspect  ratio  is  not  a  parameter  in  current  hovering  rotor  wake 
generalizations.  For  hover,  solidity  changes  through  number  of  blades,  not 
aspect  ratio,  were  found  to  be  significant,  resulting  in  ** 

important  parameter  (Ref.  11). 

Blade  Twist 

Parameters  which  affect  the  spanwise  distribution  of  blade  loading  for  a 
given  thrust  level  after  the  blade  bound  circulation  distribution.  This  in 
turn  changes  the  self  induced  flow  field  produced  by  the  trailing  wake  system. 
Thus,  it  would  be  expected  that  changes  in  such  parameters  would  alter  the 
wake  geometry  for  the  saise  flight  condition.  Blade  twist  and  taper  are  two 
design  parameters  representing  spanwise  variations  in  blade  pitch  and  chord. 
Since  the  net  effect  of  changing  these  parameters  is  similar,  this  initial 
investigation  was  limited  to  exploring  the  influence  of  blade  twist  as  a 
representative  parameter  which  significantly  alters  the  blade  spanwise  circu¬ 
lation  distribution.  For  hovering  flight,  twist  distribution  is  known  to  be  a 
primary  wake  geometry  definition  parameter. 

The  sensitivity  of  tip  vortex  geometry  to  blade  twist  was  investigated 
for  the  two-bladed  rotor  at  two  forward  flight  advance  ratios  of  0.1  and  0.2 
and  two  thrust  coefficients,  0.0032  and  0.0039.  Blade  linear  twist,  Op 
values  of  -6,  -10,  and  -lA  degrees  were  investigated.  Wake  geometry  informa¬ 
tion  available  for  hover  was  used  for  comparison  and  to  assist  in  the  under¬ 
standing  of  the  behavior  of  the  tip  vortex  distortions  with  blade  twist  varia¬ 
tions.  The  sensitivity  of  the  envelope  function  to  blade  twist  is  shown  in 
Figs.  40a  and  40b  for  the  0.0032  and  0.0039  thrust  coefficient,  respectively. 
The  conditions  which  showed  the  larger  envelope  function  sensitivity  (wake 
distortion  amplitude  sensitivity)  to  twist  were  selected  for  presentation  of 
the  shape  function  comparisons  in  Figs.  41  and  42.  In  Figs.  40  through  42  it 
is  shown  that  the  sensitivity  to  twist  varies  depending  on  the  flight  condi¬ 
tion,  but  overall  it  appears  that  twist  should  be  considered  as  a  primary 
parameter  in  future  investigations.  The  current  wake  generalization  module  is 
based  on  a  linear  twist  of  -10  degrees  for  the  representative  rotor. 

Lateral  and  Longitudinal  Tip  Vortex  Distortions 

The  investigation  of  wake  distortions  was  primarily  focused  on  axial 
distortions  of  the  tip  vortex  geometry  for  reasons  noted  earlier.  This 
section  of  the  report  will  describe  the  results  of  a  preliminary  investigation 
of  the  lateral  and  longitudinal  distort irnis  and  the  formulation  of  an  initial 
awdel  to  describe  these  in-plane  ("top  view”)  distortions.  Ax  and  Ay.  This 
model  ia  programmed  into  the  computer  module  as  an  option  to  the  user. 


As  noted  in  Figs.  5  and  6,  the  Ax  and  Ay  tip  vortex  distortions  near  the 
rotor  are  relatively  small  when  compared  with  the  axial  distortions.  As. 
Figures  43  and  44  present  plots  of  the  predicted  longitudinal  and  lateral  tip 
vortex  coordinates,  respectively,  versus  wake  age  at  several  advance  ratios 
for  one  thrust  level,  one  blade  azimuth  position,  and  the  two-bladed  represen¬ 
tative  rotor.  The  differences  noted  in  these  figures  for  the  blade  azimuth 
noted  *0*)  are  similar  for  other  blade  azimuth  positions.  Figures  45 
and  46  present  the  same  results  for  the  four-bladed  rotor. 

The  most  noticeable  differences  between  the  two-  and  four-bladed  distor¬ 
tions  are  that  the  four-bladed  geometries  show  larger  predicted  in-plane 
distortions:  and  that  at  the  lowest  advance  ratio,  0.05,  the  four-bladed 
longitudinal  distortions  are  of  the  same  order  as  the  lateral  distortions. 

For  both  the  two-  and  four-bladed  geometries  the  Ax  and  Ay  distortions  become 
small  with  increasing  advance  ratio,  and  generally  the  lateral  distortions  are 
significantly  larger  than  the  longitudinal  distortions.  Based  on  this  infor¬ 
mation,  it  was  decided  to  neglect  the  longitudinal  tip  vortex  distortions. 

This  assumption,  that  the  longitudinal  distortions  are  insignificant,  is 
justifiable  for  most  rotor  airload  calculations  because  the  regions  where  the 
longitudal  distortions  differ  significantly  froir  the  undistorted  wake  occur 
beneath  the  aft  longitudinal  axis  of  the  rotor  disk  «diere  the  distorted  wake 
geometries  are  generally  displaced  downward  from  the  rcto  after  one  wake 
revolution.  Generally,  within  the  first  wake  revolution.,  the  longitudinal  tip 
vortex  distortions  are  small,  and  therefore  the  longitudal  distortions  can  be 
neglected.  Thus,  the  undistorted  longitudinal  coordinate  equation  in  Table  VI 
is  used  to  represent  the  generalized  longitudinal  coordinates. 

The  lateral  coordinates  versus  wake  age  show  the  most  distortions  on  the 
advancing  and  retreating  sides  of  the  rotor  disk.  These  are  the  regions  where 
the  tip  vortices  shed  from  each  blade  interact  to  effectively  approximate  two 
trailing  vortices  similar  to  the  fixed  wii%  trailing  vortex  system.  As  a 
result,  the  lateral  distortions  may  be  important  for  locally  predicting  the 
blade  inflow  distributions  on  the  advancing  and  retreating  sides  of  the  rotor 
disk.  Therefore,  a  brief  investigation  was  performed  to  provide  a  preliminary 
lateral  distortion  model  as  a  correction  to  the  undistorted  lateral  coordinate 
relationship  defined  in  Table  VI. 

A  mathematical  model  which  modifies  the  lateral  coordinates  of  the 
classical  wake  was  formulated  and  applied  to  the  existing  data  base  of 
predicted  tip  vortex  geometries.  This  model  scales  ths  classical  relationship 
used  to  define  the  lateral  coordinate  to  obtain  a  distortion  model. 


tfhere  y  is  the  distorted  letersl  tip  vortex  coordinate, 

is  the  classical  ondistorted  lateral  coordinate  (see  Table  VI), 
and 

S  is  the  scaling  function; 

S  -  1  -  «y  (|y„|  -  yj^)/(l  -  y^)  for  |y„|  >  y^ 

S  -  1.0  |yj  1  Tl 


and  dy  is  a  function  of  the  fora: 

dy  -  Ay  <1  -  e*V 


Ay,  Xy  and  y^^  are  coefficients  uhich  are  flight  condition  dependent.  The 
values  for  Ay,  Xy  and  y^^  vers  deterained  by  an  iterative  least  squares  curve 
fit  to  the  predicted  data  for  each  flight  condition.  The  value  for  the 
coefficient,  yjj,  nns  ^osen  to  be  constant  for  both  the  two-  and  four-bladed 
geonetries  but  varies  as  a  function  of  advance  ratio.  The  other  coefficients 
vary  with  advance  ratio  and  blade  nuaber.  It  was  found  that  the  variations  of 
the  coefficients  with  thrust  level  gave  only  slight  differences  in  the 
resulting  curve  fitted  geoswtries,  and  thus  thrust  dependence  was  osiitted. 
These  lateral  distortion  coefficients  were  progranaed  into  the  generalised 
wake  coefficient  aodule  for  use  with  the  axial  distortions  as  an  option. 

Figure  47  is  an  exaaplc  of  the  resulting  lateral  coordinates  coapared  with  the 
original  coordinates  at  several  advance  ratios  for  one  blade  axiauth  position 
of  the  four-bladed  rotor.  Figure  48  presents  siailar  plots  for  a  thrust  level 
variation  at  four  different  blade  asisMith  positions  for  one  advance  ratio. 
Two-bladod  rotor  results,  not  shown,  deaonstrate  better  agreeaent  than  the 
four-bladed  results.  This  aay  be  due  to  the  fewer  nuiiber  of  tip  vortex  inter¬ 
actions  with  the  two-bladed  rotor. 

This  alternate  lateral  distortion  aodel  is  provided  as  an  option  to  the 
undiatorted  tip  vortex  aodel.  This  aodel  should  be  used  with  the  recognition 
that  it  represents  a  preliainary  attain  to  aodel  the  lateral  distortions. 

For  aost  first  level  generalised  wake  a^lieecions,  it  is  advised  thet  the 
ondistorted  wake  Mnerslisation  be  used  for  the  lopgitudinal  Md  lateral 
coordinates  end  gjiat  attention  be  focused  on  the  ppedoainant  axial  distor* 
tions. 


General iied  Wake  Suanary 


In  sunnary,  a  first  level  generalized  wake  geometry  has  been  developed. 
This  generalisation  uses  the  undistorted  wake  equations  of  Table  VI  as  a  base 
and  increments  the  undistorted  wake  coordinates  with  the  generalized  distor¬ 
tions.  Emphasis  was  placed  on  generalising  the  predominant  axial  tip  vortex 
distortions  normal  to  the  rotor  disk.  However,  an  optional  lateral  distortion 
generalization  is  also  provided.  The  axial  tip  vortex  distortions  are 
expressed  in  equation  form  in  Table  VII  in  terms  of  the  product  of  an  envelope 
function  and  shape  function.  The  envelope  function  describes  the  amplitude 
variation  with  wake  age  and  the  shape  function  describes  the  characteristic 
distribution  of  the  tip  vortex  distortion  locations  with  wake  age  phased  by 
blade  azimuth.  The  tip  vortex  distortion  equations  were  programmed  in  a 
generalised  wake  computer  module  for  two-  and  four-bladed  rotors.  The  coeffi¬ 
cients  for  the  envelope  function  are  dependent  on  rotor  advance  ratio  and 
thrust  coefficient.  The  coefficients  for  the  shape  function  are  dependent  on 
advance  ratio.  Rotor  disk  attitude  was  eliminated  as  a  tip  vortex  distortion 
parameter,  but  is  a  primary  parameter  in  the  basic  undistorted  portion  of  the 
tip  vortex  coordinates.  Although  tip  vortex  coordinates  can  be  generated  for 
any  combination  of  advance  ratio,  thrust  coefficient,  disk  attitude,  and 
number  of  bladea,  emphasis  in  the  generalisation  was  placed  on  advance  ratios 
between  0.05  and  0.3,  and  thrust  coefficients  between  0.0025  and  0.0075  for 
two-  and  four-bladed  rotors.  The  four-bladed  wake  coefficients  are  currently 
used  as  an  estisuite  for  other  numbers  of  blades.  Complete  generalization  for 
the  blade  design  parameters  such  as  aspect  ratio,  twist,  and  taper  ratio  were 
beyond  the  scope  of  this  initial  study.  The  generalised  wake  module  was 
developed  for  a  representative  untapered  blade  design  with  an  aspect  ratio  of 
16  and  a  linear  twist  rate  of  -10  degrees.  An  exploratory  sensitivity  study 
revealed  that  tip  vortex  distortions  are  sensitive  to  blade  twist  variations 
in  forward  flight,  but  aspect  ratio  appears  to  be  a  secondary  parameter. 

The  computer  time  requirement  for  application  of  the  generalized  wake 
nodule  to  generate  tip  vortex  coordiaatee  is  insignificant  relative  to  the 
large  time  requirement  of  a  frealy  diatertiag  wake  geometry  prediction 
analysis,  for  a  typical  case  the  difference  is  1000  to  one  in  computation 

speed. 

Valw  Onomstry  CiMrta 

A  sat  of  «dka  •aomstry  clMrts  msca  davalapad  akiali  provide  a  convaniant, 
readily  aceaasibla  aanrca  far  appranimating  rater  forward  fli^  wake  geometry 
and  idaMifylag  patent ial  laeatlaaw  af  claaa  blada-vartan  paasaga.  The  charts 
ara  praaantad  as  a  aalftantainad  mfaranca  in  falvma  XX  of  this  report. 

Charts  far  both  the  claaaical  nndiatartad  waha  modal  and  ^  ganaralisad 
distorted  wake  madal  are  inclwiad.  the  tharta  provide  the  nacaaaary  infoma- 
tian  to  datarmina  tip  vortex  eaordinataa  for  a  rmga  of  parameters  sufficient 


for  most  conventional  steady  level  flight  operating  conditions.  They  can  be 
used  for  a  variety  of  applications  which  require  the  geometric  relationship 
between  the  tip  vortices  and  spatial  locations  relative  to  the  rotor  blades 
and  other  helicopter  components.  In  particular,  blade/tip  vortex  interaction 
locations  and  rotor  wake  boundaries  can  be  rapidly  defined.  A  series  of 
sample  three-views  and  isometric  views  of  tip  vortex  patterns  for  two-  and 
four-bladed  rotors  are  also  included. 

Generalised  wake  envelope  function  and  shape  function  charts  are  included 
from  which  generalized  tip  vortex  coordinates  can  be  determined  for  the  range 
of  operating  condition  parameters  investigated.  Charts  are  presented  in 
Volume  II  for  advance  ratios  of  0.05,  0.10,  0.15,  0.20,  0.30,  and  0.40  and 
thrust  coefficients  between  0.0025  and  0.0075  for  both  two-and  four-bladed 
rotors.  Sample  envelope  function  and  shape  function  charts  are  presented 
herein  in  Figs.  49  and  50.  Throu^  interpolation  within  the  charts  to  a 
specified  operating  condition  and  by  multiplying  envelope  function  and  shape 
function  values  at  prescribed  wake  age  values,  the  user  can  determine  the 
incremental  tip  vortex  distortions,  which  when  added  to  the  undistorted  axial 
coordinates,  19^  (see  Table  VI),  provide  the  distorted  tip  vortex  axial 
coordinates.  Homographs  for  expediently  determining  rotor  inflow  ratio,  X, 
and  the  undistorted  axial  coordinates  are  included  in  the  Volume  II  charts. 
Charts  are  also  included  for  the  rapid  determination  of  the  undistorted 
lateral  and  longitudinal  coordinates  based  on  the  equations  in  Table  VI. 

Samples  from  Volusm  II  of  the  three-view  and  isometric  view  of  the  tip 
vortices  based  on  the  generalised  wake  model  and  the  undistorted  wake  i»del 
are  presented  in  Figs.  51  and  52.  A  series  of  these  tip  vortex  patterns  are 
included  in  Volume  II  for  representative  variations  of  advance  ratio,  thrust 
level,  disk  attitude,  and  two  and  four  blades.  A  representative  series  for 
varying  rotor  asimuth  position  is  also  included  in  Volume  II.  Plots  of  the 
loci  of  blade-vortex  passage  locations,  based  on  undistorted  in-plane  wake 
coordinates,  are  also  contained  in  Volume  II,  from  which  a  sample  is  presented 
in  Fig.  S3.  By  superimposing  the  rotor  blades  on  the  charts  at  selected 
azisMith  positions,  the  user  can  determine  from  the  charts  the  locations  where 
the  projections  of  the  tip  vortices  into  the  rotor  disk  plane  intersect  the 
blades.  In  reality  the  vortices  are  displaced  axially  from  the  rotor  disk, 
but  these  plots  provide  the  user  with  the  potential  locations  for  close  blade- 
vortex  passages.  Actual  close  passage  locations  are  determined  through  deter¬ 
mination  of  the  tip  vortex  axial  coordinates  for  these  locations.  Hike  age, 
p^,  values  are  indicated  on  the  plots  Vhieh  identify  the  asimuth  travel  of  a 
blade  since  it  shad  the  tip  vortex  at  the  blade  passage  point.  This  informa¬ 
tion  can  then  be  used  to  determine  the  axial  position  of  the  wake.  The  series 
of  these  plots  in  Volume  11  cover  too-  thromili  si»-bladed  rotors  for  advance 
ratios  between  0.05  and  0.4.  In  addition  to  this  polar  format  plot,  charts 


for  direct  determination  of  the  blade  radial  locations  for  blade-vortex  inter¬ 
sections  and  charts  for  the  determination  of  the  blade/tip  vortex  intersection 
angles  are  presented  in  Volume  II  in  a  cartesian  format. 

Charts  for  the  determination  of  fore  and  aft  wake  boundaries  are  included 
in  Volume  II  for  both  the  undistorted  and  distorted  wake  geometries.  A  sample 
of  the  generalized  distorted  wake  boundaries  is  presented  in  Fig.  54. 

The  wake  charts  supplement  the  wake  equations  and  generalized  wake 
computer  module  by  providing  a  source  for  expediently  acquiring  approximate 
wake  geometry  information  for  a  preselected  forward  flight  rotor  operating 
condition. 


INFLUENCE  OF  HAKE  (30METRY  ON  ROTOR 
AIRLOADS  AND  PERFORMANCE 


A  primary  objective  of  this  study  was  to  apply  theoretical  wake 
geometries  to  investigate  the  influence  of  tip  vortex  distortions  on  the  aero¬ 
dynamic  related  characteristics  of  a  rotor  in  forward  flight.  The  rotor 
analysis  described  earlier  (Fig.  9) •  was  applied  using  various  combinations  of 
the  computer  analysis  components  and  wake/inflow  models  to  compute  rotor  wake 
geometry,  inflow,  airloads,  performance,  and  blade  bending  moments.  The 
Generalized  Rotor  Performance  Analysis  (GRP)  was  used  to  provide  blade  air¬ 
loads  and  rotor  performance  predictions.  The  Normal  Modes  Aeroelastic  Rotor 
Analysis  (Y200)  was  used  to  provide  blade  bending  mosient  predictions.  The 
rotor  wake/inflow  models  used  are  listed  below: 

(1)  "uniform  inflow"  based  on  a  constant  siomentum  value  of  rotor  induced 
velocity  (calculated  in  GRP  and  Y200). 

(2)  "variable  inflotr"  based  on  "classical  undistorted  wake  geometry" 
(wake  and  induced  velocities  calculated  in  the  Prescribed  Wake  Rotor 
Inflow  Program,  F389SR). 

(3)  variable  inflow  based  on  "predicted  distorted  wake  geometry"  (wake 
calculated  by  the  Rotor  Wake  Geometry  Analysis,  F506;  induced 
velocities  calculated  by  F389SR). 

(4)  variable  inflow  based  on  a  "generalized  wake  geometry"  (wake  distor¬ 
tions  calculated  by  the  Generalised  Wake  Module;  induced  velocities 
calculated  by  F389SR). 

In  wake/inflow  models  (3)  and  (4),  the  wake  was  represented  in  the  rotor 
analysis  by  distorted  tip  vortex  filaments  and  undistorted  inboard  wake  fila¬ 
ments. 


In  accordance  with  the  above  objective,  the  intent  of  this  study  was  to 
desK>nstrate  the  sensitivity  of  the  aforeswntioned  rotor  characteristics  to  the 
selection  of  rotor  wake/inflow  model  and  to  perform  u  initial  assessment  of 
Che  new  first  level  generalised  Mke.  This  was  accomplished  through  applica¬ 
tion  to  Che  "representative  rotor"  described  previously  and  the  H-34  rotor. 

The  IhH  rotor  was  selected  as  a  rotor  for  idiich  test  data  for  blade  airloads 
and  banding  moments  were  available.  In  this  report  section,  results  for  Che 
representative  rotor  «re  presented  first  to  demonstrate  the  sensitivity  of 
predicted  performance  and  airloads  to  the  wake/inflow  models.  The  R-34 
reaulto  are  presented  nest  to  demonstrate  sensitivity  and  compare  airload  and 
bemdiiy  moMsat  prediecieoa  for  the  various  wsdia/ inflow  awdals  with  test  data. 
Xt  will  be  shewn  that  the  rotor  characteristics  are  sensitive  to  wSke 


distortions  for  a  wide  range  of  flight  conditions  that  position  the  tip 
vortices  near  the  rotor,  and  the  degree  of  this  sensitivity  and  the  associated 
accuracy  of  the  distorted  wake  models  are  dependent  on  the  blade-vortex  inter¬ 
action  model  used  to  predict  local  blade  airloads  for  close  vortex  encoun¬ 
ters. 


Representative  Rotor — Influence  on  Predicted  Rotor  Power 

The  sensitivity  of  integrated  rotor  perfonaance  to  the  choice  of  wake/ 
inflow  model  is  shown  in  Table  VIII  and  Figs.  55  and  56  for  the  two-and  four- 
blade  representative  rotor  conditions  listed  earlier  in  Table  V.  Comparisons 
of  predicted  rotor  torque  at  advance  ratios  between  0.05  and  0.3  for  constant 
rotor  thrust  coefficient  are  presented  in  Fig.  55.  Comparisons  of  predicted 
rotor  torque  at  representative  thrust  coefficients  for  0.1  advance  ratio  are 
presented  in  Fig.  56. 

It  is  shown  that  rotor  power  predictions  based  on  the  classical  un- 
diatorted  wake  model  are  consistently  higher  than  those  based  on  the  uniform 
inflow  model.  This  is  in  agreement  with  the  previous  findings  as  stated  in 
Ref.  8.  Classical  wake  variable  inflow  resulted  in  a  predicted  average  torque 
coefficient  increase  of  5  percent  and  3  percent  for  the  two-and  four-bladed 
rotors,  respectively.  The  extent  of  variation  from  the  average  is  t  3  percent 
for  both  rotors.  It  is  noted  that  these  differences  would  have  been  larger  if 
a  tip  loss  factor  had  not  been  used  with  the  uniform  inflow  model  in  the 
Generalized  Rotor  Performance  Analysis  (GRP).  The  tip  loss  aK>del  in  GRP  is  a 
parabolic  tip  loss  factor  for  which  the  lift  diminishes  parabolically  to  zero 
at  the  tip  from  2  chord  lengths  inboard  of  the  tip  (approximately  0.87R  for 
the  representative  rotor  blade). 

The  use  of  variable  inflow  based  on  the  predicted  distorted  wake  model 
and  the  generalised  wake  az>del  generally  resulted  in  power  requirements  within 
a  few  percent  of  the  undistortcd  wake  results.  However,  distinct  exceptions 
were  noted  and  no  consistent  trends  are  apparent  as  the  values  are  both  higher 
and  lower  Chan  the  classical  wake  values. 

The  generalised  wake  power  predictions  relative  to  those  based  on  the 
predicted  distorted  wake  are  similar  for  sosw  conditions  and  notably  different 
for  others,  with  no  well<-dafinod  trends  axeept  that  the  larger  differences 
generally  occur  for  the  four-bladed  rotor  conditions. 

It  will  be  shown  in  the  follmKng  discuaaiom  of  the  blade  airload  results 
that  the  mixed  swisitivity  results  for  the  influence  of  tip  vortex  distortions 
on  predicted  power  ate  related  to  the  predicted  severity  of  local  blade-vortex 
interactions  associated  with  local  tip  wortex  dietortiom  at  specific  flight 
conditions. 


56 


TABLE  VIII.  INTEGRATED  PERFORMANCE  RESULTS  FOR  THE  REPRESENTATIVE  ROTOR 
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Representative  Rotor — Wake  Geometry  and  Influence 
on  Predicted  Blade  Airloads 

To  investigate  the  sensitivity  of  blade  airload  prediction  to  the  various 
wake/inflow  models,  the  predicted  wake  geometries  and  airload  distributions 
for  the  representative  rotor  conditions  were  examined.  From  these,  the 
following  three  conditions  ««ere  chosen  for  presentation  herein: 
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For  the  representative  rotor  tip  speed  of  700  fps,  the  three  advance  ratios 
correspond  to  124,  41,  and  21  kts,  respectively.  These  three  conditions  were 
selected  to  show  conditions  representative  of  various  levels  of  wake 
influence.  Results  for  two-  and  four-bladed  rotors  and  advance  ratio 
conditions  from  0.05  to  0.3  are  included.  The  two  low  advance  ratio 
conditions  were  selected  as  conditions  for  which  the  wake  is  typically  closer 
to  the  rotor  and  significant  blade-vortex  interactions  occur.  The  0.3  advance 
ratio  condition  was  selected  as  a  relatively  wake  insensitive  condition. 
Emphasis  in  the  following  discussion  will  be  placed  on  the  0.1  condition. 

Before  concentrating  on  the  detailed  wake  and  airload  results  of  each 
condition,  presented  in  Figs.  57  to  78,  a  suanary  of  the  observations  obtained 
from  these  results  are  presented.  A  significant  influence  of  tip  vortex 
distortions  on  blade  airloads  is  noted  which  is  related  to  the  local  proximity 
of  the  tip  vortices  to  the  blades.  Airloads  predicted  by  the  classical  wake 
aK>del  exhibit  sane  higher  harmonic  content  in  the  aximuthal  variations 
relative  to  the  uniform  inflow  model.  Tip  vortex  distortions  increase  the 
harsnnic  amplitudes  relative  to  those  of  the  classical  wake  modal,  and  also 
introduce  additional  harmonic  content.  This  is  directly  attributable  to  the 
predicted  displacement  of  the  tip  vortices  closer  to  the  rotor  disk.  Harmonic 
airloads  increase  with  combinations  of  operating  condition  parameters  that 
result  in  more  numerous  blade-vortex  interactions  and/or  place  the  tip 
vortices  closer  to  the  rotor;  i.e.,  in  general,  decreasing  advance  ratio,  and 
shallow  rotor  disk  attitude.  Regarding  thrust  coefficient,  the  tip  vortex 
distortions  increase  with  thrust  coefficient;  however,  the  deflection  angle  of 
the  undistortad  wake  increases,  idiich  tends  to  compensate  the  positioning  of 
the  tip  vortices  near  the  rotor.  The  use  of  the  generalised  wake  model 
results  in  airload  distributions  that  are  characteristically  similar  to  tha 


predicted  distorted  wake  model  airloads,  but  the  local  amplitudes  are  often 
different.  This  results  from  small  but  significant  changes  in  the  relative 
position  of  the  tip  vortices  to  the  blades.  The  inverse  relationship  of 
induced  airload  to  vortex  distance,  from  the  Biot-Savart  law  used  in  the 
analysis,  can  result  in  high  airload  sensitivity.  This  was  recognized  prior 
to  this  investigation  and  will  be  discussed  in  more  detail.  In  the  absence  of 
a  validated  blade-vortex  interaction  model  for  close  tip  vortex  passages  it 
was  decided  to  use  a  simple  solid  body  rotation  vortex  core  model  and  assume 
that  the  potential  flow  relation  (Biot-Savart  law)  was  applicable  outside  of  a 
vortex  core,  O.OIR  (approximately  0.15  chord)  in  diameter.  This  blade-vortex 
interaction  model  does  not  account  for  the  complex  interactive  aerodynamics 
associated  with  a  close  passage  which  can  include  expansion  and  even  bursting 
of  the  vortex  core  and  local  suction  and  flow  separation  of  the  blade  surface 
(Refs.  6,  39,  40).  However,  it  does  demonstrate  the  airload  prediction 
sensitivity  to  the  interaction  model  and  the  requirement  for  an  improved  local 
blade-vortex  interaction  model. 


First  Sample  Representative  Rotor  Condition — Advance  Ratio 


^  condition  for  which  the  airloads  were  found  to  be  relatively  insensi- 
. iv-  '0  tip  vortex  distortions  is  the  following: 


b  -  2,  p  -  0.3,  -  0.0032,  Ojpp  -  -6* 

This  condition  corresponds  to  a  124  kt  flight  speed.  The  distorted  (general¬ 
ized  wake)  and  undistorted  tip  vortex  geometry  patterns  for  one  rotor  azimuth 
position  are  shown  in  Figs.  57  and  58,  respectively.  Comparing  the  tip  vortex 
patterns  in  these  figures,  it  is  observed  that  the  generalized  wake  and 
classical  undistorted  wake  geometries  are  similar  for  this  flight  condition, 
and  the  tip  vortex  distortions  in  Fig.  57  are  shown  to  be  generally  small  and 
occur  primarily  on  the  lateral  sides  of  the  wake.  Fig.  59  shows  the  loci  of 
points  where  tip  vortices  pass  directly  beneath  the  blades  at  each  blade 
azimuth  position. 

The  variation  of  the  spanwise  integrated  blade  loading  (total  blade  lift) 
with  blade  azimuth  is  shorn  in  Fig.  60  for  this  condition.  The  integrated 
blade  loading  for  the  three  wake  models  are  similar,  but  they  are  different 
from  the  uniform  inflow  results,  particularly  on  the  advancing  side  of  the 
rotor.  Although  generally  similar,  local  airload  differences  for  the  wake 
models  were  noted  in  the  vicinity  of  the  tip  vortex  passage  from  the  opposite 
blade  (forward  locus  of  vortex  passages  in  Fig.  59). 
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The  results  for  this  flight  condition  demonstrate  the  general  insensi¬ 
tivity  of  blade  airloads  to  tip  vortex  distortions  for  typical  steady  flight 
of  a  two-bladed  rotor  at  moderately  high  advance  ratios,  where  the  disk  is 
sufficiently  inclined  to  prevent  close  blade-vortex  passages.  However,  later 
in  this  report  it  will  be  shown  for  the  four-bladed  H-34  rotor  that  local 
airload  deviations  due  to  blade-vortex  interaction  occur  on  the  advancing  side 
for  a  similar  advance  ratio  (0.29)  and  disk  attitude  (-6  degrees)  condition. 

Second  Sample  Representative  Rotor  Condition— Advance  Ratio  "0.1 

Hake  geometries  and  predicted  blade  airloads  are  presented  in  Figs.  61 
to  72  for  the  following  41  kt  flight  condition  of  the  representative  rotor; 

b  ■  2,  u  »  0.1,  C.J.  •  0.0032,  Ojpp  ■  -3* 


Because  of  a  lower  advance  ratio  than  the  previous  condition,  the  number  of 
blade-vortex  passages  is  increased  as  sho«m  in  Figs.  61  through  65.  Distorted 
tip  vortices  (Figs.  61,  63)  are  much  closer  to  the  rotor  than  the  undistorted 
vortices  (Figs.  62,  64). 

Two  rotor  azimuth  positions  were  selected  for  a  detailed  presentation  of 
the  wake  geometry  and  airload  distributions.  For  the  first  selected  rotor 
position  the  blades  are  at  0  and  180  degrees  (Figs.  61,  62),  and  for  the 
second  selected  rotor  position  the  blades  are  at  90  and  270  degrees  (Figs.  63, 
64).  The  locations  where  the  tip  vortices  pass  under  the  blades  at  these 
blade  azimuth  locations,  identified  from  Fig.  65,  are  marked  on  Figs.  61,  63, 
65,  and  66  as  (1)  through  (11).  Close  blade-vortex  passages  can  be  identified 
from  the  axial  coordinate  plots  for  the  tip  vortex  from  each  of  the  four  blade 
azimuths  in  Figs.  66a  to  66d  where  the  eleven  blade-vortex  passages  are 
labeled  (1)  through  (11)  with  the  azimuth  location  of  the  blade  that  is 
involved  in  the  blade-vortex  passage.  It  is  noted  that  the  blade  location 
points  marked  in  Figs.  66a  to  66d  are  blade  tip  locations,  and  thus, 
considering  coning,  inboard  portions  of  the  blades  are  closer  to  the  tip 
vortices  than  indicated  by  the  siarked  points.  Close  blade-vortex  inter¬ 
actions,  within  approximately  0.03R  (one-half  chord  length),  can  be  identified 
from  Figs.  66a  to  66d  at  locations  (1),  (2),  (7),  (8),  (9),  (10),  (11). 

Regarding  the  first  rotor  position  (blades  at  0*,  180*),  two  cloae  vortex 
passages,  (1)  and  (2),  are  indicated  in  Figs.  66a  and  66c.  In  Fig.  66a  it  is 
shown  that  the  tip  vortex  from  the  downstream  blade  (at  0  degrees)  is 
predicted  to  pass  close  to  the  upstream  blade  (at  180  degrees).  The  wake  age 
(Fy  *  180*)  and  blade  radial  location  (0.68R)  for  this  blade-vortex  passage 
(1)  is  acquired  from  Fig.  65.  This  tip  vortex  passage  location  (1)  is  shown 
in  the  side  view  of  Fig.  61  to  be  very  close  to  the  blade.  This  is  the 


portion  of  the  tip  vortex  that  is  emitted  from  the  upaCream  blade  tip  and  is 

distorted  up  and  close  to  the  rotor  disk  as  it  travels  back  in  time.  This 

upward  distortion  is  consistent  with  the  experisiental  observation  of  Ref.  6  in 
which  it  is  shown  that  a  blade-vortex  impact  typically  occurs  near  180  degrees 
at  low  advance  ratios.  In  Fig.  61,  the  tip  vortex  from  the  upstream  blade 
(at  180  degrees)  is  also  shown  to  be  close  to  the  upstream  blade  at  a  radial 
location  of  0.38R  (from  Fig.  65),  and  this  is  indicated  as  blade-vortex 
passage  (2).  Blade-vortex  passages  (3)  through  (6)  are  progressively  older  in 
age  and  further  from  the  rotor  ^ich  results  in  an  absence  of  close  blade- 
vortex  passages  for  the  downstream  blade  (0  degrees  azimuth). 

Regarding  the  second  selected  rotor  position  shown  in  Fig.  63  (blades  at 
90*,  270*),  five  close  blade-vortex  passages  are  evident  in  Figs.  66b  and  66d. 

These  are  identified  in  the  figures  as  passages  (7)  throu^  (11),  and  are 

shown  in  Fig.  65  to  occur  on  the  advancing  blade  (90*)  at  0.96,  0.83,  and 
0.54R  and  on  the  retreating  blade  (270*)  at  0.93  and  0.71R.  Even  considering 
a  blade  coning  angle  of  5  degrees,  all  seven  of  the  indicated  tip  vortex 
passages  are  within  an  axial  displacemeitt  of  0.03R  of  the  blade.  The  three 
close  blade-vortex  passage  locations  at  the  90  degree  blade  azimuth  are  shown 
schematically  in  Fig.  67  for  the  three  wake  models,  and  the  close  proximity 
of  the  distorted  tip  vortices  to  the  blade,  compared  to  the  undistorted  wake, 
is  clearly  evident.  It  will  now  be  shown  that,  as  would  be  anticipated,  most 
of  the  close  blade-vortex  passages  cited  above  are  predicted  to  have  a 
significant  influence  on  blade  airloads. 

The  effects  of  the  four  wake/inflow  models  on  blade  airloading  are  shown 
in  Figs.  68a  and  68b.  These  figures  show  contours  of  the  blade  airload  (lift) 
distributions  as  the  blade  rotates  around  the  rotor  disk  for  each  specified 
blade  segment  and  azimuth  location  of  the  rotor  analysis.  Although  not  shown, 
the  loading  does  diminish  to  zero  at  the  edge  of  the  rotor  disk  (blade  tip). 
Each  case  was  trimmed  to  the  same  thrust,  disk  attitude,  and  zero  flapping. 

It  is  observed  in  these  figures  that  the  wake  distortions  make  the  character 
of  the  airload  distributions  more  irregular  and  increase  the  peak  amplitudes. 
This  is  particularly  noticeable  on  the  advancing  and  retreating  sides  of  the 
disk.  The  generalized  wake  model  produces  the  same  general  characteristic 
distributions  as  the  predicted  distorted  wake  model. 

The  azimuth  variation  of  the  total  blade  lift  resulting  from  a  spanwise 
integration  of  blade  loading  is  presented  in  Fig.  69.  The  blade  lift  varia¬ 
tion  based  on  the  two  distorted  tip  vortex  representations  is  shown  to  fluc¬ 
tuate  significantly  more  than  the  two  simpler  wake/inflow  representations. 

The  harmonic  content  from  the  generalized  wake  model  is  similar  to  that  from 
the  predicted  distorted  wake  model,  but  the  amplitudes  differ,  particularly 
near  90  degrees  blade  azimuth.  These  differences  are  attributable  to  the 
slightly  different  positions  of  the  tip  vortices  of  the  distorted  and  general¬ 
ised  wake  models  which  produce  significant  changes  in  local  blade  loading. 


The  influence  of  the  different  weke/inflow  aodela  on  the  epenwise  loading 
distributions  is  shown  in  Figs.  70a  to  70d  for  the  0,  90,  160,  and  270  degree 
blade  axiauths.  These  loading  distributions  can  be  related  to  the  wake 
geoaetries  presented  in  the  earlier  figures.  In  Fig.  70a,  the  saooth  blade 

loading  distribution  of  the  downstreaa  blade  (0*)  reflects  the  absence  of 

close  blade-vortex  interactions  for  all  wake  aodels.  In  Fig.  70c,  the  blade 
loading  distribution  of  the  upatreaa  blade  (180*)  reflects  the  very  close 

blade-vortex  passage  (1)  at  0.68R  and  the  aoderately  close  passage  (2)  at 

0.38R  for  the  tip  vortex  distortion  aodels.  Nearly  equivalent  blade-to-vortex 
distances  were  predicted  by  both  the  distorted  wake  and  generalised  wake 
aodels  (Figs.  66a  and  66c),  which  resulted  in  siailar  loading  distributions 
for  this  blade  aziauth. 

In  Fig.  70b,  for  the  advancing  blade  (90*),  the  tip  vortex  distortion 
aodels  show  a  significant  change  in  the  loading  distribution  when  coapared  to 
the  results  of  the  classical  wake  and  unifora  inflow  aodels.  The  predicted 
airloads  of  the  distorted  and  generalized  wake  aodels  reflect  the  presence  of 
the  previously  identified  close  blade-vortex  passages  (7),  (8),  and  (9)  near 
0.96,  0.83,  and  0.54R.  These  blade-vortex  passages  were  presented  scheaat- 
ically  for  the  90  degree  blade  in  Fig.  67.  The  sensitivity  of  the  predicted 
blade  airloads  to  what  aight  appear  to  be  saall  differences  in  these  tip 
vortex  locations  is  shown  in  Fig.  70b.  The  predicted  blade  airloads  are  very 
sensitive  to  relative  vortex  location  in  regions  of  close  blade-vortex 
passage.  This  sensitivity  is  also  evident  in  the  effective  angle  of  attack 
and  the  blade  drag  distributions  presented  in  Figs.  71  and  72  for  this  90 
degree  aziauth  position.  This  sensitivity  is  exeaplified  further  in  the  blade 
loading  distribution  at  270  degrees  aziauth  shown  in  Fig.  70d  where  two  close 
blade-vortex  passages,  (10)  and  (11),  were  predicted  at  0.71R  and  0.93R.  As 
was  shown  in  Fig.  66b,  the  tip  vortex  location  (11),  for  both  the  distorted 
wake  and  generalized  wake  aodels,  is  in  extreae  proxiaity  to  the  270*  blade 
near  0.93R.  The  generalized  wake  produces  large  airload  gradients  in  that 
blade  region.  The  airload  distribution  of  the  predicted  distorted  wake  aodel 
does  not  exhibit  an  extreae  variation.  This  difference  is  attributable  to  the 
finite  eleaent  blade  representation  in  the  rotor  analysis  in  which  the  induced 
velocities  and  related  airloading  are  calculated  at  the  blade  eleswnt  centers. 
With  the  distance  between  blade  eleaent  centers  larger  than  the  vortex  core 
size,  a  vortex  filaaent  with  finite  core  located  near  the  blade  eleaent  center 
can  produce  an  unbalanced  induced  influence,  the  aagnitude  of  which  is 
extreaely  dependent  on  whether  or  not  the  nearest  blade  eleaent  center  lies 
within  the  vortex  core.  For  this  condition  (Fig.  70d)  the  distorted  wake 
aodel  has  actually  positioned  the  tip  vortex  within  the  vortex  core  (point 
(11)  in  Fig.  66b),  thus  significantly  reducing  the  induced  velocity  at  this 
location  and  the  influence  on  the  airloads  coapared  with  the  generalised  wake 
aodel  results.  This  sensitivity  deiwnstrates  the  requireaent  for  attention  to 


the  blade  and  wake  aegBentation  aa  well  aa  the  aerodynaalc  repreaentation  for 
the  vortex  and  ita  influence  on  blade  airloada  in  future  blade-vortex  inter¬ 
action  aodela. 
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Third  SaBple  Rapreaentative  Rotor  Condition— Advance  Ratio  ■  0.05 

To  deaonatrate  the  application  of  the  wake/ inflow  aodela  to  a  low  apeed 
condition  with  nuaeroua  blade-vortex  paaaagea  the  following  condition  waa 
aelected: 

b  -  4,  p  -  0.05,  -  0.0064,  a„p  -  -0.5* 

Thia  condition  correaponda  to  a  21  kt  flight  apeed. 

Relative  to  the  previoua  41  kt  condition,  decreasing  the  advance  ratio 
and  increasing  the  number  of  blades  by  a  factor  of  two  results  in  an  extensive 
increase  in  the  number  of  blade-vortex  passages  aa  shown  in  Figs.  73,  74,  and 
75.  The  thrust  was  increased  and  disk  attitude  decreased  to  provide  typical 
values  consistent  with  a  low  advance  ratio  condition  of  a  four-bladed  rotor. 
The  generalised  wake  geometry,  presented  in  Pig.  73,  is  highly  deformed  from 
the  classical  wake  geometry  of  Fig.  74.  The  upward  deformation  of  the  forward 
and  lateral  portions  of  the  wake  and  downward  deformation  of  the  aft  wake  is 
typical  of  low  speed  wake  geometry.  These  characteristic  features  and  the 
displacement  of  the  tip  vortices  to  locations  near  and  above  the  rotor  disk 
have  been  observed  in  model  rotor  tests  (Ref.  6). 

The  axiauth  variation  of  spanwise  integrated  blade  loading  (total  blade 
I  lift)  for  this  condition  is  presented  in  Fig.  76.  Although  some  differences 

will  be  shown  in  the  local  blade  loading,  the  classical  wake  and  uniform 
inflow  results  are  quite  similar  on  an  integrated  loading  basis.  The  results 
I  for  the  tip  vortex  distortion  models  again  exhibit  higher  harmonics  due  to 

I  blade-vortex  passages.  A  comparison  of  the  generalised  wake  results  with  the 

I  distorted  wake  results  indicates  some  changes  in  both  the  peek  loading  loca- 

Itions  and  the  loading  amplitudes  due  to  ^nerally  small  differences  in  the  tip 
vortex  locations.  Details  of  the  spanwise  lift  distribution  are  presented  in 
Figs.  77  and  78  for  two  selected  asiauth  positions,  0  and  210  degrees. 

In  Fig.  77,  it  is  shown  that  the  air loading  of  the  downstream  blade  (0 
degrees)  is  again  relatively  insensitive  to  the  wake/inflow  model,  even  at 
this  low  speed.  This  is  associated  with  the  displacement  of  the  tip  vortices 
away  from  the  rotor  in  the  central  downatreaa  region.  The  blade  loading  for 
the  210  degree  asiauth  (Fig.  78),  was  selected  as  an  example  of  a  blade  posi¬ 
tion  where  several  distorted  tip  vortices  pass  above  the  outboard  half  of  the 
blade  (vortices  of  wake  age  92,  185,  278,  369,  end  460  degrees  in  Fig.  75), 
impinge  on  the  center  of  the  blade  (552  degrees  in  Fig.  75),  nd  pass  beneath 
the  inboerd  half  of  the  blade  (641,  773,  824,  916,  and  1012  degrees  in  Fig. 
75).  These  are  the  portions  of  the  tip  vortices  that  are  shed  from  the 
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forward  adga  of  tha  rotor  diak  and  eraval  raarward  in  a  parallel  pattam  aa 
ahown  in  Fig.  73.  A  coaplax  induced  flow  diatribution  ia  produced  at  the 
blade,  particularly  in  the  aid-apan  region  where  the  tip  vortex  filaaenta  paaa 
down  throu^  the  rotor  for  thia  condition.  Tha  blade  loading  diatribution  for 
thia  asiauth  poaition  ia  preaanted  in  Fig.  78.  The  claaaical  wake  and  unifom 
inflow  aodela  are  ahown  to  coapletely  ignore  the  influence  of  the  character- 
iatic  wake  diatortion.  The  influence  of  the  cloae  blade-vortex  paaaagea  in 
the  central  blade  region  ia  evident  in  the  tip  diatortion  wodela.  A  aigni- 
ficant  difference  in  local  odd-apan  loading  between  the  generalixed  and 
predicted  diatorted  wake  aodela  ia  evident.  Thia  difference  in  not  typical  of 
all  azxButh  poaitiona,  but  thia  210  degree  axiauth  poaition  waa  purpoaely 
aelected  to  further  exeaplify  the  requireaent  for  an  iaproved  blade-vortex 
interaction  aodel  in  the  airload  analyaia.  Siailar  to  the  findinga  for  the 
270  degree  axiauth  for  the  previoua  0.1  advance  ratio  condition,  a  vortex 
filaaent  waa  predicted  to  iapinge  on  the  blade.  The  high  aenaitivity  of  the 
airload  prediction  to  aaall  differencea  in  the  vortex  core  location  relative 
to  the  neareat  preaelected  analytical  blade  atation  reaulta  in  the  airload 
differential  for  the  two  wake  aodela.  Since  cloae  blade-vortex  paaaagea  and 
blade-vortex  encountera  are  prevalent  at  low  advance  ratioa,  awareneaa  of  the 
liaitation  of  the  current  airload  prediction  aodel  for  quantitative  loading 
accuracy  ia  iaportant,  and  improved  blade-vortex  interaction  modeling  ahould 
be  atrexaed. 

Prior  to  proceeding  to  the  reaulta  of  the  H-34  rotor  inveatigation,  for 
which  teat  data  waa  available,  the  primary  i^aervationa  from  thia  initial 
repreaentative  rotor  atudy  can  be  auamarixed. 

1.  Regarding  rotor  performance  prediction,  the  claaaical  wake  calcula¬ 
tion  of  rotor  power  for  conatant  tbruat  waa  conaiatently  higher  than 
the  uniform  inflow  calculation.  The  calculated  power  valuea  from  the 
predicted  diatorted  and  generalixed  wake  modela  were  acattered  about 
the  claaaical  wake  valuea  with  no  well-defined  trend. 

2.  Aa  previoualy  reeognixed,  the  influence  of  wake  diatortion  on  blade 
airloada  ia  directly  related  to  the  cloae  proximity  of  the  tip 
vorticaa  to  the  rotor  and  the  nua^er  of  cloae  blade-vortax  inter¬ 
act  iona.  The  wake  influence  ganerally  increaaea  with  decreaaing 
advance  ratio,  decraaaing  noae  down  tilt  of  tha  rotor  diak,  and 
incraaaing  number  of  bladaa.  Incraaaing  thruat  coefficient  increaaea 
tha  wAka  deflection  angle  but  alao  increaaea  the  waka  diatortion 
toward  the  rotor,  and  thua  ita  influence  on  blade  airloada  ia  more 
difficult  to  ganaraliae. 

3^  Relative  to  tha  claaaical  undiatortad  waka  modal,  the  forward  and 
lateral  aidea  of  tha  waka  are  diatorted  toward  the  rotor  in  forward 
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flight.  Thia  results  in  the  prediction  of  close  blede-vortex 
pesseges  which  introduce  higher  heraonic  content  and  local  spanwise 
gradients  in  blade  airloads. 

A.  The  prediction  of  tip  vortex  distortions  which  are  in  close  proxinity 
to  the  blades  at  localised  regions  of  the  rotor  results  in  a  high 
sensitivity  of  predicted  eirloede  to  saell  differences  in  tip  vortex 
geoaetry.  For  such  conditions,  differences  were  observed  in  the 
airload  predictions  of  the  generalised  and  predicted  distorted  wake 
aodels  even  though  the  widce  geoaetry  appeared  to  be  approxiaated  well 
the  generalised  wake  aodel. 

5.  The  occurrence  of  close  blade-vortex  passages  and  the  high  sensi¬ 
tivity  of  blade  airloads  to  saall  deviations  of  blade  station  to  tip 
vortex  distance  in  the  airloads  analysis,  deaonstrate  the  requirenent 
for  an  iaproved  blade-vortex  interaction  aodel  and  further  wake 
generalization  refineaent. 

H-34  Rotor~Vake  Geoaetry,  Perforaance, 

Airloads,  and  Bending  Hoaents 

The  previous  section  presented  results  illustreting  the  effect  of  tip 
vortex  distortions  on  rotor  perforaance  for  the  representative  rotor.  The 
rotor  selected  was  representative  but  fictitious  end  no  test  data  were  avail¬ 
able  for  coaparative  purposes.  This  section  presents  coaparisons  between 
aeasured  results  and  the  predictions  using  the  unifora  inflow,  classical, 
distorted,  and  generalised  wake  inflow  aodels.  Three  conditions  were  investi¬ 
gated  over  an  advance  ratio  range  froa  0.064  to  0.29  to  provide  coaperative 
results  to  investigate  the  applicability  of  distorted  wake  aethodology  for 
rotor  perforaance  and  airload  predictions.  The  four-bladed  H-34  rotor  was 
selected  for  the  coaparison  because  it  had,  at  the  tiae  of  the  study,  a  large 
base  of  test  data  readily  available  in  the  open  literature  (Refs.  30  and  31). 

Three  flight  conditions  were  selected  froa  available  flight  test  and  wind 


tunnel  data  for  the  B-34  rotor. 

Advance 

Thrust 

Disk 

Airspeed 

Ratio 

Coefficient 

Attitude 

23  kt 

0.064 

0.005 

-o.6r 

48  kt 

0.129 

0.005 

-1.13* 

110  kt 

0.29 

0.0036 

-5.7* 

Flight  test  data  were 

available  in  Ref. 

30  for  the  23 

and  48  kt  conditions. 

Hind  tunnel  test  data  of  a  full-scale  rotor  were  available  in  Ref.  31  for  the 
110  kt  condition.  Baphasis  in  the  selection  of  these  test  conditions  was 
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placed  on  ineludinf  coabinationa  of  airapaod,  diak  attitude  and  thruat 
coefficient  that  would  result  in  aitnificant  w^  interaction  effects.  Selec¬ 
tion  of  Che  0.064  and  0.129  advance  ratios  provides  a  severe  test  of  distorted 
wake  nethodology  because  of  the  predicted  occurrence  of  a  lar^s  nuaber  of 
blade/vortex  interactions.  These  23  and  48  kt  conditions  are  coaaonly 
referred  to  as  transition  conditions  where  helicopters  are  known  to  experience 
irretular  control  characteriaCics  (lef.  32)  and  hi^  vibratory  excitations  due 
Co  rotor/wake  interaction.  The  0.29  advance  ratio  condition  is  a  aoderately 
high  speed  condition  (110  kc)  for  which  soae  blade-vortex  interaction  was 
evident  in  the  available  blade  airloads  test  data  for  Che  disk  attitude  and 
lower  thrust  coefficient  of  the  selected  condition. 

The  H-34  rotor  design  characteristics  and  a  listing  of  the  test  paran- 
eters  for  each  of  Che  three  operating  conditions  are  presented  in  Table  IX. 

The  rotor  analysis  was  run  using  the  unifom  inflow*  classical  undis- 
torted  wake,  predicted  distorted  wake,  and  generalised  wake  aodels  for  all 
three  conditions.  The  Generalised  Rotor  Perfonance  Analysis  was  used  to 
predict  rotor  perfomance  and  blade  airloads  based  on  the  four  wake/inflow 
models.  The  Normal  Modes  Aeroelastic  Rotor  Analysis  was  used  to  predict  blade 
bending  moments  based  on  flexible  blade  airloads  and  the  four  wake/inflow 
models.  For  each  case  in  the  prediction  procedure  the  rotor  was  trianed  to 
the  test  thrust  coefficient  and  tip-path  plane  attitude  based  on  test  values 
of  shaft  and  flapping  angles. 

H-34  Rotor  Wake  Geometry 

Tip  vortex  geometry  pettema  at  one  rotor  asimuth  position  for  the  undis¬ 
torted  and  generalised  wake  aodels  are  presented  for  the  three  H-34  operating 
conditions  in  Figs.  79  to  84.  The  general  wake  distortion  features  are 
consistent  with  Chose  Shown  earlier  for  Che  representative  rotor  and  Chose 
demonstrated  in  test  (e.g..  Refs.  5,  6,  20).  That  is,  the  distortions 
increase  with  decreasing  advance  ratio  mid  result  in  displacing  the  forward 
and  lateral  sides  of  Che  wake  upward  and  the  rear  of  the  wake  downward.  The 
upward  distortions  move  the  vortices  near  and  even  above  the  rotor,  partic¬ 
ularly  at  the  lower  advance  ratios.  The  mmber  of  blade-vortex  passages  also 
increase  with  decreasing  advance  ratio  as  shown  in  Pigs.  85  to  87.  The 
cohbinetion  of  four-blades  and  low  advance  ratio  results  in  a  particularly 
large  concentration  of  vortex  passages  near  the  rotor  for  the  0.064  advance 
ratio  condition. 


Figure  88  shows  an  example  of  the  capability  of  the  generalised  wake 
model  to  represent  the  H-34  axial  tip  vortex  distortions  predicted  with  the 
Rotor  Hake  Geometry  Analysis.  The  general  agreement  is  good  between  the 
coordinates  of  Che  two  wAe  models.  This  la  encouraging,  considering  that  the 
representative  rotor  parameters  were  used  to  develop  the  generalised  wake. 


ZABLE  IX.  B-34  EOtOE  (SmCXBBlSTXCS 


mSICAL  CHABACIEXlSnCS  OF  IHE  H-34  EOKA 


BADIUS 

MDMBEE  OF  BLADES 
CBOBD 

HIMGE  OFFSET 
TWIST 

AIBFOIL  TTFE 


28  FT 

4 

1.366  FT 
.0357E 
>8  nBGBEES 
MODIFIED  BACA  0012 


TEST  OPERATIMG  OOMDmOMS  FOR  THIS  IMVESTIGATIOH 


CASE  1 

CASE  2 

CASE  3 

RVEREMCE 

T  (ENOTS) 

23.0 

48.0 

110.0 

nt  (FFS) 

607.0 

628.0 

650.0 

M 

.064 

.129 

.29 

Cj 

.0030 

.0030 

.0036 

*TFP  (DEB) 

>0.61 

>1.13 

>5.70 

^8  (DEC) 

-.6 

>.7 

>5.0 

*1S  (DEB) 

>.012 

>.43 

>.7 

L  (LBF) 

11600 

11600 

8252 

P  (8La6/FT»*3) 

.00254 

.00212 

.002221 

*0  (FFS) 

1116 

1116 

1116 

•75 

6.326 

5.239 

6.6 

^18 

>2.216 

>1.966 

2.4 

*18 

.345 

2.579 

6.0 

6 


i 


i. 
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However,  ee  indicated  previously,  even  setall  differences  near  Che  rotor  can  | 

have  a  significant  effect  on  airload  prediction.  The  degree  of  correlation  of  \ 

the  two  wake  models  shown  in  Fig.  88  is  generally  typical  for  advance  ratios 
above  0.1.  At  low  advance  ratios,  such  as  the  0.064  condition,  greater 
differences  can  occur  which  are  partially  attributable  Co  a  more  irregular 
distorted  wake  geometry  predicted  by  the  Rotor  Hake  Geometry  Analysis. 

H-34  Rotor  Performance 

Rotor  power  was  predicted  for  the  three  H-34  conditions  using  the  four 
wake/inflow  models  and  Che  Generalised  Rotor  Performance  Analysis  (GRP).  The 
horsepower  values  are  presented  in  Table  X.  Due  to  the  complications  of 
separating  out  the  main  rotor  power  from  other  power  requirements  (tail  rotor, 
gear  box,  etc.)  in  flight  test,  rotor  power  data  were  not  available  from  Ref. 

30  for  the  two  lower  speed  conditions.  These  two  conditions  were  selected  for 
Che  airloads  data  comparison  Co  be  presented  in  the  next  report  section. 

The  uniform  inflow  model  consistently  resulted  in  the  lowest  power 
prediction.  The  classical  undistorted  wake  power  values  are  approximately  6 
percent  higher  chan  chose  for  Che  uniform  inflow  model.  An  increase  of  3  to  5 
percent  was  predicted  in  Che  preceding  representative  rotor  study.  The  power 
predictions  of  the  two  tip  distortion  models  are  again  scattered  about  the 
classical  wake  results  at  the  two  lower  speeds  due  to  Che  occurrence  of  close 
blade-vortex  passages  and  the  sensitivity  to  the  blade-vortex  interaction 
model.  The  fact  that  the  variations  between  the  power  predictions  of  the 
three  wake  models  (classical,  distorted,  and  generalised)  are  essentially 
«tquivalent  at  Che  0.29  advance  ratio  indicates  Chat  wake  distortions  are 
generally  not  significant  for  power  prediction  at  high  advance  ratio  condi¬ 
tions  with  typical  high  nose-down  disk  attitudes  associated  with  steady  level 
flight.  At  this  condition  Che  power  predictions  using  any  of  Che  wake  models 
are  closer  to  Che  test  value,  but  generally  a  6  percent  differential  remains 
between  test  and  theory.  Based  on  past  experience,  this  differential  is  not 
typical,  and  may  be  primarily  attributable  to  difficulty  in  measuring  flapping 
values  to  Che  required  degree  of  accuracy.  The  Generalised  Rotor  Performance 
Analysis,  as  routinely  used  with  Che  classical  wake  model  to  predict  power, 
consistently  agrees  with  test  data  within  approximately  three  percent.  Thus, 
the  degree  of  correlation  for  this  one  test  point  should  not  be  considered 
typical.  Emphasis  for  comparison  with  test  data  will  be  placed  on  forth¬ 
coming  airloads  predictions. 

The  collective  and  cyclic  pitch  control  parameters  required  in  Che 
analysis  to  trim  each  condition  are  presented  in  Table  XI  and  Fig.  89  for  the 
three  operating  conditions.  Control  values  fat  the  four  wiri(e/inflow  swdels 
are  presented  along  with  the  test  values.  In  addition  to  rigid  blade  results 
from  the  GIF  analysis,  flexible  blade  results  from  the  T200  Aeroelastic  Rotor 
Analysis  ara  included  for  the  distorted  wake  model.  It  is  noted  that  test  I 
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TABLE  X.  H-34  ROTOR  HORSEPOWER 


INFLOW  MODEL 

warn 

.129 

.29 

UNIPOBM  INFLOW 

WITH  PARABOLIC 

TIP  LOSS 

635.2 

468.4 

511.4 

CLASSICAL  WARE 

681.1 

499.1 

543.8 

DISTORTED  HAKE 

645.8 

447.5 

538.0 

GENERALIZED 

HAKE 

690.6 

468.1 

542.8 

TEST  RESULTS 

NA 

NA 

579.0 

TABLE  XI.  TABLE  OF  FBEDICTEO  AMD  MEASUBED  CDMTROL  PARAMETERS 
USING  RIGID  AMD  FLEXIBLE  BLADE  RESP(»ISE  AMALYSES 
FOR  VARIOUS  HAKE/ItFLON  MODELS 


CASE  I  M  -  .064 


IMFLOW  MODEL 

®75 

Ais 

Bis 

TEST  RESULTS 

6.3 

-2.2 

0.3 

UNIFORM 

8.0 

-O.l 

1.3 

CLASSICAL  HAKE 

7.9 

-1.5 

1.1 

DISTORTED  HAKE 

8.0 

-2.1 

1.3 

OMERALIZED  HAKE 

8.2 

-2.5 

0.9 

DISTORTED  HAKE 

FLEXIBU;  BLADE 

7.4 

-2.1 

1.2 

CASE  II 

li  •  . 

129 

IMFLOH  MODEL 

®75 

Ais 

Bis 

TEST  RESULTS 

5.3 

-2.0 

2.6 

UNIFORM 

6.5 

-0.4 

2.4 

CLASSICAL  HAKE 

6.4 

-1.7 

2.3 

DISTORTED  HAKE 

6.3 

-2.0 

2.6 

GENERALIZED  HAKE 

6.3 

-1.7 

2.2 

DISTORTED  HAKE 

FLEXIBLE  BLADE 

5.8 

-1.8 

2.7 

CASE  lit  -  .29 


IMFLOH  MODEL 

*75 

Ais 

*18 

TEST  RESULTS 

6.6 

2.4 

6.0 

UNIFORM 

7.2 

2.1 

5.0 

CLASSICAL  HAKE 

7.2 

1.5 

5.0 

DISTORTED  HAKE 

7.1 

1.4 

5.0 

GENERALIZED  HAKE 

7.2 

1.5 

5.0 

DI8T0RTD  HAKE 

FLEXIBLE  BLAM 

6.9 

-1.8 

5.0 

pitch  control  angle  and  flapping  angle  aeasureaents  are  difficult  to  acquire 
accurately  on  a  full-acale  rotor,  and  accuracy  liaita  around  1  degree  are  not 
unuaual,  particularly  for  flight  teat. 
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Considering  the  expected  degree  of  accuracy  of  test  control  angles,  Fig. 
i  89  is  presented  to  denonstrate  the  outstanding  trends  rather  than  the  absolute 

I  values.  The  theoretical  collective  pitch  values,  6^^,  are  similar  for  all 

i  wake/inflow  models,  but  are  consistently  higher  than  the  test  values.  All 

I  predicted  control  values  for  the  three  wake  models  are  in  agreement  at  the  110 

I  kt  condition  (y  •  0.29)  and  within  1/2  degree  at  the  23  and  48  kt  conditions 

\  (y  >  0.064  and  0.129),  with  the  exception  of  the  lateral  cyclic,  at 

'  23  kt.  The  primary  observation  in  Fig.  89  is  the  inability  of  the  uniform 

inflow  model  to  predict  the  lateral  cyclic  pitch  requirement  for  trim  at  the 
lower  speeds  and  the  significant  improvement  in  lateral  cyclic  prediction  that 
results  from  the  wake  models,  particularly  the  two  tip  distortion  models. 
Lateral  cyclic  control.  A;,,  principally  affects  the  lateral  flapping, 
bia>  which  determines  sideward  tilt  of  the  rotor  disk.  Due  to  the  90 
degree  phase  lag  of  the  response  of  an  articulated  blade,  sideward  tilt  is 
very  sensitive  to  inflow  variations  from  front  to  rear  of  the  disk.  At  low 
forward  speed,  the  fore*aft  wake  distortion  produces  a  large  fore-aft  inflow 
variation  which  ssist  be  represented  in  any  analysis  to  predict  lateral 
flapping  and  related  cyclic  pitch.  This  ^enomenon  was  previously  demon¬ 
strated  experimentally  by  Harris  (tef.  32)  and  theoretically  by  Johnson  (Ref. 
33).  The  model  test  of  Ref.  32  indicated  a  peak  response  at  an  advance  ratio 
of  0.08  where  the  wake  distortions  and  influence  on  airloads  are  severe.  At 
high  speeds  the  predominant  inflow  is  due  to  the  flight  velocity  acting 
through  the  tilted  disk,  and  the  effect  of  wake  distortion  on  airloads  and 
flapping  is  reduced.  The  wake  influence  on  blade  airloads  is  the  subject  of 
the  following  discussion. 


The  theoretical  control  values  presented  in  Table  XI  and  Fig.  89  are 
associated  with  the  analytical  trim  procedure  used.  Because  the  geometric 
relationship  between  the  wake  and  rotor  disk  is  of  fundamental  importance  for 
this  type  of  study,  the  analysis  was  run  in  such  a  manner  as  to  trim  to  the 
awasured  rotor  disk  attitude.  To  do  this,  the  shaft  angle  and  first  harmonic 
of  blade  flapping  angle  are  prescribed  in  the  analysis.  The  response  malysis 
then  alters  the  blade  cyclic  controls  to  obtain  equilibrium  of  forces  on  the 
blades  for  the  prescribed  condition.  That  is,  the  rotor  control  values  were 
varied  in  the  computer  analysis  until  the  test  thrust  coefficient  and  longitu¬ 
dinal  tip  path  plane  attitude  (Ojpp  *  \  *  ‘‘ig)  ware  achieved.  Lateral 

flapping  trim,  b^^,  was  achieved  based  on  test  values.  An  alternate  trim 
procedure  would  have  been  to  trim  to  thrust,  propulsive  force,  and  hub 
moments.  This  was  not  selected  for  this  study  because  of  the  priswry  impor¬ 
tance  of  matching  disk  attitude  for  wake  geometry  simulation  and  the  objective 
to  study  the  influence  of  wake  distortions  in  a  controlled  manner. 
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H-34  Blade  Airloading 

Azinuth  variations  of  blade  airloads  for  the  three  selected  H-34 
conditions  at  23,  48  and  110  kt  are  presented  in  Pigs.  90  through  95.  These 
figures  consist  of  tiae  history  plots  of  local  blade  lift  force  at  specified 
blade  radial  stations  (Figs.  90a, b;  92a, b;  and  94a, b)  and  spanwise  integrated 
lift  (total  blade  lift;  Pigs.  91,  93,  and  95).  The  airload  predictions  for 
the  four  wake/inflow  nodels  are  compared  with  test  data.  It  is  noted  that  the 
test  airloads  were  reconstructed  from  harmonic  data  at  specific  blade  radial 
stations  from  Refs.  30  and  31,  and  the  predicted  airloads  were  interpolated  to 
these  blade  radial  stations.  The  flight  test  data  from  Ref.  30  for  the  23  and 
48  kt  conditions  extended  to  the  0.95R  blade  station.  The  wind  tunnel  test 
data  from  Ref.  31  for  the  110  kt  condition  extended  to  the  0.99R  blade 
atation.  The  available  test  data  had  a  limited  number  of  radial  stations  (7 
to  9)  and  harmonics  (10),  combined  with  a  limited  number  of  chordwise  pressure 
transducers  (3  to  11;  5  to  7  typical)  to  integrate  to  blade  section  lift. 

This  could  result  in  the  elimination  and/or  reduction  of  local  loading  peaks 
that  occur  within  the  discrete  data  points.  It  is,  of  course,  recognized  that 
a  discrete  representation  of  the  blade  and  wake  was  also  used  for  the 
predicted  airloads. 

To  qualitatively  separate  out  the  more  direct  wake  influences  from  the 
time  histories  of  blade  airloads,  it  is  necessary  to  use  the  uniform  inflow 
results  as  a  basis  for  comparison.  The  uniform  inflow  model  does  not  include 
a  discrete  wake  representation,  and  thus  provides  airload  variations  which  are 
related  to  variations  of  the  local  non-induced  flow  velocities,  such  as  the 
one-per-rev  variation  in  the  relative  combination  of  forward  and  rotational 
speeds  and  the  associated  blade  motions.  The  figures  show  the  well-known 
character  of  the  uniform  airload  variations.  That  is,  the  low  amplitude  first 
harmonic  content  at  low  flight  speed  and  the  higher  amplitude,  predominantly 
first  and  second  harmonic  content  at  higher  speed  are  evident.  A  local  peak 
due  to  reversed  flow  is  observed  at  the  inboard  station  in  Fig.  94a.  When 
compared  to  the  test  data,  it  is  clearly  evident  that  the  airloads  predicted 
by  the  uniform  inflow  model  are  completely  devoid  of  the  airload  fluctuations 
attributable  to  the  wake. 

The  test  airloads  have  peaka  in  local  azimuth  regions  that  are  indicative 
of  blade/wake  interaction.  These  are  observed  for  each  condition  in  Figs.  90, 
92,  and  94  where  they  are  shown  to  be  primarily  located  on  the  lateral  sides 
of  the  rotor.  The  severity  of  the  local  peaks  are  greatest  for  the  two  lower 
flight  speed  conditions  (23  and  48  kt).  However,  local  peaks,  in  addition  to 
those  related  to  the  non-induced  flow,  are  evident  at  the  110  kt  condition 
where  they  are  mainly  concentrated  on  the  advancing  side  of  the  rotor.  The 
most  obvious  difference  between  the  test  airloads  for  23  and  48  kts  is  the 
difference  in  the  azimuth  locations  of  the  high  loading  gradients  and 


associated  peaks.  This  is  believed  to  be  directly  related  to  the  azimuthal 
difference  in  blade-vortex  intersection  locations  between  the  two  speeds.  It 
is  shown  in  Figs.  80  and  82  that,  at  low  speeds,  the  tip  vorticies  shed  from 
the  upstream  edge  of  the  rotor  move  above  the  disk  and  then  down  through  it. 
The  tip  vortices  thus  impact  the  blades  where  they  cross  the  disk.  As 
aientioned  earlier,  this  behavior  has  been  verified  experimentally.  For 
example,  for  the  two-bladed  rotor  of  Ref.  6,  the  tip  vortex  from  the  opposite 
blade  impacted  each  blade  near  the  upstream  azimuth  (180  degrees)  as  shown  in 
Figs.  7  and  8.  For  the  23  kt  condition,  the  test  airloads  data  appear  to  show 
that  points  of  extreaie  blade-vortex  interaction  progressively  move  inward  on 
the  blade  from  the  blade  tip  near  60  degrees  to  approximately  0.75R  near  180 
degrees  (Fig.  90a).  A  similar  progression  takes  place  on  the  retreating  side 
of  the  rotor.  Similar  behavior  is  evident  at  the  48  kt  condition  (Fig.  92a) 
except  that  it  is  concentrated  more  toward  the  blade  tip  over  a  smaller 
azimuth  range  on  the  lateral  sides  of  the  rotor.  This  implies  that  the 
regions  of  close  blade-vortex  interaction  change  accordingly.  Also,  near  60 
and  300  degrees  the  tip  vortices  are  in  a  more  parallel  alignment  with  the 
blade.  For  close  passage  regions  near  the  tip  this  may  lead  to  large  airload 
gradients  caused  by  the  abrupt  change  in  induced  flow  direction  as  the  blade 
passes  the  vortex.  The  central  aft  region  of  the  rotor,  where  the  wake  is 
displaced  down  from  the  disk,  is  the  least  affected  by  local  wake  effects. 
However,  the  integrated  wake  influence  increases  the  total  blade  lift  in  this 
aft  region. 

Regarding  the  wake  models,  the  figures  reveal  that  the  time  histories  of 
blade  airloads  progress  from  a  smooth  low  asplitude  harmonic  variation  for  the 
classical  undistorted  wake  to  a  high  amplitude  variation  with  additional 
harmonic  content  for  the  deformed  wake  sudels.  This  behavior  is  consistent 
with  the  increased  amount  of  blade/vortex  interaction  associated  with  the  more 
complex  wake  models. 

Relative  to  uniform  inflow  results,  the  classical  wake  model  improves  the 
harmonic  airload  content  correlation  with  test  data.  Compared  to  the  test 
data,  the  classical  wake  model  considerably  underpredicts  the  amplitudes  of 
the  local  airload  fluctuations  attributable  to  the  wake.  This  is  primarily 
due  to  Che  more  distant  axial  coordinate  positioning  of  the  tip  vortices  from 
the  rotor  in  regions  where  the  distorted  wake  theory  and  limited  available 
experimental  wake  geometry  data  from  model  tests  indicate  tip  vortex  distor¬ 
tions  near  the  blades. 

For  the  conditions  investigated,  the  local  measured  airload  fluctuations 
tend  to  average  out  such  Chat,  «dien  integrated  spanwise  at  each  blade  azimuth, 
Che  correlation  of  classical  wake  total  blade  lift  with  teat  data  is  much 
improved  over  the  correlation  of  local  airloads.  Although  improved,  notable 
discrepancies  still  exist  such  as  the  40  degree  phase  shift  in  Che  negative 
peak  airloading  on  the  advancing  side. 
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The  deformed  wake  models  result  in  the  best  prediction  of  the  high 
harmonic  behavior  of  the  airloads  and  the  related  locations  of  the  local 
airload  peaks  and  gradients  with  respect  to  the  test  data.  This  supports  the 
contention  that  the  local  airload  fluctuations  are  due  to  blade-vortex  inter¬ 
action.  However,  the  msplitudes  of  the  local  airload  peaks  in  regions  of 
close  blade-vortex  passage  are  frequently  not  predicted  accurately  with  the 
deformed  wake  models. 

The  harmonic  content  of  the  airloads  from  the  generalised  wake  is  similar 
to  that  from  the  predicted  distorted  wake,  and  thus  the  local  airload  fluctua¬ 
tions  generally  occur  at  similar  azimuth  locations.  In  some  regions,  notably 
the  advancing  side  for  the  48  kt  condition,  the  generalised  wake  tends  to  damp 
out  some  of  the  excessive  amplitudes  of  the  total  blade  lift  for  the  predicted 
distorted  wake  model.  This  may  well  be  a  fortuitous  feature  of  the  general¬ 
ised  wake  model  associated  with  the  fact  that  it  is  a  curve  fit  of  the  out¬ 
standing  features  of  the  deformed  wake.  However,  as  was  noted  for  the  repre¬ 
sentative  rotor  application,  it  is  also  possible  for  the  generalized  wake 
model  to  over-accentuate  the  local  airloads  for  any  tip  vortex  distortion 
perturbation  which  places  the  vortex  closer  to  but  beyond  a  vortex  core 
distance  from  a  blade  segment  center.  Of  course  this  is  also  possible  for  any 
distorted  wake  model. 

The  over-prediction  or  under-prediction  of  the  amplitudes  of  the  local 
airload  peaks  in  regions  of  close  blade-vortex  passages  are  due  to  limitations 
of  the  close  blade-vortex  interaction  theory,  coarse  discrete  element  blade/ 
wake  modeling,  and  the  high  sensitivity  of  the  airloads  to  small  blade-vortex 
displacements  for  close  passages.  This  is  exemplified  by  examining  in  greater 
detail  the  most  obvious  airload  over-prediction  at  90  degrees  azimuth  for  the 
48  kt  condition  (Fig.  92b).  The  geometric  locations  of  the  tip  vortices 
relative  to  the  90  degree  blade  are  shown  in  Fig.  96  for  each  wake  model  and 
the  corresponding  spanwise  airload  distributions  are  presented  in  Fig.  97.  A 
rear  view  is  presented  showing  axial  and  lateral  positions  of  the  local  tip 
vortex  filament  segments  that  pass  directly  under  (or  over)  the  blade.  The 
specific  blade-vortex  passage  points,  as  observed  in  the  top  view  (Fig.  82), 
are  indicated  on  each  segment. 

From  the  geometry  presented  for  the  classical  undistorted  wake  at  the  48 
kt  condition,  it  is  apparent  that  the  predicted  tip  vortex  segments  are 
reswved  from  the  blade,  and  thus  the  related  loading  distribution  in  Fig.  97 
is  such  smoother  than  those  of  the  two  deformed  wake  models. 

For  this  condition,  the  predicted  distorted  wake  places  the  tip  vortex 
segments  very  close  to  the  blade  with  one  impacted  by  the  blade.  This  is  a 
severe  situation  for  which  no  reliable  blade-vortex  interaction  model  is  known 
Co  exist  Chat  adequately  represents  the  physics  of  the  interaction.  The  over- 
predicted  amplitude  of  Che  airload  peak  in  Fig.  92b  corresponds  to  the  peA  in 
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Fig.  97  near  0.85R  produced  by  the  clone  vortex  passage.  In  fact,  a  very 
slight  displacement  could  reduce  the  peak  to  the  test  level.  A  local  airload 
peak  is  evident  in  the  test  data  nearer  to  the  blade  tip  idtich  indicates  that 
the  strong  blade-vortex  interaction  does  occur  in  the  vicinity  of  the  tip. 

The  generalised  wake  model  for  this  condition  under-predicts  the  peak  tip 
loading.  This  was  found  to  be  due  to  the  extreme  proximity  of  the  vortex 
segment  and  the  preselected  0.85R  blade  station  in  the  computer  analysis  which 
results  in  the  interaction  being  modeled  as  a  blade-vortex  core  impact  with  an 
associated  severe  reduction  in  induced  velocity.  In  addition  to  limitations 
of  the  blade-vortex  interaction  model,  the  discrete  blade  station  selection, 
which  is  somewhat  arbitrary,  becomes  involved  in  the  problem.  For  example, 
the  aerodynamics  at  the  0.85R  station  are  modeled  to  be  representative  of  a 
finite  blade  segment.  If  the  radial  station  had  been  selected  to  be  at  0.86R, 
for  example,  the  blade  station  %nuld  have  been  outside  the  prescribed  vortex 
core  diameter  of  O.OIR  and  a  high  loading  peak  would  have  been  predicted. 

This  identifies  a  needed  refinement  for  such  close  blade-vortex  encounters. 
Similar  problems  can  result  from  coarse  discrete  azimuth  increments. 

The  differences  in  mid-span  loading  between  the  wake  models  and  test  data 
at  90  degrees  are  also  due  to  the  close  blade-vortex  interaction  modeling  and 
blade-vortex  displacement.  A  rapid  change  in  the  axial  position  of  a  vortex 
segment  near  mid-span  is  shown  to  occur  in  Fig.  96  for  the  predicted  distorted 
wake  between  azimuths  of  90  and  lOS  degrees.  This  rapid  position  change 
occurs  slightly  sooner  for  the  generalized  wake  model  (between  75  and  90 
degrees).  This,  coupled  with  the  slight  but  significant  blade-vortex 
displacement  differences  of  the  adjacent  inboard  vortex  segment,  results  in 
the  mid-span  loading  differences  in  Fig.  97.  It  is  noted  that  the  impact  of 
Che  inner  vortex  segment  on  the  blade  loading  was  much  weaker  Chan  the  impact 
of  Che  segments  at  Che  blade  tip  for  both  distorted  wake  models  because  an 
analytical  blade  station  was  not  as  near  Che  impact  point.  This  further 
exemplifies  the  problem  with  discrete  point  modeling.  If  the  inflow  points  on 
the  lifting  line  were  repositioned,  the  loading  response  could  be  altered. 

The  response  could  be  improved  or  degraded  depending  on  Che  selection  of  Che 
inflow  points. 

The  above  examples  have  served  to  downstrate  that  close  blade-vortex 
interactions  do  occur  in  forward  flight  and  that  there  is  a  high  degree  of 
sensitivity  of  predicted  blade  airloads  to  wake  geosMtry  and  blade-vortex 
interaction  modeling  in  regions  of  close  proximity  of  rotor  and  wake. 

The  observed  trend  with  advance  ratio  is  that  differences  between  the 
airloads  of  the  three  wake  models  decrease  with  increasing  advance  ratio. 

This  characteristic  is  consistent  with  the  increased  influence  of  the  non- 
induced  airflow,  the  decreased  number  of  blade-vortex  passages,  the  decreased 
amplitude  of  tip  vortex  distortions,  and  the  increased  rotor-wake  displacement 


due  to  increased  rotor  tilt  and  high  forward  speed.  However,  at  the  highest 
selected  advance  ratio  condition  (0.29),  the  peak-to-peak  amplitudes  of  the 
local  airloads  predicted  with  the  three  wake  models  are  shown  in  Fig.  9A  to  be 
subdued  on  the  advancing  side  of  the  rotor  relative  to  the  test  amplitudes. 
This  implies  that  for  this  110  kt  condition,  wake  distortions  are  still 
significant  and  the  predicted  wake  may  not  be  sufficiently  close  to  the 
advancing  side  of  the  rotor. 

In  summary,  the  test  airloads  of  the  three  selected  conditions  generally 
lie  between  the  airloads  predicted  by  the  classical  undistorted  wake  sudel  and 
the  deformed  wake  models.  All  of  the  local  airloading  results  clearly  show 
that  the  use  of  a  distorted  wake  inflow  model  is  needed  for  improved  correla¬ 
tion  over  the  uniform  inflow  results.  However,  the  results  also  indicate  that 
the  current  aerodynamic  interaction  models  still  need  to  be  refined  to  some 
degree.  The  use  of  a  uniform  inflow  model  will  not  yield  the  correct  local 
blade  air loading  or  integrated  spanwise  airloadings.  This  lack  of  correlation 
is  the  least  at  the  highest  advance  ratio,  but  the  correlation  is  still 
lacking  in  the  details  of  the  problem.  The  classical  wake  results  are  clearly 
an  improvement  over  the  uniform  inflow  results,  especially  on  an  integrated 
spanwise  load  basis.  However,  this  just  demonstrates  the  fact  that  the  inte¬ 
grated  results  can  mask  pertinent  local  features  of  any  prediction  analysis. 
The  local  airloading  results  clearly  illustrate  deficiencies  in  the  classical 
wake  results  as  compared  with  the  test  measurements.  The  distorted  wake  model 
offers  an  improvement  in  the  local  airload  prediction  capabilities,  particu¬ 
larly  in  locating  the  airload  fluctuations.  However,  the  amplitudes  are  often 
overpredicted.  The  generalised  wake  results  tend  to  mimic  the  distorted  wake 
results  as  they  should,  but  the  slight  variations  in  the  tip  vortex  locations 
between  the  two  isodels  sosietimes  result  in  significant  local  loading  amplitude 
differences.  These  local  loading  differences  can  translate  into  noticeable 
variations  between  the  models  for  the  integrated  spanwise  loading  results. 
Also,  the  differences  are  amplified  by  the  current  close  blade-vortex  inter¬ 
action  modeling,  and  the  discrete  modeling  of  the  blades  and  wake. 

H-34  Blade  Bending  Moments 

Another  measure  of  the  predictive  accuracy  of  the  wake/inflow  models  is 
the  comparison  of  measured  and  predicted  bending  moments.  Bending  moment  time 
histories  are  directly  dependent  on  the  harmonic  variation  of  the  spanwise 
distribution  of  the  airloads  along  the  blade  induced  by  the  flow  field. 

Figures  98  through  1(X)  present  the  flatwise  bending  moments  at  several  span- 
wise  locations  on  the  blades,  as  predicted  for  the  different  wake/inflow 
models,  compared  with  test  data  for  the  various  advance  ratio  conditions 
studied,  low  to  high  advance  ratio  respectively.  The  T200  Aeroelastic  Rotor 
Analysis  was  used  for  the  bending  moment  prediction.  Overall,  for  all  three 
conditions,  the  results  indicate  that  the  general  character  of  the  flatwise 
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bending  moments  is  significantly  better  predicted  using  the  predicted 
distorted  wake  inflow  model  as  opposed  to  the  classical  wake  results. 

However,  local  peaks  tend  to  be  out  of  phase  and/or  exaggerated  at  some  loca¬ 
tions  when  using  the  distorted  wake  model.  The  peak-to-peak  bending  moments 
are  generally  underpredicted  using  the  classical  wake  inflow  model  at  the 
lower  speeds  as  compared  with  the  test  measurements.  This  is  an  important 
limitation  for  stress  analysis.  The  best  correlation  in  peak-to-peak  bending 
moments  is  obtained  at  the  lowest  advance  ratio  (0.064)  using  the  distorted 
wake  model.  The  correlation  in  terms  of  peak-to-peak  variation  at  this 
advance  ratio  is  quite  good  as  can  be  seen  in  Fig.  98.  At  the  0.129  advance 
ratio  the  distorted  wake  model  predictions  still  provide  the  best  correlation 
with  teat  results  in  peak-to-peak  amplitude,  but  the  predicted  results  are  now 
lower  than  the  test  results  as  sho%m  in  Fig.  99.  At  the  highest  advance 
ratio,  0.29,  the  character  of  the  asimuth  variation  of  bending  moment  is 
predicted  with  the  distorted  wake  model,  but  the  peak-to-peak  moment  is  over¬ 
predicted,  as  shown  in  Fig.  100. 

Although  the  need  for  further  analytical  refinement  is  evident,  the 
degree  to  %ihich  the  distorted  wake  model  was  able  to  reproduce  the  measured 
bending  moment  distributions  is  encouraging  and  indicative  that,  with  future 
emphasis  on  blade-vortex  interaction  siodeling,  distorted  wake  methodology  will 
provide  an  improved  predictive  tool  for  blade  stress  analysis. 


Exploratory  Blade-Vortex  Interaction  Study 

The  accurate  prediction  of  helicopter  rotor  blade  air loading  has  been  a 
continuing  but  elusive  goal  in  the  helicopter  industry.  Of  the  numerous 
phenomena  influencing  rotor  airloads  that  require  further  research,  a  primary 
one  is  blade-vortex  interactron.  For  many  flight  regimes  the  tip  vortices 
from  the  main  rotor  blades  can  pass  sufficiently  close  to,  or  even  impact, 
following  blades  such  as  to  produce  a  complex  flow  environment  idtich 
influences  both  the  local  blade  airloading  and  the  vortex  characteristics  in  a 
currently  unpredictable  manner.  Several  investigators  have  addressed  this 
problem,  but  much  still  remains  to  be  learned  of  the  fundamentals  of  this 
phenomenon  before  it  can  be  accurately  modeled  in  airload  methodn. 

The  occurrence  and  importance  of  blade-vortex  interaction  and  the 
requirement  for  improved  modeling  in  rotor  analyses  are  documented  in  numerous 
references.  For  example,  survey  papers  on  rotor  aerodynamics,  such  as  Refs.  1 
and  34  through  36,  all  point  this  out.  Johnson  in  Ref.  37  showed  that  a 
lifting  surface  representation  for  the  blede  succeeded  in  reducing  the  over¬ 
predicted  loading  peaks.  However,  he  also  concluded  that  a  very  accurate 
vortex  model  and  wake  geometry  is  required  when  the  vortex  is  close  to  a 
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blade.  That  is,  the  vortex  structural  characteristics  should  be  accurately 
represented.  Such  features  as  vortex  core  sise,  vortex  bursting  ^propagation, 
and  local  blade  wake-tip  vortex  interaction  (effects  of  the  interaction  on 
vortex  geometry  and  blade  aerodynamics)  should  be  considered.  The  requirement 
for  such  considerations  are  consistent  with  the  findings  from  UTRC  wake 
studies  listed  in  the  References. 

The  requirement  for  improved  vortex  core  modeling  in  wake/airloads 
methods  was  addressed  by  Scully  in  Ref.  15  %iith  his  distorted  wake  analysis. 

He  found  that  the  use  of  a  sufficiently  large  tip  vortex  core  radiua  reduced 
the  distorted  wake  induced  airloads.  However,  the  vortex  core  radius  required 
to  achieve  this  result  was  much  larger  than  that  observed  experimentally. 

This  led  to  the  concept  of  including  tip  vortex-core  bursting  due  to  a  close 
blade-vortex  interaction  in  the  analytical  vortex  model.  A  hypothetical  model 
of  vortex  bursting  was  generally  successful  in  reducing  the  distorted  wake 
airloads  to  levels  approximating  those  of  measured  airloads. 

Considering  that  the  model  based  on  the  propagation  of  vortex-core 
bursting  was  speculative,  as  recognized  by  Scully,  Ham  in  Ref.  38  hypothesized 
that  local  blade  separation  induced  by  a  vortex  in  the  vicinity  of  the  blade 
limits  the  maximum  incremental  airload  due  to  close  blade-tip  vortex  inter¬ 
action.  He  performed  an  experiment  to  study  blade-vortex  interaction  by 
subjecting  the  mid-section  of  a  pressure-instrumented  airfoil,  mounted 
horizontally  and  perpendicular  to  the  wind-tunnel  airstream,  to  the  impinge¬ 
ment  of  a  vortex  moving  periodically  in  the  vertical  plane  from  a  point  below 
the  airfoil  to  a  point  above  the  airfoil  and  back  again.  This  investigation 
was  limited  to  nearly  right  angle  intersections  between  the  blade  and  airfoil. 
The  results  of  this  investigation,  reported  in  Ref.  39,  indicated  that  the 
maximum  vortex-induced  lift  coefficients  are  of  the  order  of  0.2  to  0.3.  The 
vortex- induced  peak  loadings  were  found  to  be  independent  of  the  frequency  of 
vortex  impingement  (i.e.,  vortex  approach  velocity),  proportional  to  vortex 
strength,  and  only  weakly  dependent  upon  the  airfoil  angle  of  attack  and  angle 
of  yaw.  He  noted  that  these  results  have  important  implications  with  respect 
to  predictions  of  the  unsteady  loading  on  helicopter  rotor  blades  due  to 
blade-vortex  interaction.  They  show  that  the  uximum  incremental  lift  coeffi¬ 
cients  induced  on  the  blade  by  a  returning  tip  vortex  m^  be  limited  to  values 
of  the  order  of  0.3  or  less,  and  suy  be  relatively  unaffected  by  vortex 
approach  velocity,  blade  mean  lift  coefficient,  or  blade  radial  flow. 

Although  Ham's  results  are  significant  in  identifying  local  vortex-induced 
blade  separation  as  a  potential  mechanism,  additional  fundamental  studies 
remain  to  be  conducted  to  demonstrate  this  for  a  helicopter  rotor. 
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As  has  been  noted  herein  in  the  example  conditions  investigated,  the 
local  blade  loading  is  predicted  to  be  significantly  influenced  by  the  induced 
inflow  of  a  near  tip  vortex.  As  indicated  by  comparison  with  H-34  rotor  test 
data,  this  induced  influence  is  over-predicted  for  many  of  these  encounters. 
The  analytical  formulation  for  this  investigation  was  based  on  quasi-steedy 
aerodynamics,  a  lifting  line/2D  airfoil  data  representation  of  the  blade,  and 
a  simple  vortex  model.  Lifting  surface  siodels,  when  applied  to  close  blade- 
vortex  encounters,  will  distribute  the  induced  influence  over  the  blade 
surface  rather  than  concentrate  it  at  coarsely  selected  collocation  points. 
However,  current  inviscid  lifting  surface  sMdels  are  not  sufficient  for  close 
blade-vortex  interaction  modeling  because  nonlinear  viscous  phenomena,  such  as 
local  blade  stall  due  to  vortex  induced  separation,  are  neglected.  The  simple 
vortex  siodel  used  for  the  airload  comparisons  presented  herein  is  based  on 
potential  flow  theory  for  the  induced  flow  outside  a  vortex  core  of  prescribed 
diameter  which  was  selected  without  regard  to  the  influence  of  a  close  blade 
passage  (no  vortex  core  expansion,  bursting,  etc.).  Recognizing  that  the 
developsient  of  a  consistently  accurate  and  physically  representative 
analytical  blade-vortex  interaction  model  is  the  subject  of  ongoing  and  future 
studies,  a  simple  blade-vortex  interaction  model  was  used  herein  to  demon¬ 
strate  the  limitations.  However,  an  exploratory  study  of  an  empirical  blade- 
vortex  interaction  siodel  was  also  performed  as  described  below. 

An  exploratory  blade-vortex  interaction  model  was  hypothesized  based  on 
the  above  mentioned  empirical  results  obtained  by  Ham  which  indicated  a 
limiting  change  in  the. lift  coefficient  (h  >  0.2  to  0.3)  due  to  the 
presence  of  an  isolated  vortex  near  a  fixed  wing  in  a  wind  tunnel  (Ref.  39), 
The  limit  in  the  change  in  lift  coefficient  is  attributed  to  the  flow 
separating  on  the  blade  in  the  presence  of  the  vortex.  Althou^  the  fixed 
wing/vortex  interaction,  as  tested  in  Ref.  39,  is  not  completely  representa¬ 
tive  of  unsteady  vortex  encounters  which  occur  for  a  rotor,  it  is  supportive 
of  the  hypothesis  of  including  a  limit  in  the  change  in  local  blade  section 
lift  coefficient  in  the  rotor  analysis.  Thus,  a  blade-vortex  interaction 
model  was  incorporated  in  the  Prescribed  Hake  Rotor  Inflow  Program  (F389SR) 
and  the  Generalized  Rotor  Performance  Analysis  (GRP)  which  allows  the  user  to 
prescribe  the  allowable  change  in  blade  section  lift  coefficient  per  change  in 
blade  azimuth  position.  The  limit  on  the  change  in  lift  coefficient  was 
formulated  in  Che  analysis  in  a  manner  such  that  Che  change  is  applied  only  to 
the  portion  of  the  local  lift  due  to  the  vortex  induced  flow  field  as  it 
changes  with  blade  azimuth.  The  analytical  model  does  not  identify  individual 
vortex  encounters,  but  assumes  Chat  a  rapid  change  in  the  local  induced  inflow 
which  produces  Che  lift  coefficient  change  is  primarily  attributable  to  an 
individual  vortex  encounter. 


In  a  brief  exploratory  inveatigation,  the  eapirical  blade-vortex  inter¬ 
action  model  waa  applied  to  the  0.129  advance  ratio  condition  for  the  H-34 
rotor  uaing  the  predicted  diatorted  wake  geometry  and  aeveral  different 
preacribed  limits  for  the  change  in  lift  coefficient  per  change  in  azimuth 
position  (0.1  to  0.3  per  15  degree  azimuth  increment).  The  resulting  air  load 
predictions  are  compared  with  test  data  and  the  previously  presented  results 
with  no  lift  coefficient  limit  in  Figs.  101  and  102.  As  anticipated,  applica¬ 
tion  of  the  new  blade-vortex  interaction  limit  model,  hereafter  "the  BVI  limit 
model",  did  result  in  decreasing  the  local  airload  peaks  on  the  advancing  side 
of  the  rotor  where  close  blade-vortex  passages  were  predicted  (Fig.  42).  The 
BVI  limit  of  0.2  generally  produced  predicted  airloads  which  compared  most 
favorably  with  the  teat  airloads,  particularly  the  local  airloads  in  Fig.  101. 
The  peak  amplitude  of  the  spanwise  integrated  airloading  on  the  advancing  side 
was  improved  but  still  somewhat  over-predicted  (Fig.  102).  A  BVI  limit  of  O.I 
reduced  the  local  airloads  more  severely  than  required  to  meet  the  test  levels 
on  the  retreating  side  of  the  rotor  disk.  Overall,  the  BVI  limit  model 
succeeded  in  modifying  the  amplitudes  of  the  peak  airloads  on  the  advancing 
side,  but  improvement  of  the  general  airload  distribution  is  not  evident.  In 
particular,  the  correlation  between  test  results  and  the  predicted  results  on 
Che  recreating  side  was  degraded.  Additional  brief  investigations,  in  which 
the  BVI  limit  model  was  applied  for  the  generalized  wake  geometry  and  a  lower 
advance  ratio  condition  (0.064),  produced  airload  changes  which  were  inconsis¬ 
tent  in  their  degree  of  correlation  with  test  data. 

Because  of  the  degradation  of  the  correlation  on  the  retreating  side  of 
Che  rotor  disk  with  Che  "BVI  limit  model",  further  investigation  of  the  H-34, 
0.129  advance  ratio  condition  was  perforamd  using  ocher  BVI  models.  First, 

Che  interaction  model  of  Johnson  (Refs.  17  and  37)  was  used.  This  model  is 
based  on  the  ratio  of  lifting  surface  to  lifting  line  solutions  for  an 
infinite  span  fixed  wing  model  interacting  with  an  infinite  vortex  filament  at 
variable  planform  intersection  angles  and  axial  displacements.  It  is  then 
applied  to  Che  lifting  line  solution  for  helicopter  rotors.  The  ratio 
obtained  from  Che  fixed  wing  solution  for  comparable  planform  intersection 
angles  and  axial  displacements  in  the  rotor  case  is  used  to  scale  the  local 
vortex  aegment  induced  velocity  influence.  This  model  was  incorporated  into 
the  available  inflow  analysis  and  applied  to  the  above  condition  using  the 
predicted  distorted  wake  geometry.  The  results  of  the  application  of  this 
awdel  to  the  tbtrra  case  had  very  little  effect  on  the  local  airloads.  The 
large  air loading  spike  in  Fig.  101a  for  the  distorted  wake  model  was  not 
significantly  altered.  However,  it  is  recognised  that  the  results  from  this 
limited  study  of  one  rotor  condition  are  not  necessarily  conclusive.  In  order 
to  more  accurately  account  for  the  blade-vortex  interaction  in  the  wake/air¬ 
flow/airloads  analyses  a  finer  numerical  grid  in  the  region  of  the  interaction 
may  be  necessary. 


A  final  model  was  used  based  on  introducing  a  large  vortex  core  radius 
applied  to  the  individual  vortex  segments  influence,  inside  which  the 
influence  is  assumed  to  be  aero.  This  model  could  be  considered  as  a  hypo¬ 
thetical  representation  of  a  burst  vortex.  For  example,  vortex  bursting  was 
previously  simulated  in  Ref.  38.  The  vortex  core  radius  was  assuiaed  to  be 
.03R  (approximately  one-half  chord)  to  ensure  that  the  close  blade/vortex 
interaction  occuring  near  the  90  degree  aximuth  position  at  approximately  .85R 
on  the  blade  (Fig.  96)  would  be  nullified.  It  is  recognised  that  the  vortex 
bursting  concept  as  applied  to  the  rotor  blade-vortex  interaction  has  not  been 
fully  explored  experimentally.  Also,  the  choice  of  vortex  core  radius  in  the 
wake  model  is  clearly  arbitrary  at  this  stage.  The  results  of  the  application 
of  this  model  are  shown  in  Fig.  103  for  the  airloading  and  Fig.  104  for  the 
integrated  blade  loading  which  can  be  compared  %nth  the  results  of  Fig.  101 
and  102.  This  model  improved  the  airload  correlation  by  reducing  the  local 
airloading  on  the  advancing  side  of  the  rotor  disk  without  strongly 
influencing  the  solution  on  the  retreating  side. 

In  summary,  this  brief  limited  scope  study  of  these  BVI  interaction 
models  indicates  the  need  for  continued  activity  tosard  the  development  and 
evaluation  of  a  blade-vortex  interaction  model  which  can  consistently  produce 
accurate  blade  airload  predictions  when  a  blade  is  subject  to  close  vortex 
encounters . 


81 


COMCLODIIIG  REMARKS 


1.  General i»ed  Wake  Model  -  A  firaC  level  general iaed  wake  eodel  was 
developed  for  a  helicopter  rotor  operating  in  steady  level  forward  flight 
based  on  theoretically  predicted  wake  geometries.  The  generalised  wake 
model  consists  of  equations  for  wake  geometry  in  which  tip  vortex  distor¬ 
tions  are  generalised  as  displacements  from  the  classical  undistorted  tip 
vortex  geometry.  The  equations  are  focused  on  the  predominant  distortions 
normal  to  the  rotor  disk.  In  the  generalised  tip  vortex  geometry  equa¬ 
tions,  the  locations  of  points  on  a  tip  vortex  are  defined  in  terms  of  tip 
vortex  age  and  blade  asimuth,  and  the  coefficients  of  the  distortion 
functions  in  the  equations  are  dependent  on  rotor  advance  ratio,  thrust 
coefficient,  and  number  of  blades.  It  was  shown  that  rotor  disk  attitude, 
a  recognised  parameter  for  the  undiatorted  wake  coordinates,  is  not  a 
significant  parameter  for  the  incremental  tip  vortex  distortions.  For 
typical  operating  conditions,  for  this  initial  generalisation,  the  rotor 
operating  condition  parameters  were  emphasised,  and  a  representative 
untapered  blade  design  with  a  linear  twist  of  -10  degrees  and  a  blade 
aspect  ratio  of  16  was  used.  A  limited  study  indicated  that  tip  vortex 
geometry  is  sensitive  to  twist  variations  hut  insensitive  to  aspect  ratio 
variations.  Rotor  wake  geometry  generalisation  for  forward  flight  should 
proceed  beyond  this  initial  first  level  generalisation  to  extend  the  range 
of  operating  and  blade  design  parameters  and  to  refine  the  generalisation 
as  experimental  data  and  improved  wake  methodology,  particularly  for 
blade-vortex  interaction,  become  available. 

2.  Generalised  Wake  Program  Module  -  The  generalised  wake  equations  for 
forward  flight  were  programswd  in  a  generalised  wake  subroutine  which  can 
be  used  to  expediently  generate  wake  coordinates  for  use  with  airflow  and 
airloads  prediction  computer  prograsM.  Numerical  interpolation  procedures 
are  provided  in  this  generalised  wake  module  to  acquire  wake  coordinates 
for  user  selected  combinations  of  advance  ratio,  thrust  coefficient,  disk 
attitude,  and  number  of  blades.  Although  tip  vortex  coordinates  can  be 
generated  for  eny  combination  of  these  parameters,  smphesis  in  the  gener- 
elisatiott  was  placed  on  advance  ratios  between  0.05  and  0.3  and  thrust 
coefficients  between  0.0025  and  0.0075  for  two-  and  four-bladed  rotors. 

The  computer  time  requirement  for  application  of  the  generalised  wake 
module  to  generate  tip  vortex  coordinates  is  insignificant  relative  to  the 
large  requirement  of  a  freely  distorting  wake  geometry  prediction 
analysis. 

3.  Wake  Geometry  Charts  -  A  set  of  wake  geometry  charts  were  developed  which 
provide  a  convenient,  readily  accessible  source  ^or  approximating  rotor 
forward  flight  wake  geometry  and  identifying  potential  locations  of  close 
blade-vortex  passage.  The  cherts  are  presented  as  a  self-contained 


reference  in  Volume  II  of  thie  report.  Charts  for  both  the  classical 
undistorted  wake  model  and  the  generalized  distorted  wake  model  are 
included.  The  charts  provide  the  necessary  infonaation  to  determine  tip 
vortex  coordinates  for  a  range  of  paraswters  sufficient  for  most  conven¬ 
tional  steady  level  flight  operating  conditions.  They  can  be  used  for  a 
variety  of  applications  which  require  the  geometric  relationship  between 
the  tip  vortices  and  spatial  locations  relative  to  the  rotor  blades  and 
other  helicopter  components.  In  particular,  blade/tip  vortex  interaction 
locations  and  rotor  wake  boundaries  can  be  rapidly  defined.  Sample  three- 
views  and  iaosietric  views  of  tip  vortex  patterns  for  two-  and  four-bladed 
rotors  are  also  included. 


Tip  Vortex  Geometry  and  Influence  on  Airloads.  Performance,  and  Bending 
Moments  -  Tip  vortex  geometry  and  the  sensitivity  of  predicted  rotor 
airloads  and  related  rotor  performance  and  blade  bending  moments  to 
various  rotor  inflow  and  wake  geometry  models  were  desMnstrated  for 
forward  flight  conditions.  The  inflow/wake  models  consisted  of  uniform 
rotor  inflow  and  variable  inflow  based  on  classical  undistorted  wake 
geometry,  predicted  distorted  wake  geometry,  and  the  newly  developed 
generalized  wake  geometry.  The  principal  findings  of  this  analytical 
study  are: 

•  Although  sufficient  experiswntal  data  is  not  available  for  a  point-by¬ 
point  comparison  of  measured  and  predicted  forward  flight  wake 
geometry,  the  charactliristic  wake  distortion  features  observed  from 
test  results  are  present  in  the  predicted  and  generalized  tip  vortex 
geometries.  Relative  to  the  classical  undistorted  wake  model,  the 
forward  and  lateral  sides  of  the  wake  are  predicted  to  be  distorted 
toward  the  rotor.  This  results  in  the  prediction  of  close  blade-vortex 
passages  for  typical  low  and  moderate  flight  speed  conditions  which 
introduce  high  harmonic  content  and  severe  local  spanwise  gradients  in 
blade  airloads  which  are  not  predicted  with  the  uniform  inflow  model 
and  the  undistorted  wake  model. 

•  Advance  ratio,  thrust  coefficient,  disk  attitude  and  number  of  blades 
determine  the  classical  undistorted  wake  geometry  in  forward  flight 
based  on  momentum  considerations.  Tip  vortex  distortions  from  the 
undistorted  wake  geometry  were  predicted  to  be  most  sensitive  to 
advance  ratio  variations  and  are  insensitive  to  disk  attitude  varia¬ 
tions.  Unlike  the  hovering  condition,  for  which  the  tip  vortex  coordi¬ 
nates  are  blade  loading  (C^/b)  dependent,  the  forward  flight  wake  is 
predicted  to  be  more  disk  loading  (Cj)  dependent.  The  amplitude  of 
the  tip  vortex  distortions  are  proportional  to  throat  coafficient,  but 
the  characteristic  shape  of  the  principal  distortions  are  insensitive 
to  thrust  coafficient. 


The  influence  of  wake  distortion  on  blade  airloads  is  related  to  the 
degree  of  proximity  of  the  tip  vortices  to  the  rotor  and  the  nuoiber  of 
close  blade-vortex  passages.  For  steady  level  flight,  the  wake 
influence  on  blade  airloads  generally  increases  with  decreasing  advance 
ratio,  decreasing  disk  attitude  and  increasing  number  of  blades. 
Increasing  thrust  coefficient  increases  the  undistorted  wake  deflection 
angle  away  from  the  rotor  but  also  increases  the  distortions  of  the  tip 
vortices  towards  the  rotor,  and  thus  the  influence  of  thrust 
coefficient  on  blade  airloads  is  not  necessarily  as  consistent  as  the 
other  parameters.  The  wake  distortions  and  related  airload  influence 
are  oust  complex  at  low  advance  ratios  (e.g.,  0.05)  where  the 
difficulty  of  accurate  wake  geometry  and  airloads  prediction  is 
compounded  by  the  large  number  of  close  blade-vortex  passages  and  the 
movement  of  the  tip  vortices  above  and  then  down  through  the  rotor  disk 
resulting  in  blade-vortex  impingement. 

The  prediction  of  tip  vortex  distortions  which  are  in  close  proximity 
to  the  blades  at  localised  regions  of  the  rotor  results  in  a  high 
sensitivity  of  predicted  airloads  to  small  differences  in  tip  vortex 
geometry.  For  such  conditions,  differences  were  observed  between  the 
airload  predictions  of  the  generalized  and  predicted  distorted  wake 
models  even  when  the  predicted  wake  geometry  was  well  approximated  by 
the  generalized  wake  model.  The  occurrence  of  close  blade-vortex 
passages,  and  the  high  sensitivity  of  blade  airloads  to  small  devia¬ 
tions  of  blade  station  to  tip  vortex  distance  in  the  airloads  analysis, 
demonstrate  the  requirement  for  an  improved  blade-vortex  interaction 
model  and  further  refinement  of  wake  geometry  prediction  and  general¬ 
ization  techniques.  Fundamental  experimental  and  theoretical  studies 
of  blade-vortex  interaction  involving  vortex  structure  and  blade 
surface  flow  separation  influences  are  required  to  provide  a  physically 
accurate,  and  dependable  blade-vortex  interaction  model  for  incorpora¬ 
tion  in  rotor  wake  and  airloads  smthodology.  Attention  must  also  be 
given  to  the  elimination  of  analytical  modeling  inaccuracies  due  to 
discrete  blade  segment  a»deling  and  discrete  time  steps  in  the  vicinity 
of  close  blade-vortex  passages. 

Regarding  rotor  performance  prediction,  the  classical  undistorted  wake 
calculation  of  rotor  power  for  constant  thrust  was  consistently  higher 
than  the  uniform  inflow  calculation.  Previous  investigations  have 
indicated  that  this  increase  induced  power  results  in  improved  correla¬ 
tion  with  test  data.  The  calculated  power  values  from  the  predicted 
distorted  and  generalised  wake  models  were  scattered  about  the 
classical  wake  values  with  no  well  defined  trend.  The  consideration  of 
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wake  diatortions  is  necessary  for  the  prediction  of  rotor  control 
angles  at  low  speed.  Although  the  need  for  further  analytical  refine¬ 
ment  is  evident,  the  degree  to  idiich  the  distorted  wake  model  was  able 
to  reproduce  the  measured  airload  and  bending  moisent  distributions  is 
encouraging  and  indicative  that,  with  future  emphasis  on  blade-vortex 
interaction  modeling,  distorted  wake  methodology  has  the  potential  to 
provide  an  accurate  predictive  tool  for  blade  airload  and  stress 
analysis. 
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FIGURE  7.  SMOKE  FLOW  VISUALIZATION  OF  ROTOR  WAKE  FOR  HOVER  AND  LOW 
SPEED  FORWARD  FLIGHT.  (REPRODUCED  FR(»1  REF.  6) 
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FIGURE  8.  COMPARISON  OF  MEASURED  ROTOR  WAKE  BOUNDARIES  AND  TIP  VORTEX 
LOCATIONS  FOR  HOVER  AMD  LOW  SPEED  FLIGHT  (REPRODUCED  FROM  REF.  6) 
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FIGURE  11.  REPRESENTATION  OF  BLAK  AND  HAKE  BY  BOUND 
AND  TRAILING  V(»TEX  SEGMENTS 
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nCUU  15.  OOMPAKISON  OF  EXPEtlHENTiO.  AMO  UNDISTORTO)  WAKE  PATTEniS— MODEL  ROTOR 
SIMOLATED  90  KT,  10,000-LB  LIFT  (REPRODUCED  FROM  REF.  5) 


SUPERPOSITION  OF  AXIAL  WAKE  LOCATION  FOR  VARIOUS  BLADE  AZIMUTH  POSITION 
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FIGURE  26.  EFFECT  OF  ROTOR  TIF-FASH  FLAKE  ANCLE  OH  MODEL  ROTOR 
LOADING  AND  WAKE  GEOHEHIY  (REPRODUCED  FROM  REF.  5) 
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riGUKE  28.  EFFECT  OF  THRUST  COEFFICIENT  ON  MODEL  ROTOR  LOADING 
AND  HAKE  GEOMETRY  (RDRODOCED  FROM  REF.  5) 
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FIgura  aa  Variationc  of  tlio  Qonoraliaod  Wako  Envalopo  Function  with  Advance  Rotor 
and  Thruat  Level  for  the  Two  Bladed  Repreeentative  Rotor 
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Figura  31.  Variations  of  tha  Qanarallzad  Waka  Envaiopa  Funotlon  with  Advanea  Rotor 
and  Throat  Laval  for  tha  Four  Bladad  Rapraaantativa  Rotor 
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PHASED  WAKE  AGE,  7w  (DEQ) 

Figur*  32.  Vaitotion  of  tho  Qononillsotf  Wako  Shop#  Funetion  with  Advaneo  Rotor  for 
tha  Two>Bladod  Rapraaantathra  Rotor 
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FICURE  34.  VAKIATION  OF  THE  GENERALIZED  WAKE  SHAPE  FUNCTION  WITH  THRUST 
LEVEL  FOR  THE  TWO-BLADED  REPRESENTATIVE  ROKHl  (M-O.l) 
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FICUU  33.  VABIATIOR  OF  THE  GEMIEALIZED  WAKE  SHAPE  FUNCTION  WITH  THRUST 
LEVEL  FOR  THE  FOUR-ELADED  REPRESENTATIVE  ROTOR  (N*0.1) 
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BLADES  FOR  THE  SAME  THRUST  AND  ADVAIfCE  RATIO  (C^-0. 00321 ,M-0. 05) 
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IZED  CEOHETRIC  SHAPE  PUNCTION 
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FIGURE  39.  COMPARISON  OF  THE  ENVELOPE  FUNCTION  AND  NORMALIZED  GEOMETRIC  SHAPE  FUNCTION  FOR  TWO  ROTOR 
SOLIDITIES  (ASPECT  RATIO)  AT  THE  SAME  CONDITION  (C  -0.0039 ,M=>0. 1)  FOR  THE  TWO  BLADED  ROTOR 
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rZCDIE  41.  COMPARISON  OT  THE  MORNALXZEO  GEOMETRIC  SHAPE  PUMCTXON  FOR  TWO 
DIFFERENT  TWIST  CONDZTXOMS  AT  ONE  FLIGHT  CONDITION  (CfPO.OOSS,  M-0. 1) 


C  SHAPE  FUNCTION  NMHALIZED  CEOHETRIC  SHAPE  FUNCTION 


FieURB  42.  COMPAUSON  Of  THE  MOIMALXZtD  GBOHETMC  WE  FOMCnOH  FOE  1WO 
DIFFEEENT  TWIST  COMOITIOMS  AT  ONE  FUG8T  COMOmOM  (Ct  ■0.032»M^.2) 
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43.  DISTORTED  AND  CLASSICAL  TIP  VORTBC  LONGITUDINAL  COORDINATE  VERSUS  WAKE  ACE 
FOR  VARIOUS  ADVANCE  RATIOS  FOR  THE  REPRESENTATIVE  ROTOR  (Cj-0.0032,  b-2.0) 
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PIGUBE  49a.  GENERALIZED  HAKE  ENVELOPE  FUNCTION  CEARTS  FOR 
FOUR  BLADES  (ft"  .10) 
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riCDRE  51.  PROJECTION  AND  ISOMETRIC  VIEWS  FOR  GENERALIZED  DISTORTED 
TIP  VORTEX,  VARYING  ADVANCE  RATIO  (V/flR  ■  .10) 
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FIGOSE  59.  POTEMTIAL  BLADE/TIP  VCffiTEX  INTERSECTION  PLOT  FOR  TWO  BLADES 
(PCnUAR  FORMAT),  n  -  .30 


AMTIWrUMII 


UMDISTORTEO  WAKE 
WAKE  GENEItALIZATION> 

ACTUAL  CALCULATED  DISTORTED  WAKE 
BLADE  TIP  LOCATIONS 


270  BLADE 


POSITION  OF  IHE  TIP  VORTEX  FROM  HIE  BLADE  AT 


b-2,  M-O.l,  •''b-lBO",  C^.00321,  “TW— 3.0 


POSITION  OP  THE  TIP  VORTEX  PROM  THE  BLADE  AT  -  180 


AXIAL  POSITION  Of  THE  TIP  V(»tTEX  PROH  THE  BLADE  AT  270 


- FUDICTED  DISTORTED  UMCE 

- CUSSICAL  AM)  CEHERALIZED  UARB 
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riCUIE  70b.  COMPABISON  OP  PREDICTED  SPANWISE  LIFT  DISTRIBUTION  FOR  VARIOUS  INFLOW  MODELS  AT  -  90 
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FIGURE  75.  POTENTIAL  BLADE/TIP  VORTEX  INTERSECTION  PLOT  FOR  FOUR  BLADES 
(POLAR  FORMAT).  It  -  .05 
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16  Abstrar^An  analyticif' Invest igation  to  generalise  »ake  geometry  of  a  helicopter 
rotor  in  steady  level  forward  flight  and  to  demonstrate  the  influence  of  i^hc  deforma¬ 
tion  in  the  ptediction  of  rotor  airloads  and  performance  is  described.  YeJ&rrt  i. 

—  In-j^olumd'i^  a  first  level  generalised  wake  model  ^e  pfeeewte^eliieli  ir~based 
on  theoretically  predicted  tip  vortex  geometries  for  a  selected  representative  blade 
design.  The  tip  vortex  distortions  are  generalised  in  equation  form  as  displacegwnts 
from  the  classical  undistorted  tip  vortex  geometry  in  terms  of  vortex  age,  blade 
asimuth,  rotor  advance  ratio,  thrust  coefficient,  and  number  of  blades.  These 
equations  were  programmed^  e  eempmtee  medwi^to  provide  distorted  wake  coordinates 
at  very  low  cost  for  use  in  rotor  airflow  and  airloads  prediction  analysea.  The 
sensitivity  of  predicted  rotor  airloads,  performance,  and  blade  bending  moments  to  the 
modeling  of  the  tip  vortex  distortion  are  demonstrated  for  low  to  moderately  high 
advance  ratios  for  a  representative  rotor  and  the  H-34  rotor.  Comparisons  with  H-34 
rotor  test  data  demonstrate  the  effects  of  the  classical,  predicted  distorted,  and  the 
newly  developed  generalised  wake  models  on  airloads  and  blade  bending  moments.  *Wie ***  — 
use  of  distorted  wake  models  results  in  the  occurrence  of  numerous  blade-vortex  inter- 
Tetions  on  the  forward  and  lateral  sides  of  the  rotor  disk.  The  significance  of  these 
interactions  is  related  to  the  number  and  degree  of  proximity  to  the  blades  of  the  tip 
vortices.  The  correlation  obtained  with  the  distorted  wake  models  (generalised  end 
predicted)  is  encoursging.j^wever,  the  reaulting  high  sensitivity  of  the  predicted 
aitloads  to  small  deviationk  in  tip  vortex  position  deaKMiatrate  the  requirement  for 
improved  blade-vortex  interanion  modeling. 

A  sec  of  wake  geometry  charts  are  presented  in  Voliac  Zl  to  provide  a 
convenient,  readily  accessible  source  for  approximating  rotor  forward  flight  wake 
geometry  and  identifying  wake  boundaries  and  locations  of  blade-vorten  passage. 
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